Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,732,141 G→C 5.0% intergenic (+18/‑60) xylF → / → xylG D‑xylose transporter subunit/D‑xylose ABC transporter dual domain ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,732,1410GC 5.0% 426.8 / 11.5 139intergenic (+18/‑60)xylF/xylGD‑xylose transporter subunit/D‑xylose ABC transporter dual domain ATPase
Reads supporting (aligned to +/- strand):  ref base G (70/62);  new base C (3/4);  total (73/66)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-02

Only 100 of 145 total aligned reads displayed.

GCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTGG  >  NC_000913/3731995‑3732289
                                                                                                                                                  |                                                                                                                                                    
gCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAg                                                                                                                                                   >  1:1344066/1‑150 (MQ=255)
  aGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTATGCGTTACGCCCCAGCGCGGAGCg                                                                                                                                                 >  1:1098380/1‑150 (MQ=255)
    aTACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCggg                                                                                                                                               >  2:450334/1‑150 (MQ=255)
     tACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTTAGCGTTACGCCCCAGCGCGGAGCgggg                                                                                                                                              >  1:624972/1‑150 (MQ=255)
     tACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGATGTAAGCGTTACGCCCCAGCGCGGAGCgggg                                                                                                                                              >  2:1429711/1‑150 (MQ=255)
        cacaCTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCg                                                                                                                                           >  1:431520/1‑150 (MQ=255)
         acacTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAATAGAGCGAGCTGTAAGCGATACGCCCCACCGCGGAGCGGGGGCGt                                                                                                                                          >  1:411883/1‑150 (MQ=255)
           acTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTgt                                                                                                                                        >  1:631055/1‑149 (MQ=255)
           acTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTCACGCCCCCGCTCCGCGCTGGGGCGtaa                                                                                                                                        <  2:370253/150‑3 (MQ=255)
            cTGAATAATGGCCTGAACGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTCAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGtaac                                                                                                                                       <  1:489396/150‑4 (MQ=255)
            cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCGtaac                                                                                                                                       <  2:1362670/150‑4 (MQ=255)
              gAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAttt                                                                                                                                     >  2:728160/1‑150 (MQ=255)
              caataatGGCCTGACAGATGTCCCCACCCGCCTCCTGACACCGCTCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAACGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAttt                                                                                                                                     <  2:518342/149‑1 (MQ=255)
                 taatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctc                                                                                                                                  <  2:1344066/150‑1 (MQ=255)
                  aatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGACTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctct                                                                                                                                 <  1:701234/150‑1 (MQ=255)
                    tGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTcc                                                                                                                               <  2:1399985/150‑1 (MQ=255)
                    tGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTcc                                                                                                                               <  2:878219/150‑1 (MQ=255)
                     ggCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGCTTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCa                                                                                                                              >  1:1058659/1‑150 (MQ=255)
                     ggCCTGAAAGATGTCCCCTCCCCCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGCGCGAGCTGTAATCGTAACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCa                                                                                                                              >  2:1651528/1‑150 (MQ=255)
                      gCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCAt                                                                                                                             >  2:740555/1‑150 (MQ=255)
                          taaaGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCCCAAAGAGCGCGAGCTGTAAGCGTTCCGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCcg                                                                                                                         <  2:214679/149‑1 (MQ=255)
                          gAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGCTTCCACAAAGAGAGCGAGCTCTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGcag                                                                                                                         >  1:279683/1‑148 (MQ=255)
                          gAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCcg                                                                                                                         <  2:1280852/150‑1 (MQ=255)
                            aaGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGCGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCcgcg                                                                                                                       >  2:1369206/1‑150 (MQ=255)
                             aGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGAGAGCTGTAAGCGTTACGCCCAAGCGCGGAGAGGGGGCGTGATTTCTCTCCATGCCGCGt                                                                                                                      >  1:737888/1‑150 (MQ=255)
                              gATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtg                                                                                                                     >  2:1688935/1‑150 (MQ=255)
                                tGTCCCCTCCCGCCTCCTGACACCGATCGATTTGAATAAAAAGAACATCAACGATACGGTAATTAAAGCCGGATTCCAGAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaa                                                                                                                   <  1:1575173/150‑1 (MQ=255)
                                 gTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaat                                                                                                                  <  2:1320309/150‑1 (MQ=255)
                                 gTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaat                                                                                                                  <  1:538666/150‑1 (MQ=255)
                                  tCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTGCTCTCCATGCCGCGtgaatg                                                                                                                 >  1:301922/1‑150 (MQ=255)
                                  tCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATAAAAGACGGAATCCACAAAGAGAGCGAGCTGAAAGCGTTACGCCCCAGCGCGCAGCGGGGCCGTGATTTCTCTCAAAGCCGCGtgaatg                                                                                                                 >  2:544604/1‑150 (MQ=255)
                                   ccccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatga                                                                                                                >  2:1006906/1‑150 (MQ=255)
                                      cTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAAtt                                                                                                             <  1:358415/150‑1 (MQ=255)
                                         ccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGc                                                                                                          >  2:1153219/1‑150 (MQ=255)
                                          cGCCCCCTGATACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCACGTGAATGAATTGGCt                                                                                                         <  2:45631/150‑1 (MQ=255)
                                             ctcctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCAATACAGCGTGAATGAATTGGCTTAg                                                                                                      >  1:634651/1‑150 (MQ=255)
                                                 tGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAGAGAGAGCGCGCTGTAAGCGTTACGCCCCAGAGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                  >  2:174399/1‑150 (MQ=255)
                                                 tGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                  <  1:1724838/150‑1 (MQ=255)
                                                 cgACACCGATCGATGTGAATAACAACAACCTCAAAGATACGGTAATTAAAGACGGATACCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTgg                                                                                                  <  2:323341/149‑1 (MQ=255)
                                                     acCGATCGATGTGAATAAAAACAACATAAAAGATACGGTAATTAACGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGTGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTc                                                                                              <  2:1458556/150‑1 (MQ=255)
                                                      ccGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCg                                                                                             >  2:365346/1‑150 (MQ=255)
                                                        gatcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGtt                                                                                           <  2:1384609/150‑1 (MQ=255)
                                                            gatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGc                                                                                       >  1:381208/1‑150 (MQ=255)
                                                            gatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCGTAGGTGGAGTCGTTATGc                                                                                       >  1:1564242/1‑150 (MQ=255)
                                                             ttgtgAATAAACACCACATCAACGATACGGTAATTACAGACGGCTTCCACACCGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                      <  1:1471170/149‑1 (MQ=255)
                                                             atGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGGAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                      <  1:1114521/150‑1 (MQ=255)
                                                              tgtgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCt                                                                                     <  1:945446/150‑1 (MQ=255)
                                                              tgtgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCt                                                                                     >  1:1654271/1‑150 (MQ=255)
                                                                 gAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGTCGCGTGAATGAATTGGGTTAGGTGGAGTCGTTATGCCtttt                                                                                  >  2:482793/1‑148 (MQ=255)
                                                                   aTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCt                                                                                >  2:1246329/1‑150 (MQ=255)
                                                                    aaaaaaCACCATCAAAGATACCGTAATTAAAGACGGATTCCACACAGAGACCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTCGAGTCGTTATGCCTTATCTa                                                                               <  1:1272068/149‑1 (MQ=255)
                                                                     aaaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTAc                                                                              >  2:1011648/1‑150 (MQ=255)
                                                                      aaaaCATCATCAAAGTTACGGTAATTAAAGTCGGTTTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACt                                                                             <  2:411883/150‑1 (MQ=255)
                                                                        aacaacATCAAAGATACGGTAATTACAGACGGCTTCCCCAAAGAGAGCGTGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTg                                                                           <  2:1333342/150‑1 (MQ=255)
                                                                              aTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATga                                                                     >  2:1040120/1‑150 (MQ=255)
                                                                              aTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATga                                                                     <  2:750352/150‑1 (MQ=255)
                                                                              aTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTAATTTAAATga                                                                     >  2:1422406/1‑150 (MQ=255)
                                                                                cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGATATgaag                                                                   >  2:265326/1‑150 (MQ=255)
                                                                                cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaag                                                                   >  1:248686/1‑150 (MQ=255)
                                                                                 aaaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                  <  1:1703749/150‑1 (MQ=255)
                                                                                    gATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACa                                                               <  2:475444/150‑1 (MQ=255)
                                                                                    gATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACa                                                               <  1:1135158/150‑1 (MQ=255)
                                                                                      tACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACAtt                                                             >  2:1498799/1‑150 (MQ=255)
                                                                                          gTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCa                                                         <  1:1490161/150‑1 (MQ=255)
                                                                                          gTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCa                                                         >  1:1021658/1‑150 (MQ=255)
                                                                                             aTTAAAGACGTATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCaaaa                                                      <  2:937529/150‑1 (MQ=255)
                                                                                             aTTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCaaaa                                                      >  1:350026/1‑150 (MQ=255)
                                                                                              ttAATGACGGAGTCTACAACGACAGCGCGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGTCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGGCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAAc                                                     <  2:281192/150‑1 (MQ=255)
                                                                                              ttAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATTAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAAc                                                     <  2:138992/150‑1 (MQ=255)
                                                                                                   gACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGTAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCg                                                <  1:1051430/150‑1 (MQ=255)
                                                                                                       gATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAg                                            <  1:700093/150‑1 (MQ=255)
                                                                                                       gATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAg                                            >  1:55879/1‑150 (MQ=255)
                                                                                                       gATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTCGCTTAGGTGGAGTCGTTATGCCTTATCTCCTTGAAATGAAGCACATTACCAACACCTTCGGCAg                                            >  1:55267/1‑150 (MQ=255)
                                                                                                        aTTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgt                                           >  1:1053437/1‑150 (MQ=255)
                                                                                                          tCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgtgt                                         <  1:686531/150‑1 (MQ=255)
                                                                                                           ccACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGtgtg                                        >  1:1450970/1‑150 (MQ=255)
                                                                                                               aaaGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAgg                                    <  1:1611458/150‑1 (MQ=255)
                                                                                                                    gagcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGCtt                               >  1:250842/1‑150 (MQ=255)
                                                                                                                    gagcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGAtt                               <  1:1731790/150‑1 (MQ=255)
                                                                                                                     agcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTg                              >  2:507401/1‑150 (MQ=255)
                                                                                                                         agcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATaa                          >  1:1049719/1‑150 (MQ=255)
                                                                                                                          gcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAAc                         >  2:1330502/1‑150 (MQ=255)
                                                                                                                          gcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAAc                         <  1:740555/150‑1 (MQ=255)
                                                                                                                             gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                      >  2:1193652/1‑150 (MQ=255)
                                                                                                                             gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                      >  2:716711/1‑150 (MQ=255)
                                                                                                                             gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                      <  2:431520/150‑1 (MQ=255)
                                                                                                                             gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                      <  1:450334/150‑1 (MQ=255)
                                                                                                                               aaGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTg                    <  2:647159/150‑1 (MQ=255)
                                                                                                                                 gCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCt                  <  2:1098380/150‑1 (MQ=255)
                                                                                                                                   atcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGCACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg                >  1:589384/4‑150 (MQ=255)
                                                                                                                                    tcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGc               >  2:636401/3‑150 (MQ=255)
                                                                                                                                     caCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCg              >  2:1586146/2‑150 (MQ=255)
                                                                                                                                       cGCCCCATCGCGGAGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTTGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt            <  1:689363/150‑1 (MQ=255)
                                                                                                                                       cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTCTCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt            >  1:1180647/1‑150 (MQ=255)
                                                                                                                                           ccAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCATCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGaa        >  1:571950/1‑150 (MQ=255)
                                                                                                                                                 cgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAGCAAAACCTTCGGCAGTGTGAAGGCGTTTGATAACGTCTGCTTGCGGTTGAATGCTgg  <  1:1498799/150‑1 (MQ=255)
                                                                                                                                                 cgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTgg  >  2:1600403/1‑150 (MQ=255)
                                                                                                                                                 cgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTgg  <  1:1246329/150‑1 (MQ=255)
                                                                                                                                                 cgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTgg  <  1:473230/150‑1 (MQ=255)
                                                                                                                                                 cgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGGCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTgg  <  2:1058659/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCAGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTGG  >  NC_000913/3731995‑3732289

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