Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,732,141 G→C 6.2% intergenic (+18/‑60) xylF → / → xylG xylose ABC transporter periplasmic binding protein/xylose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,732,1410GC 6.2% 323.2 / 8.2 96intergenic (+18/‑60)xylF/xylGxylose ABC transporter periplasmic binding protein/xylose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (39/51);  new base C (2/4);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

AAAGCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTG  >  NC_000913/3731992‑3732288
                                                                                                                                                     |                                                                                                                                                   
aaaGCAGATACCACACTGAATAATGGCCGGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGcgcg                                                                                                                                                     >  2:470882/1‑150 (MQ=255)
 aaGCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCgg                                                                                                                                                    <  1:867752/150‑1 (MQ=255)
   gCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAg                                                                                                                                                  >  2:843192/1‑150 (MQ=255)
   gCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAg                                                                                                                                                  >  2:1222326/1‑150 (MQ=255)
      gATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCgg                                                                                                                                               >  2:351363/1‑150 (MQ=255)
        tACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCgggg                                                                                                                                             <  1:1190801/150‑1 (MQ=255)
          ccACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGc                                                                                                                                           >  1:1187005/1‑150 (MQ=255)
           caccCTCAATAAGCGCCTGAAAGATGCCCCCTCCCGCCTCCTGCCACCGATCGAGGTGAATCAAAACAACATCAAAGATCCGGTAATTAAAGACGTATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCg                                                                                                                                          <  2:310468/150‑1 (MQ=255)
            acacTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGt                                                                                                                                         >  1:982915/1‑150 (MQ=255)
             cacTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTg                                                                                                                                        <  1:726267/150‑1 (MQ=255)
             cacTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTg                                                                                                                                        <  2:581322/150‑1 (MQ=255)
               cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTTAATAAAAACAACATCAAAGATACGGTAATTAAAGACGTATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAt                                                                                                                                      <  2:1335395/150‑1 (MQ=255)
               cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGtaac                                                                                                                                      <  1:1226818/150‑4 (MQ=255)
               cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAt                                                                                                                                      >  1:631206/1‑150 (MQ=255)
               cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAt                                                                                                                                      >  2:520041/1‑150 (MQ=255)
                 gAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTCACGCCCCCGCTCCGCGCTGGGGCGtaacgc                                                                                                                                    <  2:656681/150‑6 (MQ=255)
                  aataatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtc                                                                                                                                   <  2:452878/150‑1 (MQ=255)
                  aataatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtc                                                                                                                                   <  1:887534/150‑1 (MQ=255)
                    taatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctc                                                                                                                                 >  1:386122/1‑150 (MQ=255)
                    taatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctc                                                                                                                                 <  1:1093336/150‑1 (MQ=255)
                      atGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctctc                                                                                                                               <  1:1312208/150‑1 (MQ=255)
                      atGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctctc                                                                                                                               <  2:737101/150‑1 (MQ=255)
                            tGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGcc                                                                                                                         >  1:249754/1‑150 (MQ=255)
                               aaGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCcgcg                                                                                                                      <  1:984874/150‑1 (MQ=255)
                               aaGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCcgcg                                                                                                                      <  1:843192/150‑1 (MQ=255)
                                    gTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaat                                                                                                                 <  1:925008/150‑1 (MQ=255)
                                    gTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaat                                                                                                                 <  1:1390616/150‑1 (MQ=255)
                                    gTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaat                                                                                                                 <  1:1385758/150‑1 (MQ=255)
                                      ccccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatga                                                                                                               >  1:463192/1‑150 (MQ=255)
                                       cccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGGGCTTTCTCTCCATGCCGCGtgaatgaa                                                                                                              >  2:1244321/1‑150 (MQ=255)
                                       cccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGAAACGGTAATTAAAGAAGGATTCCACAAAGAGAGCGCGCTGTAAGCGTTACGCCCCACCGCGGAGCGGGGGCGTCATTTCTCTCCATGCCGCGtgaatgaa                                                                                                              >  2:461626/1‑150 (MQ=255)
                                           cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTgg                                                                                                          <  2:839813/150‑1 (MQ=255)
                                           cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTCACGCCCCCGCTCCGCGCTGGGGCGtaacgct                                                                                                                                   <  2:595902/125‑7 (MQ=38)
                                           cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTCACGCCCCCGCTCCGCGCTGGGGCGtaacgct                                                                                                                                   >  1:595902/1‑119 (MQ=38)
                                            ccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGc                                                                                                         >  1:1315757/1‑150 (MQ=255)
                                                   ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                  >  1:490119/1‑150 (MQ=255)
                                                   ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                  <  2:492284/150‑1 (MQ=255)
                                                   ctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTg                                                                                                  <  1:1324407/150‑1 (MQ=255)
                                                    tGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                 <  2:761085/150‑1 (MQ=255)
                                                    tGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                 >  2:396576/1‑150 (MQ=255)
                                                    tGACACCGATCGATGTGAATAAAAACAACATAAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                 >  2:1052745/1‑150 (MQ=255)
                                                      acacCGATCGATGTTAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAg                                                                                               <  2:1315757/150‑1 (MQ=255)
                                                         ccGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCg                                                                                            >  2:1232594/1‑150 (MQ=255)
                                                         ccGATCGATGTGAATAAAAACAAAATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGAGGGGGCGTGATTTCTCTCCATGCAGAGTGAATGAATTGGCTTAGGTGGAGTCg                                                                                            >  1:33726/1‑150 (MQ=255)
                                                          ccagcgatGTGAAAAACGACAACATCAAAAATAGGGTAATTAAAGACGGATTCCACAAAGAGAGAGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGGCGCGTGAATGAATTGGCTTAGGTGGAGTCGt                                                                                           >  2:121968/5‑150 (MQ=255)
                                                           gatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGtt                                                                                          >  1:329231/1‑150 (MQ=255)
                                                           gatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGtt                                                                                          <  2:1091276/150‑1 (MQ=255)
                                                                atGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                     <  2:1289700/150‑1 (MQ=255)
                                                                  gtgAATAAAAACAACATCAAAGATACGGTAATTAAAGCCGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCtt                                                                                   >  2:257665/1‑150 (MQ=255)
                                                                         aaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACt                                                                            <  1:879655/150‑1 (MQ=255)
                                                                         aaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACt                                                                            <  2:631206/150‑1 (MQ=255)
                                                                          aaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACtt                                                                           <  2:751819/150‑1 (MQ=255)
                                                                           aacaacATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          <  2:643867/150‑1 (MQ=255)
                                                                           aacaacATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          <  1:1362147/150‑1 (MQ=255)
                                                                           aacaacATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          >  1:1380797/1‑150 (MQ=255)
                                                                           aacaacATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                          <  2:463192/150‑1 (MQ=255)
                                                                             caacaTCAAAGATACGGTAATTAAAGCCGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGaa                                                                        <  2:813121/150‑1 (MQ=255)
                                                                              aacaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGaaa                                                                       <  2:1187005/150‑1 (MQ=255)
                                                                                caTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATg                                                                     >  1:722338/1‑150 (MQ=255)
                                                                                 aTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATga                                                                    >  2:1335439/1‑150 (MQ=255)
                                                                                   cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaag                                                                  >  2:1197096/1‑150 (MQ=255)
                                                                                   cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaag                                                                  <  1:278035/150‑1 (MQ=255)
                                                                                   cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaag                                                                  <  2:531819/150‑1 (MQ=255)
                                                                                      aGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAAc                                                               <  2:982915/150‑1 (MQ=255)
                                                                                      aGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAAc                                                               <  1:1052745/150‑1 (MQ=255)
                                                                                       gATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGa                                                                         <  1:1034938/139‑1 (MQ=255)
                                                                                       gATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGa                                                                         >  2:1034938/1‑139 (MQ=255)
                                                                                           cGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAc                                                          <  1:1169/150‑1 (MQ=255)
                                                                                           cGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAc                                                          <  1:1348824/150‑1 (MQ=255)
                                                                                                aTTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCaaaa                                                     <  2:722338/150‑1 (MQ=255)
                                                                                                aTTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATAAAGAACATTACCaaaa                                                     >  1:687252/1‑150 (MQ=255)
                                                                                                      gACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCg                                               <  1:172246/150‑1 (MQ=255)
                                                                                                        cGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGc                                             <  2:418248/150‑1 (MQ=255)
                                                                                                             tCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgtgt                                        <  2:401846/150‑1 (MQ=255)
                                                                                                              ccACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGtgtg                                       <  1:470882/150‑1 (MQ=255)
                                                                                                              ccACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGtgtg                                       >  1:184602/1‑150 (MQ=255)
                                                                                                              ccACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGtgtg                                       >  2:595334/1‑150 (MQ=255)
                                                                                                              ccACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGtgtg                                       >  1:382741/1‑150 (MQ=255)
                                                                                                              ccACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGtgtg                                       >  1:872239/1‑150 (MQ=255)
                                                                                                               cacaAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAAAATTACCAAAACCTTCCGCAGTGTGa                                      >  1:573631/1‑150 (MQ=255)
                                                                                                                acaAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGag                                     <  1:648032/150‑2 (MQ=255)
                                                                                                                   aatatagcgagtTTTAAGCGTTACGGCCCAGCGCGGAGCGGGGGCGTGATTTCTATCCATGCCGCGTGAATGAATTGGGTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGc                                  <  1:619695/145‑1 (MQ=255)
                                                                                                                       gagcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGAtt                              >  2:1018984/1‑150 (MQ=255)
                                                                                                                        agcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTg                             <  1:257665/150‑1 (MQ=255)
                                                                                                                         gcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTTAAGGCGATTGa                            <  2:413245/150‑1 (MQ=255)
                                                                                                                           gagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATa                          <  2:490119/150‑1 (MQ=255)
                                                                                                                           gagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATa                          <  1:595334/150‑1 (MQ=255)
                                                                                                                              cTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACg                       >  2:1382305/1‑150 (MQ=255)
                                                                                                                               tGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTg                             <  1:211323/143‑1 (MQ=255)
                                                                                                                               tGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTg                             >  2:211323/1‑143 (MQ=255)
                                                                                                                                     aatcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTATCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCtt                >  1:1291272/5‑150 (MQ=255)
                                                                                                                                      gTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg               <  1:520041/150‑1 (MQ=255)
                                                                                                                                      gTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg               <  1:734390/150‑1 (MQ=255)
                                                                                                                                       ttACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGc              >  1:1427786/1‑150 (MQ=255)
                                                                                                                                       ttACGCCCCAGCCCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGATTGAATTGGCTTAGGTGGAGTAGTTATCCCTAATCTACTTGAAATGAAGAACATAACCAAAACCTTCGCCAGTGTGAAGGCGATTGAATACGTCTGCTTGc              >  2:740843/1‑150 (MQ=255)
                                                                                                                                          cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt           >  1:922109/1‑150 (MQ=255)
                                                                                                                                          cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt           <  1:1222326/150‑1 (MQ=255)
                                                                                                                                               cAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAAt      >  1:329210/1‑150 (MQ=255)
                                                                                                                                                   gcgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCGACTGGAAATGAAGAACATTACCAAACCCATCGGCAGTGTGAAGGCGATTGCTAACGTCTGCTTGCGGTTGAATGCTg  >  1:1384914/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAAGCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTG  >  NC_000913/3731992‑3732288

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑