Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 66,465 T→C 5.6% R14R (AGA→AGG)  USA300HOU_0062 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073066,4650TC 5.6% 235.7 / 3.7 90R14R (AGA→AGG) USA300HOU_0062hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base C (2/3);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

CATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCCTTCTT  >  CP000730/66326‑66603
                                                                                                                                           |                                                                                                                                          
cATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTACTCTGTATTAAGTTGATTGTCTCGTTCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCctc                                                                                                                                           <  2:204660/141‑1 (MQ=255)
cATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCctc                                                                                                                                           <  1:97215/141‑1 (MQ=255)
cATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCctc                                                                                                                                           <  1:134714/141‑1 (MQ=255)
   gtgtAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTtg                                                                                                                                        >  2:739833/1‑141 (MQ=255)
     gtAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCCCTtgtg                                                                                                                                      >  2:118824/1‑141 (MQ=255)
     gtAATTGACAACATATTATCATCTTCATGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCCCTtgtg                                                                                                                                      >  2:502734/1‑141 (MQ=255)
     gtAATTGACAACATATTATCATCTTATTGTAGTTTGAAAATGTGTAGGTATTCTGTATTACGTTGATTGTCTCGATCTTTGACCATGTTCCTATGTAAGCTTTGACGGTCTGGAGCTAGTTGTTCACTCATGCccctactg                                                                                                                                      >  2:300843/1‑137 (MQ=255)
          tGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTa                                                                                                                                 >  2:132333/1‑141 (MQ=255)
          tGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAAACTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTa                                                                                                                                 >  1:555160/1‑141 (MQ=255)
            acaacaTATTATCATCTTCTGGTAGTTTGAAAATGTGTAGGTAATCTGTCTTAAGTTGATTGTCTCGATCTTTTACCATGTTGCAGAGTAAGATTTAAAGGTCTGGAGATAGTTGTTAACTAATGGATCTtgtagtgtttc                                                                                                                               >  1:507087/1‑133 (MQ=37)
               acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCgat                                                                                                                            <  2:557612/141‑1 (MQ=255)
               acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCgat                                                                                                                            >  1:435105/1‑141 (MQ=255)
               acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCgat                                                                                                                            >  2:196754/1‑141 (MQ=255)
               acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTAGTGATATATCgat                                                                                                                            >  1:21458/1‑141 (MQ=255)
               acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCACGATCTTTTACCATGTTCCACAGTAAGATTTAAAGGTCTAGAGATCTGTCTTCCCTAATCCCTCTTGTGACGTAtggct                                                                                                                            >  2:313871/1‑137 (MQ=25)
                      atcatcTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGAtttt                                                                                                                     >  2:718740/1‑141 (MQ=255)
                         atcTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCAt                                                                                                                  <  1:95134/141‑1 (MQ=255)
                               ttGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCAAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTAtt                                                                                                            >  2:619410/1‑141 (MQ=255)
                                 gTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATGGAATGATTTTCATACTAtttt                                                                                                          >  2:125246/1‑141 (MQ=255)
                                 gTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTAtttt                                                                                                          >  2:547661/1‑141 (MQ=255)
                                 gTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACCAATGCCGCTTGTGATGTATCGATTGATTTTCATACTAtttt                                                                                                          >  2:479612/1‑141 (MQ=255)
                                   aGTTTGAAAATGTTTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTATCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTTTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTtcc                                                                                                        <  2:226215/141‑1 (MQ=255)
                                   aGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTtcc                                                                                                        >  2:606506/1‑141 (MQ=255)
                                         aaaaTGTGTAGGCAATCTGGATTAAGTTGATTGTCTCGATCTTTTTCCATGTTCCAAAGTAATATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCAtt                                                                                                  <  2:455558/141‑1 (MQ=255)
                                             tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATTTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  2:689063/141‑1 (MQ=255)
                                             tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  2:22067/141‑1 (MQ=255)
                                             tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  1:647794/141‑1 (MQ=255)
                                             tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  1:749569/141‑1 (MQ=255)
                                             tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              >  1:91211/1‑141 (MQ=255)
                                             tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  1:180006/141‑1 (MQ=255)
                                              gtgtAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGCTTGATTTTAATACTATTTTCGTCCCTTTTGCt                                                                                             >  1:267010/1‑141 (MQ=255)
                                              gtgtAGGTAATATGTATTAAGTTGATTGTCTCGATCTTTTAACATGTTCCAAAGTAAGATTTGAAGGTCTAGCGATAGTTGTTCACTAATGCCTCTTGTGATGTATGGATTGATTTACTTAATGTCTTCCTCCATTTTgat                                                                                             >  2:37737/1‑139 (MQ=37)
                                                 tAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCtttt                                                                                          <  1:472388/141‑1 (MQ=255)
                                                 tAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCtttt                                                                                          <  1:727097/141‑1 (MQ=255)
                                                   ggTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTtttct                                                                                        >  2:291286/1‑141 (MQ=255)
                                                        tCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCat                                                                                   <  2:82844/141‑1 (MQ=255)
                                                         cTGTCTTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTTAGATTTGAATGTTTATTTAGAGTTTTTCACTAATGCCTCTTGGTATGTATTGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCatg                                                                                  <  2:369174/141‑1 (MQ=255)
                                                          tGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCatga                                                                                 >  1:668213/1‑141 (MQ=255)
                                                          tGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATGGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTTatga                                                                                 >  2:621178/1‑141 (MQ=255)
                                                            tATTAAGTTGATTGTCCCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGGTGTATCGATGGGTTTTCTTACTATTTTCCTCAAGTTTGCTTTTCTTTCatgatg                                                                               >  1:303945/1‑141 (MQ=37)
                                                              ttAAGTTGATTGTCTCGATCTTTTCCCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTc                                                                             <  1:492717/141‑1 (MQ=255)
                                                              ttAAGTTGATTGTCTCGATCCTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTc                                                                             <  1:323847/141‑1 (MQ=255)
                                                                aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATCCTATTTTCCTCCATTTTGTTTCTCTTTCATGATGTCaa                                                                           <  1:452142/141‑1 (MQ=255)
                                                                aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCaa                                                                           <  1:629682/141‑1 (MQ=255)
                                                                aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCaa                                                                           <  1:251508/141‑1 (MQ=255)
                                                                aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCaa                                                                           <  1:312320/141‑1 (MQ=255)
                                                                 aGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAAt                                                                          <  1:36629/141‑1 (MQ=255)
                                                                       ttGTCTCGATCTTTTACCATGTTCCAAATTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTc                                                                    >  1:512617/1‑141 (MQ=255)
                                                                       ttGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTc                                                                    >  1:344528/1‑141 (MQ=255)
                                                                          tctcGATCTTCTACCCTGTTCCAAAGTAAGATTTGAAGGCCCAGAGATAGTTGTTCACTAATGCCTTTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGtt                                                                 <  2:3103/141‑1 (MQ=25)
                                                                            tcGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTaa                                                               <  1:739833/141‑1 (MQ=255)
                                                                             cGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAAt                                                              >  2:418898/1‑141 (MQ=255)
                                                                                tCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATgct                                                           >  2:32402/1‑139 (MQ=255)
                                                                                tCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATACCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCCCTTCGTTAATgcc                                                           >  2:319409/1‑139 (MQ=255)
                                                                                    ttACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGCCTGTa                                                       >  2:314580/1‑141 (MQ=255)
                                                                                     tatcATGTTCCAAAGTAAGATTTGAGGGTCTATAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTaa                                                      <  1:784975/138‑1 (MQ=255)
                                                                                     tACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTCATGACTGTaa                                                      >  1:552731/1‑141 (MQ=255)
                                                                                     tACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTaa                                                      <  1:287156/141‑1 (MQ=255)
                                                                                     tACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTaa                                                      <  1:210896/141‑1 (MQ=255)
                                                                                      tccATGTTCCAAAGTAAGCTTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAAc                                                     <  1:161249/140‑1 (MQ=255)
                                                                                       ccATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACa                                                    <  2:707425/141‑1 (MQ=255)
                                                                                           gTTCCAAAGTAGGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGatt                                                >  2:95137/1‑140 (MQ=255)
                                                                                           gTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGatt                                                >  1:233273/1‑140 (MQ=255)
                                                                                             tCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATAtt                                              >  1:466535/1‑141 (MQ=255)
                                                                                              ccAAAGTAAGATTTTAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGacagtt                                             <  2:559517/141‑6 (MQ=255)
                                                                                              ccAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATAttt                                             >  1:717923/1‑141 (MQ=255)
                                                                                               ccaaGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTtg                                            <  2:377693/139‑1 (MQ=255)
                                                                                               cAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTtg                                            <  2:377669/141‑1 (MQ=255)
                                                                                                  atatagcttttAAGGTCTAGCGATAGTTGTTCACTAATGCCTCTTGTGGTTTTTCGATTGATTTCCATCCTTTTTTCATCCATTTTTCTTTTCTCTTATAATTTCAATAACTTAGTCTATGACTGGATCAGATATTTGTGc                                         <  2:273788/134‑1 (MQ=255)
                                                                                                  aGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTATTTCATGATGTCAATCACTTCGTTAATGACTGTAACCGATATTTGTGc                                         >  1:421806/1‑141 (MQ=255)
                                                                                                    tAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCa                                       <  2:36895/141‑1 (MQ=255)
                                                                                                     aaGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATAACTTCGTTAATGACTGTAACAGATATTTTTGcccc                                      >  2:266008/1‑139 (MQ=255)
                                                                                                       gATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTATTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACGTCGTTCATGACTGTAACAGATTTTTTGACCACtt                                    >  2:425040/1‑141 (MQ=25)
                                                                                                           tGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCGATCACTTCGTTAATGACTGTAACAGATATTTGTGCAACTTTGAt                                >  1:150832/1‑141 (MQ=255)
                                                                                                           tGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGAt                                >  1:762738/1‑141 (MQ=255)
                                                                                                             aaGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATcc                              <  1:354420/141‑1 (MQ=255)
                                                                                                               ggTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCaa                            <  2:3412/141‑1 (MQ=255)
                                                                                                                gtttAGCGATATGTGTTCACTAACGCCTTTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           <  2:555160/138‑1 (MQ=37)
                                                                                                                gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           <  2:495259/141‑1 (MQ=255)
                                                                                                                gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           <  1:132333/141‑1 (MQ=255)
                                                                                                                gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           >  2:31069/1‑141 (MQ=255)
                                                                                                                gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAAAGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           >  2:360957/1‑141 (MQ=255)
                                                                                                                  cTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGCCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAatta                         <  2:605857/141‑1 (MQ=255)
                                                                                                                          aGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTAGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATgat                 >  2:464241/1‑139 (MQ=255)
                                                                                                                          aGTTGTTCACAAATGCCGCTTGTGATGTATCGATGGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTCACAGATATTTGTGCCACTTTGATCCAATTATTCATGGt                 >  2:275144/1‑141 (MQ=255)
                                                                                                                             tgttCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTCGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTAt              >  1:298654/1‑141 (MQ=255)
                                                                                                                              gttCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTGTTTTGCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTAtt             >  2:608161/1‑141 (MQ=255)
                                                                                                                                   cTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCCTGGTTATTCcttc        >  2:351518/1‑141 (MQ=255)
                                                                                                                                    tAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcc       <  2:346501/141‑1 (MQ=255)
                                                                                                                                     aaTGCCCCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcct      <  1:502734/141‑1 (MQ=255)
                                                                                                                                     aaTGCCCCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcct      <  1:118824/141‑1 (MQ=255)
                                                                                                                                     aaTGCCCCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcct      <  1:300843/141‑1 (MQ=255)
                                                                                                                                        gCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTGATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCcttct   >  1:521960/1‑141 (MQ=255)
                                                                                                                                         ccTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCcttctt  <  2:435105/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCCTTCTT  >  CP000730/66326‑66603

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑