Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,370,709 G→C 8.0% P110A (CCG→GCG) ‡ 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,7090GC 8.0% 368.5 / 17.9 112P110A (CCG→GCG) ‡nanEputative N‑acetylmannosamine‑6‑phosphate epimerase
Reads supporting (aligned to +/- strand):  ref base G (55/48);  new base C (3/6);  total (58/54)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 120 total aligned reads displayed.

GTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGG  >  NC_000913/3370561‑3370858
                                                                                                                                                    |                                                                                                                                                     
gTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACgcggg                                                                                                                                                      <  1:1219769/150‑1 (MQ=255)
 tGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACgcgggc                                                                                                                                                     <  2:1132640/150‑1 (MQ=255)
  ggTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACgcgggcg                                                                                                                                                    <  1:927862/150‑1 (MQ=255)
    tataGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggt                                                                                                                                                  >  1:127599/1‑150 (MQ=255)
     ataGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtc                                                                                                                                                 >  1:233848/1‑150 (MQ=255)
     ataGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtc                                                                                                                                                 >  2:775429/1‑150 (MQ=255)
      taGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcg                                                                                                                                                >  2:206035/1‑150 (MQ=255)
       aGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcgg                                                                                                                                               <  2:9641/150‑1 (MQ=255)
        gCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcggt                                                                                                                                              >  1:1158414/1‑150 (MQ=255)
        gCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcggt                                                                                                                                              <  1:1190278/150‑1 (MQ=255)
        gCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcggt                                                                                                                                              >  1:675486/1‑150 (MQ=255)
           aGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcc                                                                                                                                           <  1:1266581/150‑1 (MQ=255)
            gAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcct                                                                                                                                          <  1:676169/150‑2 (MQ=255)
             aaaGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcctg                                                                                                                                         <  2:1375324/150‑3 (MQ=255)
             aaaGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcctg                                                                                                                                         <  2:498850/150‑3 (MQ=255)
               aGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcctgtt                                                                                                                                       <  2:665320/150‑5 (MQ=255)
                 cGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAAt                                                                                                                                     <  1:68814/150‑1 (MQ=255)
                    aGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGc                                                                                                                                  <  2:177118/150‑1 (MQ=255)
                    aGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGc                                                                                                                                  <  2:930291/150‑1 (MQ=255)
                    aGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGc                                                                                                                                  <  1:1569710/150‑1 (MQ=255)
                      tGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATCGCGa                                                                                                                                <  2:46435/150‑1 (MQ=255)
                        ccAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATa                                                                                                                              >  2:85368/1‑150 (MQ=255)
                         cAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATaa                                                                                                                             <  1:770932/150‑1 (MQ=255)
                                  cGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCgcgc                                                                                                                    >  1:213587/1‑150 (MQ=255)
                                       cccAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGcgcc                                                                                                                   >  1:1336511/1‑146 (MQ=255)
                                       cccAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGcgcc                                                                                                                   <  2:1336511/146‑1 (MQ=255)
                                        ccAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCt                                                                                                              >  2:811143/1‑150 (MQ=255)
                                        ccAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCt                                                                                                              >  2:405569/1‑150 (MQ=255)
                                          aGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGCGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTgc                                                                                                            >  1:487609/1‑150 (MQ=255)
                                           gCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTgcg                                                                                                           <  2:613089/150‑1 (MQ=255)
                                             ttttGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGcc                                                                                                         <  1:904484/150‑1 (MQ=255)
                                                  gCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgcgc                                                                                                    >  2:386558/1‑150 (MQ=255)
                                                   cATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCa                                                                                                   >  1:1404441/1‑150 (MQ=255)
                                                     tGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATc                                                                                                 >  1:1295091/1‑150 (MQ=255)
                                                      gCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCa                                                                                                <  2:1006929/150‑1 (MQ=255)
                                                      gCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCa                                                                                                <  2:948518/150‑1 (MQ=255)
                                                        cAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTc                                                                                                                                        >  1:392043/1‑108 (MQ=255)
                                                        cAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTc                                                                                                                                        <  2:392043/108‑1 (MQ=255)
                                                         aGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACa                                                                                             <  1:902057/150‑1 (MQ=255)
                                                          ggCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTACCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACAt                                                                                            >  1:1016805/1‑150 (MQ=255)
                                                           gCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATc                                                                                           >  1:1271369/1‑150 (MQ=255)
                                                            ccGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCt                                                                                          <  2:127599/150‑1 (MQ=255)
                                                                cTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCaa                                                                                      <  1:1482605/150‑1 (MQ=255)
                                                                 ttCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAat                                                                                     <  2:1190070/150‑1 (MQ=255)
                                                                     ggCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGGCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAg                                                                                 <  2:192379/150‑1 (MQ=255)
                                                                       cGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGc                                                                               >  2:386745/1‑150 (MQ=255)
                                                                        gTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGcc                                                                              >  1:345344/1‑150 (MQ=255)
                                                                          tGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGt                                                                            <  2:1295091/150‑1 (MQ=255)
                                                                           gAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTg                                                                           <  1:1092059/150‑1 (MQ=255)
                                                                            aGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGa                                                                          >  1:1083449/1‑150 (MQ=255)
                                                                                gTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCg                                                                      >  2:1120323/1‑150 (MQ=255)
                                                                                gTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCg                                                                      <  2:1271369/150‑1 (MQ=255)
                                                                                   ggTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTAc                                                                   >  2:1503919/1‑150 (MQ=255)
                                                                                     tCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCg                                                                 <  2:767127/150‑1 (MQ=255)
                                                                                        tCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGgag                                                              <  2:163316/150‑1 (MQ=255)
                                                                                         cGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGgaga                                                             >  2:973908/1‑150 (MQ=255)
                                                                                         cGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGgag                                                              <  1:32029/149‑1 (MQ=255)
                                                                                         cGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGgag                                                              >  2:32029/1‑149 (MQ=255)
                                                                                          gCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGaa                                                            >  2:1623437/1‑150 (MQ=255)
                                                                                             aGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAAtcc                                                         >  1:692160/1‑150 (MQ=255)
                                                                                              gTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAg                                                                                                       >  1:1494595/1‑103 (MQ=255)
                                                                                              gTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAg                                                                                                       <  2:1494595/103‑1 (MQ=255)
                                                                                                aaaCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCCCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAAtcctcc                                                      >  2:1099430/1‑150 (MQ=255)
                                                                                                  aCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCTGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAg                                                    >  1:1298597/1‑150 (MQ=255)
                                                                                                  aCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAg                                                    <  1:145641/150‑1 (MQ=255)
                                                                                                        aTGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATcgcg                                              <  1:55986/150‑1 (MQ=255)
                                                                                                          ggTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGtt                                            <  1:812313/150‑1 (MQ=255)
                                                                                                            tGAATACGTTCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCTTCAATGGCTATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGATAATCCTCCAGATCGCGTTTc                                          <  1:419018/150‑1 (MQ=255)
                                                                                                             gAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTca                                         <  2:1383603/150‑1 (MQ=255)
                                                                                                                  cGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAAtt                                    <  1:1120323/150‑1 (MQ=255)
                                                                                                                  cGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAAtt                                    >  2:883533/1‑150 (MQ=255)
                                                                                                                    tGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAAtgca                                  <  1:386558/150‑4 (MQ=255)
                                                                                                                    tGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTcc                                  >  1:694389/1‑150 (MQ=255)
                                                                                                                    tGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTcc                                  <  1:26885/150‑1 (MQ=255)
                                                                                                                      ccAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCaa                                <  2:1158414/150‑1 (MQ=255)
                                                                                                                       cagcagCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAAGTCCaat                               >  2:84289/1‑150 (MQ=255)
                                                                                                                       cagcagCGTTTCAACAGGCACCGGACGCGGGCGGGCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCaat                               <  2:146486/150‑1 (MQ=255)
                                                                                                                        agcagcGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCaata                              >  1:1473495/1‑150 (MQ=255)
                                                                                                                           agcGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATc                           >  2:1576533/1‑150 (MQ=255)
                                                                                                                              gTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGc                        >  2:965707/1‑150 (MQ=255)
                                                                                                                              gTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGc                        <  2:692160/150‑1 (MQ=255)
                                                                                                                               tttCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCa                       >  2:114271/1‑150 (MQ=255)
                                                                                                                                 tCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACg                     >  1:138421/1‑150 (MQ=255)
                                                                                                                                  cAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGAGCAAGCGCATCAGCATGTTCAAAATAGGCCGTGATGGGGAACGGAGGATCCTGGATAGGGGGTTTGACATGTCCAATAATCGGCAAGc                    >  2:1199389/1‑150 (MQ=255)
                                                                                                                                   aaCAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCt                   >  1:805968/1‑150 (MQ=255)
                                                                                                                                     cAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTca                 >  2:1515511/1‑150 (MQ=255)
                                                                                                                                     acggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTca                 >  1:999935/3‑150 (MQ=255)
                                                                                                                                      aGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcac                <  1:1503919/150‑1 (MQ=255)
                                                                                                                                       ggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcacc               >  1:505428/1‑150 (MQ=255)
                                                                                                                                        gCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGGTTCACAATTCCAATAATCGGCCCGCTcacca              >  1:188058/1‑150 (MQ=255)
                                                                                                                                            cGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCa                                                     >  1:82278/1‑107 (MQ=255)
                                                                                                                                            cGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCa                                                     <  2:82278/107‑1 (MQ=255)
                                                                                                                                             ggACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGCATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCCCCGCa         >  2:661409/1‑150 (MQ=255)
                                                                                                                                             ggACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCCCCACCGCa         >  2:1361552/1‑150 (MQ=255)
                                                                                                                                               aCGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCAcg       <  2:406774/150‑1 (MQ=255)
                                                                                                                                                 gcgggcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATcgcg                                              <  1:470590/109‑1 (MQ=255)
                                                                                                                                                 gcgggcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCAcgcg     <  2:58218/150‑1 (MQ=255)
                                                                                                                                                 gcgggcggTCGCTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATcgcg                                              >  2:470590/1‑109 (MQ=255)
                                                                                                                                                  cgggtggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCGTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGt    >  1:820632/1‑150 (MQ=255)
                                                                                                                                                    ggcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTgg  >  1:1589020/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGG  >  NC_000913/3370561‑3370858

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: