Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,7810GT29.8% 75.0 / 16.4 57R197R (CGC→CGA) pitBmetal phosphate:H(+) symporter PitB
Reads supporting (aligned to +/- strand):  ref base G (11/29);  new base T (16/1);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 3.72e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTA  >  NC_000913/3135543‑3136025
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          
gCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑a                                                                                                                                                                                                                                              <  1:308900/251‑1 (MQ=255)
 cTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑ATTA‑CCGTTGCCCTTTTTTTTTTTGCgtgattccggg                                                                                                                                                                                                                                                >  1:38355/1‑240 (MQ=255)
         ggTAGAGTTCGCAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGAATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAACACGCCACGCCCGCCGCGGAAACAATCAGCGCACTACGCGTCCAGAATGGCGG‑TTTA‑TGGTTTCAATTTTTCTTTTTTCATTCTTCC‑ggcctgcggggg                                                                                                                                                                                                                                      >  2:306944/1‑241 (MQ=255)
            aGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTTCGGTCTTCC‑GG‑AAAGCgggggat                                                                                                                                                                                                                                   >  1:533661/1‑245 (MQ=255)
                          gTAACGGCATCGCGGGTTCGGGAAATTTCAGAGCCGGACGCATTCATAGTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTGGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCttttt                                                                                                                                                                                                                  <  2:105215/251‑1 (MQ=255)
                             aCGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑GTTT‑CGTTTGCCTTTTTTCTTTTTGCGGTCTTCC‑GG‑GATGCGGGGAATaagggtacccttcttcgt                                                                                                                                                                                                                 >  2:93251/1‑234 (MQ=255)
                                 cATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTAACGTTT‑CCTTTTTTTTTTTTGAGAACTTCC‑GG‑AATGCGGTGAATAAGGT‑CAC‑GGTTTTTTGTcccc                                                                                                                                                                                                           >  1:390767/1‑250 (MQ=255)
                                       gggtacgggtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTTTTTTTTTGATTTTCC‑GG‑GAGTCGGTGAATACGGT‑‑ACAGCTTTTTTGTCCCGCTcccg                                                                                                                                                                                                     >  1:93294/1‑249 (MQ=255)
                                              ggtaATTTCATAGCCTGACTCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTc                                                                                                                                                                                              <  2:143232/251‑1 (MQ=255)
                                               gtaATTTCATAGCCGGACGCATTCAGAGTGACGACGAAGCCAGCAGGGGCAAGCCCCTCCAGTACCAGCATTACCAGGCCGATCCCTGTTGGTCCGTCGTGCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCTCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTGGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGGGAATACGGT‑CAC‑GCTTTGTCGTCCCGCTCCAGTAGCGTCg                                                                                                                                                                                             <  2:136894/251‑1 (MQ=255)
                                                 aaTTTCATAGCCGGACTCATTCATATTGACGACGAAGCCAGCAGGGGCAAGCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCATCCCGCTCCAGTAGCGTCgca                                                                                                                                                                                           <  2:425773/251‑1 (MQ=255)
                                                  aTTTCATAGCCGGACGCATTCATATTGACGGCGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCgcag                                                                                                                                                                                          <  1:38070/251‑1 (MQ=255)
                                                    ttCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTTCGATCTTCC‑GG‑GATGCGGGGGATACGGG‑CAC‑GCTTTTTTGTCCCGCTCCAGTAGCGTCgcagcc                                                                                                                                                                                        >  1:160868/1‑250 (MQ=255)
                                                      cATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAgg                                                                                                                                                                                      <  1:486652/251‑1 (MQ=255)
                                                            cGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CCTTTTCCTTTTTTTTTTTTGCGATCTTCC‑CG‑GATGCGGTGAATACGGT‑CAC‑GCTTTTTTGTCCCGCCCCCGTAGCGTCGCAGCCGGAAAATc                                                                                                                                                                                >  1:184005/1‑251 (MQ=255)
                                                             ggACGCATTCATATTGACGACGAAGCCAGCAGGGGCAAGCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCa                                                                                                                                                                               <  2:75763/251‑1 (MQ=255)
                                                               aCGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAgg                                                                                                                                                                             <  2:160868/251‑1 (MQ=255)
                                                                 gCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGGTTGCCTTTTTTTTTTTTGCGATTTTCC‑GG‑AATGCGGGGAATAAGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTCGCGTCGCAGCAGGAATATccagcc                                                                                                                                                                           >  1:476079/1‑246 (MQ=255)
                                                                         aatgaagacgaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGGATGACTTTTTTCTTTTTTCGACCCTCC‑GGGAAAG‑GGTGAATACGGG‑CAA‑GTTATTTAGGCCCGCTCCAGTAGCGTCTCAGCAGGAATAACAGGGCTCCCGCa                                                                                                                                                                    >  2:150729/7‑250 (MQ=255)
                                                                          ttGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGGTTTC‑GGTTGGCCTTTTTTTTTTTTTGGA‑CTTCCCGG‑AA‑GCGGTGGAAGCGGG‑CAC‑CCTTTTTCGTCCCGCTCCCGGAGCGTCGCCGCAGGAATATCAGGccccccccc                                                                                                                                                                    >  1:493595/1‑246 (MQ=255)
                                                                           tgccgacgaAGCCAGCAGGGGCAATCCCCTCCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑t                                                                                                                                                                 <  2:138149/247‑1 (MQ=255)
                                                                             atcgcgaAGCCCGCAGCTGCAATCCCCACCATTTCCAGCATTACCAGGCCGAGCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCATCGGAAACAATCATCTCAATACGCGTCCAGAATTGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCTTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGAc                                                                                                                                                              <  2:295342/247‑1 (MQ=255)
                                                                                acgaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCCTTTTTTTTTTTTCGGTCTTCC‑GG‑GATGCGGGGGAAAAGGT‑CAC‑GCTTTTTTGTCCCCCTCCCGGAGGGTCGCCGCAgggaaaacaggccccccccacacgac                                                                                                                                                              >  1:338362/1‑225 (MQ=255)
                                                                                 cgaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATACCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑GTTA‑AGGTTGGCTTTTTTTTTTTTGCGGTCTTTC‑GG‑AATGCGGGGCAAACGGT‑CAC‑CCGTTTTTGTCCCCCTCCAGGAGCGTCGCAGGAGGGATATCAGGGCTCCCCCaaaggccaa                                                                                                                                                            >  1:295342/1‑243 (MQ=255)
                                                                                    aGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTTCGATCTTCC‑GG‑AATGCGGGGGATACGGT‑CAC‑GCTTTTTCGTCCCCCTCCAGTCGCGTCGCAGCAGGGAAATCAAGCCTccccccaagaaccggccc                                                                                                                                                        >  1:173866/1‑236 (MQ=255)
                                                                                           agGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTTTTTTTTCGATCTTCC‑GG‑GATGCGGGGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCCCTAGCGGCGCAGCAGGGATATCAGGGCCCCCCCAA‑TGAcccagcccacggt                                                                                                                                                  >  1:484540/1‑238 (MQ=255)
                                                                                             gggCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTTCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCCGGCCTCCCGCAA‑TGACCCGGCCGCCGATAggg                                                                                                                                              >  1:96853/1‑251 (MQ=255)
                                                                                                               ccAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCgaggag                                                                                                                            <  1:477420/251‑1 (MQ=255)
                                                                                                                      tACCAGGCCGATCCCGTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGaaagtt                                                                                                                      <  1:162324/250‑4 (MQ=255)
                                                                                                                       aCCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTTCGATCTTCC‑GG‑AATGCGGGGAATACGGG‑CAC‑GCCTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGGATATCAGGCCTCCCGCCA‑AGACCCGGCCGACGGGAGGGGAAACAAACACCGAGGAGAAAAattt                                                                                                                    >  1:40734/1‑251 (MQ=255)
                                                                                                                       aCCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAAtttt                                                                                                                    <  1:37611/251‑1 (MQ=255)
                                                                                                                                aTCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGGGAATACGGT‑CAC‑GCTTTTTTGGCCCCCCCCCGTAGCGTCGCCGCAGGGAAATCAGGCCCCCCGCAAATG‑CCCGGCCGACGGTAGGGGAAACAATCCACGCGGGGAAAATTTTGGTC‑ACTTc                                                                                                          >  1:182079/1‑251 (MQ=255)
                                                                                                                                          gTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGttt                                                                                                <  1:77034/251‑1 (MQ=255)
                                                                                                                                             cgtcgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTTCGATCTTCC‑GG‑AATGCGGGGGATACGGT‑CAC‑GCTTTTTTGTCCCGCCCCAGTAGCGTCGCAGCAGGGATATCAGGCCTCCCGCCA‑TGACCAGGCCGACGAAAGGGGAAACAATCAGCGGGGGGAAAATTTTGGGC‑ACTTCACGCAGGGTTAAc                                                                                             >  1:233317/1‑251 (MQ=255)
                                                                                                                                              gtcgtTCGCGCCGTGCGAAAACGCCACGCCCGCATCTGAAACAATCAGCGCAATACGCGTCCAGAATTGCGG‑TTTA‑CGTTTGCCTTTTTTCGTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACg                                                                                            <  1:60045/251‑1 (MQ=255)
                                                                                                                                              gtcggtcgCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAAGCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑ACTGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTGCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACg                                                                                            <  2:395741/251‑1 (MQ=255)
                                                                                                                                                cgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGca                                                                                          <  1:401812/251‑1 (MQ=255)
                                                                                                                                                cgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGc                                                                                           <  1:196448/250‑1 (MQ=255)
                                                                                                                                                 gtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGTTCCC‑CCTTTTT‑GTTCCGCTCCCGTAGCGCCGCCGCAGGAAAATCAGGGCTCCCCCCA‑TGACCCGGCCGACGATAGGGGAAACAAATAGCGAGGGGAAAATTTTgggcaactctcgccggctttaaacac                                                                                          >  1:64369/1‑227 (MQ=255)
                                                                                                                                                 gtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGcat                                                                                         <  1:489522/251‑1 (MQ=255)
                                                                                                                                                 gtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGcat                                                                                         <  2:230021/251‑1 (MQ=255)
                                                                                                                                                         cGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTTTTTTTTCGATCTTCC‑GG‑GATGCGGTGGATACGGT‑CCC‑GCTTTTTCTTACCGCCCCACAAGCGGCGCCGCCGGAATATTAGGCCTCCCCCCA‑TGACCCGGCCGACGATAGGGGGAACAAACAGCGAGGCGAAAATTTTGGGC‑ACCTCACGCAGGGTTAACGcagccaacccc                                                                                 >  2:383133/1‑246 (MQ=255)
                                                                                                                                                                    gccacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑AGTTTGCCTTTTTTCTTTTTTAGATCTTTC‑GG‑AATGCGGTGAATACGGG‑CAC‑GCTTTTTCGTCCCCCCCCAGTAGCGGCGCAGCAGGAATAAAAGGCCTCCCGCAA‑TGACCCGGCCGACGGTAGGGGAAAAAATCAGCGAGGAGAAAATTTTGGTC‑ACCTCACGCAGGTTTAAAGCATCCATCACCGATGAGCCggg                                                                      >  1:483047/1‑250 (MQ=255)
                                                                                                                                                                     ccacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGT‑TTTG‑CGTTTGCCTTTTTTCTTTTTTCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTGTAACGCATCCATCGCCGATGAGCCGGtt                                                                     <  2:503553/251‑1 (MQ=255)
                                                                                                                                                                              cAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGGGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAAAATTAGGCCTCCCCCAA‑TGACCAGGCCGACGAAAGGGGAAACAAACAGCGGGGAGAAAAATTTGGGC‑ACTTCACGCAGGGTTAACGCATCCATCACCGCTGAACCGGTTTACAAcgcg                                                            >  1:173621/1‑251 (MQ=255)
                                                                                                                                                                                 cGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTg                                                         <  1:410543/251‑1 (MQ=255)
                                                                                                                                                                                 cGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑AGGTTGCCTTTTTTCTTTTTTCGATCTTTC‑GG‑GATGCGGTGAATAAGGG‑CAC‑GCTTTTTTGTCCCCCCCCCTTAACGGCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGAACAGGGCGCCGGGAGGGGGAACAAAAAGCGGGGGGAAAAATTTTGTG‑ACTTCTCGCAGGTTTAAAACATCCATCACCGCTGCGCCGGTTataacggcggtg                                                         >  2:168359/1‑240 (MQ=255)
                                                                                                                                                                                    aaaCAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGtt                                                      <  1:481596/251‑1 (MQ=255)
                                                                                                                                                                                                             gAATGGCGG‑TTTT‑CGTTTGCCTTTTTTCTTTTTTCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCCGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCCACGATAGGGGAAAAAAACAGCGAGGAGAAAATTTTGGTCCACTT‑ACGCAGGTTTAAACCCTCCATCACCGATTAGCCGGTTAACAGCGCGGTGGTTAAAACGACGcccgaaatcgccccac                             >  2:498237/1‑237 (MQ=255)
                                                                                                                                                                                                             gAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCaa                             <  1:168359/251‑1 (MQ=255)
                                                                                                                                                                                                              aaTGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAAt                            <  2:533661/251‑1 (MQ=255)
                                                                                                                                                                                                                          a‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCg                 <  1:453172/251‑1 (MQ=255)
                                                                                                                                                                                                                               ttGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAAGACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAAc              <  2:184005/250‑1 (MQ=255)
                                                                                                                                                                                                                                tGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACt             <  1:93251/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTgg           >  1:244298/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tttttCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGGCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAg         <  2:517154/249‑1 (MQ=255)
                                                                                                                                                                                                                                           tttttGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTa  >  1:358629/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          
GCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGG‑TTTA‑CGTTTGCCTTTTTTCTTTTTGCGATCTTCC‑GG‑AATGCGGTGAATACGGT‑CAC‑GCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAA‑TGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTC‑ACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTA  >  NC_000913/3135543‑3136025

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