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→AAA) pitBmetal 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.

GGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGT  >  NC_000913/3135530‑3136024
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
ggCAGTTCAGGATGCTGGTGCAGGTAGTGTTCGAAGTTGGTAAAGGAATCGCGGTTACGGATAAATTCATAGCCGGACGCATTCTGATTGACGACGAATCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCTCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCttttt                                                                                                                                                                                                                                                          <  2:239440/251‑1 (MQ=255)
   aGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCACCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGTCGACGAAGCCAGCAGGGGCAACCCCCACCAGTACCAGCATTACCAGTCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCTCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCg                                                                                                                                                                                                                                                       <  1:390523/251‑1 (MQ=255)
        aGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCtt                                                                                                                                                                                                                                                  <  1:404118/251‑1 (MQ=255)
         ggATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGGTTATG‑T‑TTGC‑CCTttttttttttttgggcttt                                                                                                                                                                                                                                                 >  2:557681/1‑242 (MQ=255)
             gCTGTTGCAGGTAGTGTGCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGGCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCATCCgg                                                                                                                                                                                                                                              <  1:152346/250‑1 (MQ=255)
             gCTGTTGCAGGTAGAGTTCAAAGTTGGTAACGGCATCGCGGGTGCGGGTACTTTCATAGCCGGACGGATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGCCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAACCAATCAGCGCAATACGCGTCCAGAATGGCGGATTTTG‑TATTGCTCTTTTTTTTTTTTGcttcttctgg                                                                                                                                                                                                                                               >  2:103599/1‑242 (MQ=255)
             gCTGTGGCATGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTTCGATCTTCCGGa                                                                                                                                                                                                                                             <  2:36162/251‑1 (MQ=255)
                                      ggTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTAG‑T‑TTGC‑CCTTTTTCTTTTTTCGGTCTTCCGGAATGCGGGGGATA‑AGG‑TCCCGCtttt                                                                                                                                                                                                                  >  1:110341/1‑251 (MQ=255)
                                         aaCGGCATCCCGGGGACGGGTAATTTCATAGCCGACCGAATTCATATTGACGACGAAGCCAGCAGGGTCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTc                                                                                                                                                                                                                <  2:523792/250‑1 (MQ=255)
                                            ggCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGATTTCG‑T‑TTTG‑CTTTTTTTTTTTTTTGCTCTTTCTGAGTGCGGTGAATT‑CGG‑GTACGCtttttttctc                                                                                                                                                                                                            >  2:247637/1‑246 (MQ=255)
                                                 cgcgGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCt                                                                                                                                                                                                        <  1:588527/250‑1 (MQ=255)
                                                        acgggtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCATCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGc                                                                                                                                                                                                <  2:306965/251‑1 (MQ=255)
                                                            gtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCg                                                                                                                                                                                            <  1:131997/251‑1 (MQ=255)
                                                                  tCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTATG‑T‑TTTC‑CTTTTTTTTTTTTTCGAGCTTTCGGAATGCGGTGAATA‑AGG‑TCACGCTTTTTCGTCCCGCTCCCGAAGCGTCg‑ccgc                                                                                                                                                                                       >  1:191188/1‑246 (MQ=255)
                                                                          ggACGCATTCATATTGACGACGAAGCCAGCAGGGGCAAGCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCgg                                                                                                                                                                                                                                              <  2:636392/189‑1 (MQ=255)
                                                                                aTTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTATG‑T‑TTTC‑CTTTTTTTTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TGACGCTTTTTTGTCCCCCTCCCGTAGCGTCG‑CAGCCGGAATATCAGGCCTc                                                                                                                                                                       >  1:382111/1‑251 (MQ=255)
                                                                                aTTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAG‑G‑TTGC‑CTTTTTTTTTTTTTCGATCTTCCGGAATGCGGGGGATA‑CGG‑TCACGCTTTTTTGTCCCGCTCCAGTAGCGGCG‑CAGCAGGAATATCAGGCCTc                                                                                                                                                                       >  1:329592/1‑251 (MQ=255)
                                                                                   aatatTGACGACGAAGCCAGACGGGGCAATCCCTCCCAGTCCTAGCATTACCAGGCCGATCCCGTTTTGTCCGTCGTTCGCGCCGGGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTGTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTGCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGTATATCAGGCCTCCCg                                                                                                                                                                    <  2:416305/250‑1 (MQ=255)
                                                                                      atTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACA‑G‑TTGC‑CTTTTTTTTTTTTTTGTTCTTCCGGAATACGGGGGATA‑CGGGCCAC‑CTTTTTCGTCCCCCTCCAGGAGCGGCG‑CAGCAGGGATATCAgacatcccccc                                                                                                                                                                  >  1:233949/1‑240 (MQ=255)
                                                                                                         ggggCAACCCCCCCCAGTACCAGCATTACCAGGCCGAGCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAgg                                                                                                                                              <  2:329042/251‑1 (MQ=255)
                                                                                                             cAATCCCCACCAGTACCAGCATTACCAGGCCGACCCCTGTTTGGCCGTCGTTCGCGCCGTGCGAAAACGCCACGGCCGCAGCGGAAACAATCAGCGCAATTCGCGGCCATAATGGCGGTTTCCG‑T‑TTGC‑TGTTTTTCTTTTTGCGAGCTTCCGGAATGCGGTGAATA‑CGG‑TCCCGCTTTTTCGTCCCGCTCCAGTAGCGGCG‑CAGCAGGAATATCAGGCCTCTCGCAATGACCAGGCCGACGATAGGGGaa                                                                                                                                          <  1:357036/251‑1 (MQ=255)
                                                                                                                         gTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTAG‑T‑TTGC‑CTTTTTTCTTTTTGCGCTCTTTCGGGATGCGGTGGATA‑CGG‑TTACGCTTTTTTGTCCCCCCACAGTCGCGTCG‑CAGCAGGAATATCTGGGCCCCCCCAATTACCCAGCCGACGATAGGGGGAA‑Caaaaaaag                                                                                                                               >  2:452531/1‑243 (MQ=255)
                                                                                                                          tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGGT‑TTGC‑‑TTTTTTTTTTTTGCGATCTTCCTGAATGAGGGGAATA‑CGG‑TCACGCTTTTTTGTCCCCCTCCCCGAGCGTCG‑CAGCAGGGAAATCAAGGCTCCCCCAAAGACCCGGCCGACGATAGGGGAAAACAA‑CAGCgagg                                                                                                                            >  1:428004/1‑251 (MQ=255)
                                                                                                                          tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTTG‑CTTTTTTTTTTTTTTGGTCTTTCGGGATGCGGGGAATA‑CGG‑TCAATCTTTTTCGTCCCGCCCCCGTAGCGTCG‑CCGCAGGGATAATAGGCCGCCCGCCATGACCAGGCCGACCGCAGGGGAAA‑AAAACAGCgagg                                                                                                                            >  1:279067/1‑251 (MQ=255)
                                                                                                                                      cAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTCG‑T‑TTGC‑CTTTTTTTTTTTTTCGATCTTCCGGGATGCGGTGGATA‑CGG‑TCACGCTTTTTCGTCCCGCCCCCGTAGCGTCG‑CAGCAGGGATATCAGGCCCCCCGCCATGACCCGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTgg                                                                                                                >  1:118958/1‑251 (MQ=255)
                                                                                                                                       aGGCCTATCCCTTTTTGGCCGCCTTCCGCGCCTTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGt                                                                                                               <  2:118958/251‑1 (MQ=255)
                                                                                                                                                 cTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACA‑T‑TTTC‑CCTTTTTCTTTTTTAGATCTTCCGGAAAGCGGTGAAAAACGG‑‑CACGCTTTTTTGTCCCCCTCCAGTAGCGTCGGCGGCAG‑AATATCAGGACCCCCCCAATGACCACGCCGACGCCACGGGGAA‑CAAAAAACGAGGAGAAAATTTTTGTccactcacgc                                                                                                     >  1:316804/1‑242 (MQ=255)
                                                                                                                                                   ttttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCa                                                                                                    <  1:233536/250‑1 (MQ=255)
                                                                                                                                                    tttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑G‑TTTC‑CTTTTTTTTTTTTGCGGTCTTCCGGAATGCGGTGGATA‑CGG‑TCACGCTTTTTCGTCCCCCCCCAGGAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCCAAGACCCGGCCGGCGATTGGGGGAA‑CAAACAGCGAGGGGAAAAATTTTGTCGCTTCACGCAgg                                                                                                  >  1:620486/1‑251 (MQ=255)
                                                                                                                                                     ttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTTTTTTTTTGATCTTCCGGGATGCGGTGGATa                                                                                                                                                                                                                                 >  2:255574/1‑127 (MQ=255)
                                                                                                                                                     ttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATa                                                                                                                                                                                                                                 <  1:255574/127‑1 (MQ=255)
                                                                                                                                                      tGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGt                                                                                                 <  1:554203/250‑1 (MQ=255)
                                                                                                                                                                          cGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATg                                                                            <  2:329592/251‑1 (MQ=255)
                                                                                                                                                                            aaaaCGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGAt                                                                             <  1:165479/248‑1 (MQ=255)
                                                                                                                                                                               acgccacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAG‑T‑TTGC‑CTTTTTTTTTTTTTCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CCGCAGGAAAATTAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGATTAAAGCATCCATCACCggtggaccg                                                                       >  1:44031/1‑243 (MQ=255)
                                                                                                                                                                                   cacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTg                                                                                                                 <  2:596083/205‑1 (MQ=255)
                                                                                                                                                                                           cAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTTCGATCTTCCGGAATGCGGGGGATA‑AGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGCTAGGGGGAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTTACGCATCCATCACCGAAGAGCCCGGTAACACcgcg                                                           >  1:408853/1‑251 (MQ=255)
                                                                                                                                                                                               ggAAACAATCAGCGCAATACGCGTCCAGAAAGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGCCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTgg                                                       <  2:335513/251‑1 (MQ=255)
                                                                                                                                                                                                    cAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTaaa                                                  <  2:620486/251‑1 (MQ=255)
                                                                                                                                                                                                        cAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGa                                              <  1:247637/251‑1 (MQ=255)
                                                                                                                                                                                                           cgcAATACGCGTCCAGAATGGCGGTTTTCG‑T‑TTTC‑CTTTTTTCTTTTTGCGCTCTTCCGGAATGCGGGGGATA‑CGG‑TCACGCTTTTTTGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGGCTCCCGCCATGAACAGGCCGACGATAGGGGGAA‑CAATCAGCGGGGAGGAAAATTTGGGCACTTCACGCAGGTTTAACGCATCCCTCTACGATGAGCCGGTTAACAGCGCGTTGGGTAAACCCATGc                                           >  1:462639/1‑251 (MQ=255)
                                                                                                                                                                                                                   gcgTCCAGAAGGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATc                                   <  1:103599/251‑1 (MQ=255)
                                                                                                                                                                                                                         aGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTGCGTCCCGCGCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCa                             <  2:408853/251‑1 (MQ=255)
                                                                                                                                                                                                                             tGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCg                                                                               >  1:121069/1‑197 (MQ=255)
                                                                                                                                                                                                                             tGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGGTTAACGCATCCATCACCg                                                                               <  2:121069/197‑1 (MQ=255)
                                                                                                                                                                                                                              ggCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCa                         >  1:142868/1‑250 (MQ=255)
                                                                                                                                                                                                                               gCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAg                       <  1:452531/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGgt                     >  1:515063/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cG‑T‑TTGC‑CTTTTTTCTTTTGGCGATCTACCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCg                <  1:176938/250‑1 (MQ=255)
                                                                                                                                                                                                                                          t‑ttGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAAc             >  1:351703/1‑251 (MQ=255)
                                                                                                                                                                                                                                               c‑cTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTgg          >  1:67330/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   tttttCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAACATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATTAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACCGGAGGc      >  1:555938/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        cTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGAAGAGCCGGTTAACAGCGCGTTGGTTAAACCGCTGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGt  >  2:589771/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
GGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGT  >  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: ‑