Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,7750CT23.5% 92.0 / 11.9 51K199K (AAG→AAApitBmetal phosphate:H(+) symporter PitB
Reads supporting (aligned to +/- strand):  ref base C (11/28);  new base T (12/0);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 8.52e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGT  >  NC_000913/3135530‑3136024
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
ggCAGTTCAGGATGCTGGTGCAGGTAGTGTTCGAAGTTGGTAAAGGAATCGCGGTTACGGATAAATTCATAGCCGGACGCATTCTGATTGACGACGAATCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCTCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCttttt                                                                                                                                                                                                                                                          <  2:239440/251‑1 (MQ=255)
   aGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCACCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGTCGACGAAGCCAGCAGGGGCAACCCCCACCAGTACCAGCATTACCAGTCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCTCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCg                                                                                                                                                                                                                                                       <  1:390523/251‑1 (MQ=255)
        aGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCtt                                                                                                                                                                                                                                                  <  1:404118/251‑1 (MQ=255)
         ggATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGGTTATGTTTGCCCTttttttttttttgggcttt                                                                                                                                                                                                                                                 >  2:557681/1‑242 (MQ=255)
             gCTGTTGCAGGTAGTGTGCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGGCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCATCCgg                                                                                                                                                                                                                                              <  1:152346/250‑1 (MQ=255)
             gCTGTTGCAGGTAGAGTTCAAAGTTGGTAACGGCATCGCGGGTGCGGGTACTTTCATAGCCGGACGGATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGCCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAACCAATCAGCGCAATACGCGTCCAGAATGGCGGATTTTGTATTGCTCTTTTTTTTTTTTGcttcttctgg                                                                                                                                                                                                                                               >  2:103599/1‑242 (MQ=255)
             gCTGTGGCATGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGa                                                                                                                                                                                                                                             <  2:36162/251‑1 (MQ=255)
                                      ggTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTAGTTTGCCCTTTTTCTTTTTTCGGTCTTCCGGAATGCGGGGGATAAGGTCCCGCtttt                                                                                                                                                                                                                  >  1:110341/1‑251 (MQ=255)
                                         aaCGGCATCCCGGGGACGGGTAATTTCATAGCCGACCGAATTCATATTGACGACGAAGCCAGCAGGGTCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTc                                                                                                                                                                                                                <  2:523792/250‑1 (MQ=255)
                                            ggCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGATTTCGTTTTGCTTTTTTTTTTTTTTGCTCTTTCTGAGTGCGGTGAATTCGGGTACGCtttttttctc                                                                                                                                                                                                            >  2:247637/1‑246 (MQ=255)
                                                 cgcgGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCt                                                                                                                                                                                                        <  1:588527/250‑1 (MQ=255)
                                                        acgggtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCATCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGc                                                                                                                                                                                                <  2:306965/251‑1 (MQ=255)
                                                            gtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCg                                                                                                                                                                                            <  1:131997/251‑1 (MQ=255)
                                                                  tCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTATGTTTTCCTTTTTTTTTTTTTCGAGCTTTCGGAATGCGGTGAATAAGGTCACGCTTTTTCGTCCCGCTCCCGAAGCGTCgccgc                                                                                                                                                                                       >  1:191188/1‑246 (MQ=255)
                                                                          ggACGCATTCATATTGACGACGAAGCCAGCAGGGGCAAGCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCgg                                                                                                                                                                                                                                              <  2:636392/189‑1 (MQ=255)
                                                                                aTTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTATGTTTTCCTTTTTTTTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTGACGCTTTTTTGTCCCCCTCCCGTAGCGTCGCAGCCGGAATATCAGGCCTc                                                                                                                                                                       >  1:382111/1‑251 (MQ=255)
                                                                                aTTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAGGTTGCCTTTTTTTTTTTTTCGATCTTCCGGAATGCGGGGGATACGGTCACGCTTTTTTGTCCCGCTCCAGTAGCGGCGCAGCAGGAATATCAGGCCTc                                                                                                                                                                       >  1:329592/1‑251 (MQ=255)
                                                                                   aatatTGACGACGAAGCCAGACGGGGCAATCCCTCCCAGTCCTAGCATTACCAGGCCGATCCCGTTTTGTCCGTCGTTCGCGCCGGGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTGTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTGCGTCCCGCTCCAGTAGCGTCGCAGCAGGTATATCAGGCCTCCCg                                                                                                                                                                    <  2:416305/250‑1 (MQ=255)
                                                                                      atTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACAGTTGCCTTTTTTTTTTTTTTGTTCTTCCGGAATACGGGGGATACGGGCCACCTTTTTCGTCCCCCTCCAGGAGCGGCGCAGCAGGGATATCAgacatcccccc                                                                                                                                                                  >  1:233949/1‑240 (MQ=255)
                                                                                                         ggggCAACCCCCCCCAGTACCAGCATTACCAGGCCGAGCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAgg                                                                                                                                              <  2:329042/251‑1 (MQ=255)
                                                                                                             cAATCCCCACCAGTACCAGCATTACCAGGCCGACCCCTGTTTGGCCGTCGTTCGCGCCGTGCGAAAACGCCACGGCCGCAGCGGAAACAATCAGCGCAATTCGCGGCCATAATGGCGGTTTCCGTTTGCTGTTTTTCTTTTTGCGAGCTTCCGGAATGCGGTGAATACGGTCCCGCTTTTTCGTCCCGCTCCAGTAGCGGCGCAGCAGGAATATCAGGCCTCTCGCAATGACCAGGCCGACGATAGGGGaa                                                                                                                                          <  1:357036/251‑1 (MQ=255)
                                                                                                                         gTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTAGTTTGCCTTTTTTCTTTTTGCGCTCTTTCGGGATGCGGTGGATACGGTTACGCTTTTTTGTCCCCCCACAGTCGCGTCGCAGCAGGAATATCTGGGCCCCCCCAATTACCCAGCCGACGATAGGGGGAACaaaaaaag                                                                                                                               >  2:452531/1‑243 (MQ=255)
                                                                                                                          tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGGTTTGCTTTTTTTTTTTTGCGATCTTCCTGAATGAGGGGAATACGGTCACGCTTTTTTGTCCCCCTCCCCGAGCGTCGCAGCAGGGAAATCAAGGCTCCCCCAAAGACCCGGCCGACGATAGGGGAAAACAACAGCgagg                                                                                                                            >  1:428004/1‑251 (MQ=255)
                                                                                                                          tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTTGCTTTTTTTTTTTTTTGGTCTTTCGGGATGCGGGGAATACGGTCAATCTTTTTCGTCCCGCCCCCGTAGCGTCGCCGCAGGGATAATAGGCCGCCCGCCATGACCAGGCCGACCGCAGGGGAAAAAAACAGCgagg                                                                                                                            >  1:279067/1‑251 (MQ=255)
                                                                                                                                      cAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTCGTTTGCCTTTTTTTTTTTTTCGATCTTCCGGGATGCGGTGGATACGGTCACGCTTTTTCGTCCCGCCCCCGTAGCGTCGCAGCAGGGATATCAGGCCCCCCGCCATGACCCGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTgg                                                                                                                >  1:118958/1‑251 (MQ=255)
                                                                                                                                       aGGCCTATCCCTTTTTGGCCGCCTTCCGCGCCTTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGt                                                                                                               <  2:118958/251‑1 (MQ=255)
                                                                                                                                                 cTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACATTTTCCCTTTTTCTTTTTTAGATCTTCCGGAAAGCGGTGAAAAACGGCACGCTTTTTTGTCCCCCTCCAGTAGCGTCGGCGGCAGAATATCAGGACCCCCCCAATGACCACGCCGACGCCACGGGGAACAAAAAACGAGGAGAAAATTTTTGTccactcacgc                                                                                                     >  1:316804/1‑242 (MQ=255)
                                                                                                                                                   ttttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCa                                                                                                    <  1:233536/250‑1 (MQ=255)
                                                                                                                                                    tttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGGTTTCCTTTTTTTTTTTTGCGGTCTTCCGGAATGCGGTGGATACGGTCACGCTTTTTCGTCCCCCCCCAGGAGCGTCGCAGCAGGAATATCAGGCCTCCCGCCAAGACCCGGCCGGCGATTGGGGGAACAAACAGCGAGGGGAAAAATTTTGTCGCTTCACGCAgg                                                                                                  >  1:620486/1‑251 (MQ=255)
                                                                                                                                                     ttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTTGATCTTCCGGGATGCGGTGGATa                                                                                                                                                                                                                                 >  2:255574/1‑127 (MQ=255)
                                                                                                                                                     ttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATa                                                                                                                                                                                                                                 <  1:255574/127‑1 (MQ=255)
                                                                                                                                                      tGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGt                                                                                                 <  1:554203/250‑1 (MQ=255)
                                                                                                                                                                          cGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATg                                                                            <  2:329592/251‑1 (MQ=255)
                                                                                                                                                                            aaaaCGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGAt                                                                             <  1:165479/248‑1 (MQ=255)
                                                                                                                                                                               acgccacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAGTTTGCCTTTTTTTTTTTTTCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCCGCAGGAAAATTAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGATTAAAGCATCCATCACCggtggaccg                                                                       >  1:44031/1‑243 (MQ=255)
                                                                                                                                                                                   cacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTg                                                                                                                 <  2:596083/205‑1 (MQ=255)
                                                                                                                                                                                           cAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGAATGCGGGGGATAAGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGCTAGGGGGAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTTACGCATCCATCACCGAAGAGCCCGGTAACACcgcg                                                           >  1:408853/1‑251 (MQ=255)
                                                                                                                                                                                               ggAAACAATCAGCGCAATACGCGTCCAGAAAGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGCCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTgg                                                       <  2:335513/251‑1 (MQ=255)
                                                                                                                                                                                                    cAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTaaa                                                  <  2:620486/251‑1 (MQ=255)
                                                                                                                                                                                                        cAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGa                                              <  1:247637/251‑1 (MQ=255)
                                                                                                                                                                                                           cgcAATACGCGTCCAGAATGGCGGTTTTCGTTTTCCTTTTTTCTTTTTGCGCTCTTCCGGAATGCGGGGGATACGGTCACGCTTTTTTGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGGCTCCCGCCATGAACAGGCCGACGATAGGGGGAACAATCAGCGGGGAGGAAAATTTGGGCACTTCACGCAGGTTTAACGCATCCCTCTACGATGAGCCGGTTAACAGCGCGTTGGGTAAACCCATGc                                           >  1:462639/1‑251 (MQ=255)
                                                                                                                                                                                                                   gcgTCCAGAAGGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATc                                   <  1:103599/251‑1 (MQ=255)
                                                                                                                                                                                                                         aGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTGCGTCCCGCGCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCa                             <  2:408853/251‑1 (MQ=255)
                                                                                                                                                                                                                             tGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCg                                                                               >  1:121069/1‑197 (MQ=255)
                                                                                                                                                                                                                             tGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGGTTAACGCATCCATCACCg                                                                               <  2:121069/197‑1 (MQ=255)
                                                                                                                                                                                                                              ggCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCa                         >  1:142868/1‑250 (MQ=255)
                                                                                                                                                                                                                               gCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAg                       <  1:452531/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGgt                     >  1:515063/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cGTTTGCCTTTTTTCTTTTGGCGATCTACCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCg                <  1:176938/250‑1 (MQ=255)
                                                                                                                                                                                                                                          tttGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAAc             >  1:351703/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ccTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTgg          >  1:67330/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   tttttCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAACATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATTAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACCGGAGGc      >  1:555938/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        cTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGAAGAGCCGGTTAACAGCGCGTTGGTTAAACCGCTGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGt  >  2:589771/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
GGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGT  >  NC_000913/3135530‑3136024

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