Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,370,8750AC31.9% 138.1 / 10.6 70G54G (GGT→GGGnanEputative N‑acetylmannosamine‑6‑phosphate 2‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (35/12);  new base C (0/22);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 9.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCA  >  NC_000913/3370629‑3371123
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         
gggatttgaacagtgtgtctgagCCAGTAACCCGTGATGGTGAATACGGGCCAGCAGCGTTTCAACTGGCACCGGTCGCGGGCGGGCGGTGCCGTCAAGGGGGATAATGCCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCTCGCTCACCAACGCACGCGTGGCTTGCAGATTTGCCACACCtt                                                                                                                                                                                                                                                       <  2:69397/229‑1 (MQ=255)
   cGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCa                                                                                                                                                                                                                                                     >  1:177088/1‑250 (MQ=255)
       gAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGc                                                                                                                                                                                                                                                <  1:349573/250‑1 (MQ=255)
                 ccATCGCCAGTAACCCGTGAGGGTGAATACGTGCCAGAAGCTTTTCAACAGGCACCGGACGCGGGCTGGCGGTGCCGTCAATGGCCATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAAct                                                                                                                                                                                                                                     <  2:5036/250‑2 (MQ=255)
                      gCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGcgcc                                                                                                                                                                                                                                >  2:16164/1‑251 (MQ=255)
                       ccAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCCCCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGcgccc                                                                                                                                                                                                                               <  2:67807/251‑1 (MQ=255)
                           ttaaCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCt                                                                                                                                                                                                                           >  1:203970/3‑251 (MQ=255)
                           tAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCt                                                                                                                                                                                                                           >  1:203995/1‑251 (MQ=255)
                            aaaCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTg                                                                                                                                                                                                                          <  1:601055/251‑1 (MQ=255)
                                   gatggtaattccgtcccAGCAGCTTTTCAACAGGCCCCGGAGGGGGGGGGTCGGTGCCTTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCCCCGCACGCGTGGCTTGCAGATTTGCCACCCCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGcc                                                                                                                                                                                                                   <  1:120972/236‑1 (MQ=255)
                                      ggTGAATGCGGGCCAGCAGCGTTTCAACAGGCCCGGGGCGGGGGCGGTCGGTGCCGTAAATGGCGATAATGCCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGCCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACCCCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCt                                                                                                                                                                                                                <  1:466842/251‑1 (MQ=255)
                                             acggcccggcagcGTTCCAACGCGCCCCGGCCGCGGGCGTGCGGTGCGGTTAAGGGCATTAATGTCCGCCCCCGCCCGCGCCAGCGCATCAACATCTTCAAAAAAGGCCGGGATGCGTCCCGGAGAATCCTCCAGATCGCGTTCCCCAATTCCAATAAGCGTCGCGCCCCCCACCGCCCGCGTGGCTTGCAGATTTCCCACCCCTTCAATGCGAATGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGcc                                                                                                                                                                                                          <  1:363511/242‑1 (MQ=255)
                                                                 cagccccgggccgggggcggTCGGGGCCGTAAAGGGCGATAATGCCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGt                                                                                                                                                                                     <  1:16164/237‑1 (MQ=255)
                                                                   ggCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGttt                                                                                                                                                                                   >  2:313971/1‑251 (MQ=255)
                                                                     cACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTAt                                                                                                                                                                                 >  1:390221/1‑251 (MQ=255)
                                                                                   gatcggtGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGAGTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTACGGGTTTATCGAGCGGGCTGTcc                                                                                                                                                                   >  2:88494/3‑251 (MQ=255)
                                                                                    gtcggtGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTccc                                                                                                                                                                  >  2:170116/1‑250 (MQ=255)
                                                                                       ggtGCCTTCAAGGGCGATATTGCCCGCGCCCGCCGGCGCCAGCGCATCAACATCTTCAATATAGGCCTTGATGCGTCCCGGAGACCCCTCCAGATCGCGTTTCCCAATTCCAATAACCGGCCCGCTCCCCACCGCACGCGTGGCTTGCAGATTTGCCACCCCTTCAATGCGAATGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGaa                                                                                                                                                               <  2:29842/251‑1 (MQ=255)
                                                                                         tGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAAcc                                                                                                                                                             >  2:106666/1‑251 (MQ=255)
                                                                                                   gggattatttcccggccccgccTGCGCCAGCACACAAACATCTACAATAGGGGCCGTGAGGGGTCCGGAAAAACCCTCCAGATCGGGTTTATCAATTCCAAGAATCGCCCCGCTCACCACCGCAGGCGTGGCTTGCAGATTTGCCACCCCTACAATGCGAATGGCAGCCGCTCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAgg                                                                                                                                                   <  1:522373/235‑1 (MQ=255)
                                                                                                             cccgatccccccTGCGCCAGCGCATCAACATCTCCAATATAGCCCGTGATGCGTCCCGGAGAACCCTCCAGATCGCGTTTCACAATCCCAATAACCGGCCCGCTCACCACCGCACGCGTGGCTTGCAGATTTCCCCCCCCTTCAATGCGAAGGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAg                                                                                                                                         <  2:36105/241‑1 (MQ=255)
                                                                                                              ccGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTTGCCGGAACCGGCTGGCAGGAGACAATCAgg                                                                                                                                        >  2:446887/1‑251 (MQ=255)
                                                                                                                  ccccgccTGCGCCAGCGCACCAACATCTTCATTATAGCCCGTGGTGCGACCCGGAGAATCCACCAGATCGTTTTTCACAATTCAAATAATCGGCCCGCTCACCACCGCACGCTTGGCTTGCAGATTTGCCACCCCTTCAAGGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCAc                                                                                                                                    <  2:66707/250‑1 (MQ=255)
                                                                                                                     cgccTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAg                                                                                                                                         >  1:36105/1‑243 (MQ=255)
                                                                                                                     cgccTGCGCCAGCGCATAAACATCTTCAAGATAGGCCTTGATGCTGACGGAAAACCCCTCCAGACCGGTTTCCACAATCCCAAAAACCGGCACCCTCCCCCCCGCCCGCGTGTGTTGCAGATTTCCCACCCCTACAAGGCGAATGTCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGt                                                                                                                                 <  2:111375/251‑1 (MQ=255)
                                                                                                                      accTGCGCCAGCGCACAACAACCTTAAAAATAGCCCGTGATCCGTAGCGGATAACCCCCCAGCTCTGGTTTCGCAATTCCAATACCCGGCACGCCAACCACCGCACGCGTGGCTTGAAGTTTTCCCCCACCTTCAAGGCGAAGGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGtt                                                                                                                                <  2:86600/250‑1 (MQ=255)
                                                                                                                            gcCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTagcagc                                                                                                                          >  1:342909/1‑251 (MQ=255)
                                                                                                                               ngcgcATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGAt                                                                                                                       >  2:270693/2‑251 (MQ=255)
                                                                                                                               agcccataaaCATCTTAAATATAGGCCGTGATGCGTACCGGAGAACCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCCCCCCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACTGTTAGCAGCGAt                                                                                                                       <  1:446887/243‑1 (MQ=255)
                                                                                                                                  gcATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCCGGCAGGAGACAATCAGGCCACCGTTCGCAGCGAgttt                                                                                                                    >  1:351781/1‑251 (MQ=255)
                                                                                                                                         cATCTTCAATATAGGCCTTGATGCGTACCGGAAAACCCTCCAGATCGCTTTTCAAAATTCAAATAGTCGGCACCCCCACCACCGCACGCGTGGCTTGCAGATTTCCCACCCCTTAAAGGCGAATGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATcc                                                                                                             <  1:558518/251‑1 (MQ=255)
                                                                                                                                         aactctaaaatatagccctggATGCGTCCCGAGGAATCCTCCAGACCGCGTTTCACAATTCCAAAAACCGGCACGCCCCCCACCGCCCGCGTGGCTTGCGGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATc                                                                                                              <  2:252346/231‑1 (MQ=255)
                                                                                                                                                    aggggCCGGGATGCGTACCGGAGAATCCCCCAGAGCGCGTTTCACATTCCCAAAAACCGGCACCCCCACCACCGCGCGCGTGGTTTGCAGATTTCCCCCCCCTTAAATGCGAAGGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCaa                                                                                                   <  2:177088/248‑1 (MQ=255)
                                                                                                                                                      aGGCCGTGATGCGTACCGGAGAATCCTACAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGCGACAATCAGGCCACCGGTAGCAGCGATTTTTTGATTCAGTTGTGCCAGt                                                                                                 >  2:115840/1‑251 (MQ=255)
                                                                                                                                                       ggCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGt                                                                                                 >  1:44722/1‑250 (MQ=255)
                                                                                                                                                       ggCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCACCGATTTTTTGATCCAGTTGTGCAAGt                                                                                                 >  2:301995/1‑250 (MQ=255)
                                                                                                                                                           gTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACATTTCCAATAATCGGCACGCTCACCCCCGCACGCGTGGCTTGCAGATTTGCCACCCCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGa                                                                                            <  1:406328/251‑1 (MQ=255)
                                                                                                                                                                gCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAACCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGacatac                                                                                       <  1:106666/251‑1 (MQ=255)
                                                                                                                                                                  gaaCCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGacatacat                                                                                     >  1:212800/3‑251 (MQ=255)
                                                                                                                                                                   tACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATc                                                                                    >  2:237742/1‑251 (MQ=255)
                                                                                                                                                                      cGGAAAATCCTCCAGATCGCGTTTCACAATTCCAAAAATCGGCCCGCCCACCACCGCACGCGTGGCTTGCAGATTGGCCACCCCTTAAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTc                                                                                 <  1:420384/251‑1 (MQ=255)
                                                                                                                                                                        gggaaaaccctccggcgcgcgTTTCCCAATTCCAATAATCGGCCCGCTCACCCCCGCCCGCGTGGCTTGCAGATTTGCCACCCCTTCAATGCGAATGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTc                                                                                 <  2:244654/234‑1 (MQ=255)
                                                                                                                                                                            gagcataatgtatcgcatggcacaATTCCAATAATCGGCACGCTCACCCCCGCACGCGTGGCTTGCAGATTTGCCACCCCTTCAATGCGAATGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCt                                                                              <  1:534367/230‑1 (MQ=255)
                                                                                                                                                                              tcctccAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGc                                                                          >  1:165189/1‑251 (MQ=255)
                                                                                                                                                                              tcctccAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGc                                                                          >  1:402882/1‑251 (MQ=255)
                                                                                                                                                                                ctccAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGCATGACAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAAACAATCAGGCCACCGGTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTACCCTTAg                                                                           >  2:20252/1‑248 (MQ=255)
                                                                                                                                                                                ctccAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCg                                                                         >  1:594927/1‑250 (MQ=255)
                                                                                                                                                                                  ccaaactcgtgTTTCCCAATCCCAATAATCGCCCCCCTCCCCCCCCCACGCGGTGCTTGCAGTTTTGCCCCCCCTTCAAGGCGAATGGCACCCGCCCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGan                                                                       <  2:276492/241‑2 (MQ=255)
                                                                                                                                                                                            gTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGGGGCTTACAGATTTGCAACACCTTCAATGCGAATGGCAACCGCGCCCGCGTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGATGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGCTAGCAGCGATTTTTTTATCCCGATGTGCACGTAAAGACATACATCATCCCTTAGCGaatggtacgggg                                                             >  2:400182/1‑241 (MQ=255)
                                                                                                                                                                                                   caATTCCCAAAAATCGGCACGCCCCCCCCCGCCCGCGGGGCTGGCAGTTTTGCACCCCCTCCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAg                                                        <  2:44722/251‑1 (MQ=255)
                                                                                                                                                                                                     aTTCCAATAACCGGCCCCCTCACCACCGCCCGCGTGGCTTGCAGATTTGCCCCCCCTTCAATGCGAATGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGAcc                                                     <  1:237742/251‑1 (MQ=255)
                                                                                                                                                                                                        ccAATAATCGGCACGCTCCCCACCGCACGCGGGGCTTGCAGATTTGCCCCCCCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTTCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGCACATAGACCggg                                                  <  1:431690/251‑1 (MQ=255)
                                                                                                                                                                                                         cAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGc                                                 >  1:537110/1‑251 (MQ=255)
                                                                                                                                                                                                                cgggCACGCTCCCCCCCGCACGCGTGGTTTGCAGTTTTGCCCCCCCTGAAATGCAAATGGCAACCGCTCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAgg                                           <  1:525377/249‑1 (MQ=255)
                                                                                                                                                                                                                ccGGCACGCCCCCCACCGCCCGCGTGGTTTGCAGATTTGCCTCCCCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAgg                                           <  1:313971/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGAt                                         >  1:130102/1‑251 (MQ=255)
                                                                                                                                                                                                                   gcccGCTCCCCCCCGCACGCGTGGCTTGCAGATTGGCCACACCTTAAATGCGAAGGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGAtt                                        <  1:212596/248‑1 (MQ=255)
                                                                                                                                                                                                                   gCACGCTCACCACCGCACGCGTGGCTTGCATATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGAtt                                        >  2:106046/1‑251 (MQ=255)
                                                                                                                                                                                                                    cACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGAtt                                        >  2:120488/1‑250 (MQ=255)
                                                                                                                                                                                                                          cccccccGCGCGCGTGGCTTGAAGATTTCCCCCACCTTCAAGGAGAAGGCCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACtttt                                 <  2:332114/246‑1 (MQ=255)
                                                                                                                                                                                                                          caccccCGCACGCGTGGCTTGCAGTTTTGCCCCACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACtt                                   <  1:134516/249‑1 (MQ=255)
                                                                                                                                                                                                                            ccaccGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTgg                               >  1:53107/1‑251 (MQ=255)
                                                                                                                                                                                                                             caccGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGt                              >  2:446807/1‑251 (MQ=255)
                                                                                                                                                                                                                              accGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGtt                             >  1:497410/1‑251 (MQ=255)
                                                                                                                                                                                                                               ccGCACGCGGGGCTTGCAGATTTGCCACCCCTTCAATGCGAATGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGttt                            <  1:279195/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cccgcgGGGCTTGAAGATTCGCCACCCCTCAAGTGCGAGGGGCACCCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGa                         <  2:130102/249‑1 (MQ=255)
                                                                                                                                                                                                                                       gTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGGTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTaa                     >  2:303972/1‑250 (MQ=255)
                                                                                                                                                                                                                                        tGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAAt                   >  2:298096/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTCCCGCTAATGACATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGttttt            >  1:343904/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTg           >  1:245234/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         acacCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGCAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGcc   >  2:528439/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          cacCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCa  >  2:539505/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          cacCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGGTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATCCATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTACCTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCa  >  2:85255/1‑250 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         
CGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCA  >  NC_000913/3370629‑3371123

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