Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,370,714 T→C 5.7% D108G (GAC→GGC) ‡ nanE ← putative N‑acetylmannosamine‑6‑phosphate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,370,7140TC 5.7% 484.6 / 13.1 162D108G (GAC→GGC) ‡nanEputative N‑acetylmannosamine‑6‑phosphate epimerase
Reads supporting (aligned to +/- strand):  ref base T (81/69);  new base C (6/3);  total (87/75)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

Only 100 of 178 total aligned reads displayed.

TATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGC  >  NC_000913/3370565‑3370863
                                                                                                                                                     |                                                                                                                                                     
tataGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggt                                                                                                                                                       >  1:2480872/1‑150 (MQ=255)
  taGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcg                                                                                                                                                     >  1:1279207/1‑150 (MQ=255)
    gCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAACCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcggt                                                                                                                                                   >  1:2186232/1‑150 (MQ=255)
     ccAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTg                                                                                                                                                  <  1:2156931/150‑1 (MQ=255)
     ccAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTg                                                                                                                                                  <  1:1496240/150‑1 (MQ=255)
      cAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGc                                                                                                                                                 <  2:1799677/150‑1 (MQ=255)
       aGAAAGCGCAGTGCCAATACTTTCGGCGCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcc                                                                                                                                                <  2:1405887/150‑1 (MQ=255)
        gAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCg                                                                                                                                               <  1:1456512/150‑1 (MQ=255)
        gAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCg                                                                                                                                               <  2:11998/150‑1 (MQ=255)
          aaGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTc                                                                                                                                             <  2:389302/150‑1 (MQ=255)
          aaGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcctgt                                                                                                                                             <  1:2301566/150‑4 (MQ=255)
            gCGTAGTGCCAAAAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGAAGGTGAATACGTGCCAGCAGCGTTTCCACAGGCACCGGCCGCGGGCGGTCGGTCCCGTCaa                                                                                                                                           >  2:2355492/1‑150 (MQ=255)
              gTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATg                                                                                                                                         <  2:829082/150‑1 (MQ=255)
              gTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATg                                                                                                                                         <  2:20970/150‑1 (MQ=255)
               tAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATgg                                                                                                                                        >  1:406496/1‑150 (MQ=255)
                aGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGc                                                                                                                                       <  1:1194785/150‑1 (MQ=255)
                aGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGc                                                                                                                                       <  2:1432161/150‑1 (MQ=255)
                 gTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCg                                                                                                                                      >  1:1266627/1‑150 (MQ=255)
                     cAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGCCGCGGGCGGTCGGTGCCGTCAATGGCGATaa                                                                                                                                  >  2:1295243/1‑150 (MQ=255)
                     cAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATaa                                                                                                                                  >  1:1553959/1‑150 (MQ=255)
                     cAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATaa                                                                                                                                  >  2:1479413/1‑150 (MQ=255)
                          atTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTcc                                                                                                                             >  2:1438108/1‑150 (MQ=255)
                           tttCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCcg                                                                                                                            <  2:575972/150‑1 (MQ=255)
                           tttCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCcg                                                                                                                            <  2:384477/150‑1 (MQ=255)
                              cGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCgcgc                                                                                                                         >  2:1804247/1‑150 (MQ=255)
                                  tCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGcgcccgc                                                                                                                     >  1:125922/1‑150 (MQ=255)
                                  tCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGcgcccgc                                                                                                                     >  2:2406379/1‑150 (MQ=255)
                                       gCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCTTGATGGTGAATACGTGCCAGCAGCGTTTCAAGAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTgcg                                                                                                                <  2:877556/150‑1 (MQ=255)
                                       gCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTgcg                                                                                                                >  1:932941/1‑150 (MQ=255)
                                           ttGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAg                                                                                                            <  1:2170429/150‑1 (MQ=255)
                                            tGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgc                                                                                                           <  2:1450826/150‑1 (MQ=255)
                                             ggCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgcg                                                                                                          >  2:1137063/1‑150 (MQ=255)
                                             ggCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgcg                                                                                                          >  1:371466/1‑150 (MQ=255)
                                             ggCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgcg                                                                                                          >  1:219661/1‑150 (MQ=255)
                                                 tGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATc                                                                                                      <  2:162013/150‑1 (MQ=255)
                                                   ccAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCaa                                                                                                    >  2:327216/1‑150 (MQ=255)
                                                    cAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGACAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAAc                                                                                                   <  1:100967/150‑1 (MQ=255)
                                                     aGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACa                                                                                                  <  2:1906713/150‑1 (MQ=255)
                                                     aGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACa                                                                                                  >  1:748978/1‑150 (MQ=255)
                                                         cGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGGCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCtt                                                                                              <  2:1719673/150‑1 (MQ=255)
                                                          gTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTc                                                                                             <  1:1227399/150‑1 (MQ=255)
                                                            cTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCaa                                                                                           <  2:117572/150‑1 (MQ=255)
                                                            cTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCaa                                                                                           >  2:428916/1‑150 (MQ=255)
                                                            cTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCaa                                                                                           <  1:1780446/150‑1 (MQ=255)
                                                              tCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAata                                                                                         >  1:2034914/1‑150 (MQ=255)
                                                               ccGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCACATCAACATCTTCAatat                                                                                        >  2:2228899/1‑150 (MQ=255)
                                                                 ggCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAg                                                                                      >  1:2386881/1‑150 (MQ=255)
                                                                  gCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAgg                                                                                     <  1:885623/150‑1 (MQ=255)
                                                                   cGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGc                                                                                    <  1:1969464/150‑1 (MQ=255)
                                                                    gTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGcc                                                                                   <  2:1235770/150‑1 (MQ=255)
                                                                         gCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGAt                                                                              <  1:2032453/150‑1 (MQ=255)
                                                                         gCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGAt                                                                              <  2:2480872/150‑1 (MQ=255)
                                                                          cAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATg                                                                             <  1:984491/150‑1 (MQ=255)
                                                                           aGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGc                                                                            >  1:1541905/1‑150 (MQ=255)
                                                                               ggTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTAc                                                                        >  1:1116938/1‑150 (MQ=255)
                                                                                  cATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCgg                                                                     <  1:1140124/150‑1 (MQ=255)
                                                                                  cATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCgg                                                                     >  1:2375907/1‑150 (MQ=255)
                                                                                   aTCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGga                                                                    >  2:510031/1‑150 (MQ=255)
                                                                                    tCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGgag                                                                   <  2:1274630/150‑1 (MQ=255)
                                                                                     cGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGgaga                                                                  >  2:550609/1‑150 (MQ=255)
                                                                                      gCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGaa                                                                 >  2:1108155/1‑150 (MQ=255)
                                                                                      gCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCACCGCATAAACAACTTCAAAAAAGGCCGTCATGCGTACCGGAGaa                                                                 >  2:1997016/1‑150 (MQ=255)
                                                                                       ccAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGGCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAAt                                                                <  1:1842685/150‑1 (MQ=255)
                                                                                         aGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAAtcc                                                              >  2:2046589/1‑150 (MQ=255)
                                                                                          gTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAAtcct                                                             >  1:1578710/1‑150 (MQ=255)
                                                                                             aaCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCa                                                          >  1:1543150/1‑150 (MQ=255)
                                                                                              aCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAg                                                         <  2:125922/150‑1 (MQ=255)
                                                                                                     tGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGt                                                  >  2:2185969/1‑150 (MQ=255)
                                                                                                        tGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTc                                               >  1:455471/1‑150 (MQ=255)
                                                                                                         gAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTca                                              <  1:2227691/150‑1 (MQ=255)
                                                                                                             aCGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAAt                                          <  2:1543150/150‑1 (MQ=255)
                                                                                                              cGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAAtt                                         >  2:1038788/1‑150 (MQ=255)
                                                                                                              cGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAAtt                                         >  2:617154/1‑150 (MQ=255)
                                                                                                                 gCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCa                                      >  1:899647/1‑150 (MQ=255)
                                                                                                                  ccAGCAGCGTTTCAACAGGCACCGGACGCTGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCaa                                     <  1:569256/150‑1 (MQ=255)
                                                                                                                  ccAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCaa                                     <  1:1801353/150‑1 (MQ=255)
                                                                                                                          gTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACTGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGc                             >  2:959071/1‑150 (MQ=255)
                                                                                                                            cattgacggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCAc                           >  1:150594/8‑150 (MQ=255)
                                                                                                                             tCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACg                          <  2:2310110/150‑1 (MQ=255)
                                                                                                                             attgacggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCGTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGGTTCCCAATTCCAATAATCGGCACg                          >  2:759884/7‑150 (MQ=255)
                                                                                                                              cAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGc                         <  2:1916101/150‑1 (MQ=255)
                                                                                                                              cAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGc                         >  2:1512368/1‑150 (MQ=255)
                                                                                                                              cAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGc                         >  2:1338878/1‑150 (MQ=255)
                                                                                                                              cAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGc                         >  1:85803/1‑150 (MQ=255)
                                                                                                                                 cAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcc                      >  2:886885/1‑149 (MQ=255)
                                                                                                                                 cAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTca                      <  1:706463/150‑1 (MQ=255)
                                                                                                                                 acggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGTCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCCCAATTCCAATAATCGGCACGCTca                      >  1:2111227/3‑150 (MQ=255)
                                                                                                                                 acggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcc                      >  1:597195/3‑149 (MQ=255)
                                                                                                                                  cggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCGTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcac                     >  1:1015986/2‑150 (MQ=255)
                                                                                                                                  aGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcac                     >  1:2218170/1‑150 (MQ=255)
                                                                                                                                       ccGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCACCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCaccaccc                >  2:2168376/1‑149 (MQ=255)
                                                                                                                                          gACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGGATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCGCACCACCGCAc             >  2:1737406/1‑150 (MQ=255)
                                                                                                                                              cgggcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGt         >  2:7786/1‑150 (MQ=255)
                                                                                                                                                ggcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTgg       <  1:1479413/150‑1 (MQ=255)
                                                                                                                                                 gcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGc      <  1:703570/150‑1 (MQ=255)
                                                                                                                                                   ggtcggtGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCtt    >  2:371319/1‑150 (MQ=255)
                                                                                                                                                   ggtcggtGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCtt    >  1:1708949/1‑150 (MQ=255)
                                                                                                                                                    gtcggtGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTg   <  1:864119/150‑1 (MQ=255)
                                                                                                                                                    gtcggtGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGATTCCCAATTCCAATAATCGGCACGCTCCCCACCGCAAGCGTGGCtag   >  2:1058787/1‑148 (MQ=255)
                                                                                                                                                    gtcggtGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCAAGCGCATCAACATCTTCAAAATAGGAAGAGATGCGTACCGGAGAATCCTCCAGATCCCGGTTCAACGATCCAATAATCCGCACGATCCCCCCCGCAAGAGTGGCATg   >  2:832588/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGC  >  NC_000913/3370565‑3370863

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