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→CGApitBmetal 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.

GCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTA  >  NC_000913/3135543‑3136025
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          
gCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGa                                                                                                                                                                                                                                              <  1:308900/251‑1 (MQ=255)
 cTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGATTACCGTTGCCCTTTTTTTTTTTGCgtgattccggg                                                                                                                                                                                                                                                >  1:38355/1‑240 (MQ=255)
         ggTAGAGTTCGCAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGAATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAACACGCCACGCCCGCCGCGGAAACAATCAGCGCACTACGCGTCCAGAATGGCGGTTTATGGTTTCAATTTTTCTTTTTTCATTCTTCCggcctgcggggg                                                                                                                                                                                                                                      >  2:306944/1‑241 (MQ=255)
            aGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGGTCTTCCGGAAAGCgggggat                                                                                                                                                                                                                                   >  1:533661/1‑245 (MQ=255)
                          gTAACGGCATCGCGGGTTCGGGAAATTTCAGAGCCGGACGCATTCATAGTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTGGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCttttt                                                                                                                                                                                                                  <  2:105215/251‑1 (MQ=255)
                             aCGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGGTTTCGTTTGCCTTTTTTCTTTTTGCGGTCTTCCGGGATGCGGGGAATaagggtacccttcttcgt                                                                                                                                                                                                                 >  2:93251/1‑234 (MQ=255)
                                 cATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAACGTTTCCTTTTTTTTTTTTGAGAACTTCCGGAATGCGGTGAATAAGGTCACGGTTTTTTGTcccc                                                                                                                                                                                                           >  1:390767/1‑250 (MQ=255)
                                       gggtacgggtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTTGATTTTCCGGGAGTCGGTGAATACGGTACAGCTTTTTTGTCCCGCTcccg                                                                                                                                                                                                     >  1:93294/1‑249 (MQ=255)
                                              ggtaATTTCATAGCCTGACTCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTc                                                                                                                                                                                              <  2:143232/251‑1 (MQ=255)
                                               gtaATTTCATAGCCGGACGCATTCAGAGTGACGACGAAGCCAGCAGGGGCAAGCCCCTCCAGTACCAGCATTACCAGGCCGATCCCTGTTGGTCCGTCGTGCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCTCAATACGCGTCCAGAATGGCGGTTTACGTTGGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGGGAATACGGTCACGCTTTGTCGTCCCGCTCCAGTAGCGTCg                                                                                                                                                                                             <  2:136894/251‑1 (MQ=255)
                                                 aaTTTCATAGCCGGACTCATTCATATTGACGACGAAGCCAGCAGGGGCAAGCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCATCCCGCTCCAGTAGCGTCgca                                                                                                                                                                                           <  2:425773/251‑1 (MQ=255)
                                                  aTTTCATAGCCGGACGCATTCATATTGACGGCGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCgcag                                                                                                                                                                                          <  1:38070/251‑1 (MQ=255)
                                                    ttCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGGATGCGGGGGATACGGGCACGCTTTTTTGTCCCGCTCCAGTAGCGTCgcagcc                                                                                                                                                                                        >  1:160868/1‑250 (MQ=255)
                                                      cATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAgg                                                                                                                                                                                      <  1:486652/251‑1 (MQ=255)
                                                            cGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACCTTTTCCTTTTTTTTTTTTGCGATCTTCCCGGATGCGGTGAATACGGTCACGCTTTTTTGTCCCGCCCCCGTAGCGTCGCAGCCGGAAAATc                                                                                                                                                                                >  1:184005/1‑251 (MQ=255)
                                                             ggACGCATTCATATTGACGACGAAGCCAGCAGGGGCAAGCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCa                                                                                                                                                                               <  2:75763/251‑1 (MQ=255)
                                                               aCGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAgg                                                                                                                                                                             <  2:160868/251‑1 (MQ=255)
                                                                 gCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGGTTGCCTTTTTTTTTTTTGCGATTTTCCGGAATGCGGGGAATAAGGTCACGCTTTTTCGTCCCGCTCCAGTCGCGTCGCAGCAGGAATATccagcc                                                                                                                                                                           >  1:476079/1‑246 (MQ=255)
                                                                         aatgaagacgaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGGATGACTTTTTTCTTTTTTCGACCCTCCGGGAAAGGGTGAATACGGGCAAGTTATTTAGGCCCGCTCCAGTAGCGTCTCAGCAGGAATAACAGGGCTCCCGCa                                                                                                                                                                    >  2:150729/7‑250 (MQ=255)
                                                                          ttGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGGTTTCGGTTGGCCTTTTTTTTTTTTTGGACTTCCCGGAAGCGGTGGAAGCGGGCACCCTTTTTCGTCCCGCTCCCGGAGCGTCGCCGCAGGAATATCAGGccccccccc                                                                                                                                                                    >  1:493595/1‑246 (MQ=255)
                                                                           tgccgacgaAGCCAGCAGGGGCAATCCCCTCCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAAt                                                                                                                                                                 <  2:138149/247‑1 (MQ=255)
                                                                             atcgcgaAGCCCGCAGCTGCAATCCCCACCATTTCCAGCATTACCAGGCCGAGCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCATCGGAAACAATCATCTCAATACGCGTCCAGAATTGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCTTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGAc                                                                                                                                                              <  2:295342/247‑1 (MQ=255)
                                                                                acgaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCCTTTTTTTTTTTTCGGTCTTCCGGGATGCGGGGGAAAAGGTCACGCTTTTTTGTCCCCCTCCCGGAGGGTCGCCGCAgggaaaacaggccccccccacacgac                                                                                                                                                              >  1:338362/1‑225 (MQ=255)
                                                                                 cgaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATACCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGGTTAAGGTTGGCTTTTTTTTTTTTGCGGTCTTTCGGAATGCGGGGCAAACGGTCACCCGTTTTTGTCCCCCTCCAGGAGCGTCGCAGGAGGGATATCAGGGCTCCCCCaaaggccaa                                                                                                                                                            >  1:295342/1‑243 (MQ=255)
                                                                                    aGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGAATGCGGGGGATACGGTCACGCTTTTTCGTCCCCCTCCAGTCGCGTCGCAGCAGGGAAATCAAGCCTccccccaagaaccggccc                                                                                                                                                        >  1:173866/1‑236 (MQ=255)
                                                                                           agGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTCGATCTTCCGGGATGCGGGGAATACGGTCACGCTTTTTCGTCCCGCTCCCCTAGCGGCGCAGCAGGGATATCAGGGCCCCCCCAATGAcccagcccacggt                                                                                                                                                  >  1:484540/1‑238 (MQ=255)
                                                                                             gggCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCCGGCCTCCCGCAATGACCCGGCCGCCGATAggg                                                                                                                                              >  1:96853/1‑251 (MQ=255)
                                                                                                               ccAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCgaggag                                                                                                                            <  1:477420/251‑1 (MQ=255)
                                                                                                                      tACCAGGCCGATCCCGTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGaaagtt                                                                                                                      <  1:162324/250‑4 (MQ=255)
                                                                                                                       aCCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGAATGCGGGGAATACGGGCACGCCTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGGATATCAGGCCTCCCGCCAAGACCCGGCCGACGGGAGGGGAAACAAACACCGAGGAGAAAAattt                                                                                                                    >  1:40734/1‑251 (MQ=255)
                                                                                                                       aCCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAAtttt                                                                                                                    <  1:37611/251‑1 (MQ=255)
                                                                                                                                aTCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGGGAATACGGTCACGCTTTTTTGGCCCCCCCCCGTAGCGTCGCCGCAGGGAAATCAGGCCCCCCGCAAATGCCCGGCCGACGGTAGGGGAAACAATCCACGCGGGGAAAATTTTGGTCACTTc                                                                                                          >  1:182079/1‑251 (MQ=255)
                                                                                                                                          gTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGttt                                                                                                <  1:77034/251‑1 (MQ=255)
                                                                                                                                             cgtcgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGAATGCGGGGGATACGGTCACGCTTTTTTGTCCCGCCCCAGTAGCGTCGCAGCAGGGATATCAGGCCTCCCGCCATGACCAGGCCGACGAAAGGGGAAACAATCAGCGGGGGGAAAATTTTGGGCACTTCACGCAGGGTTAAc                                                                                             >  1:233317/1‑251 (MQ=255)
                                                                                                                                              gtcgtTCGCGCCGTGCGAAAACGCCACGCCCGCATCTGAAACAATCAGCGCAATACGCGTCCAGAATTGCGGTTTACGTTTGCCTTTTTTCGTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACg                                                                                            <  1:60045/251‑1 (MQ=255)
                                                                                                                                              gtcggtcgCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAAGCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGACTGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTGCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACg                                                                                            <  2:395741/251‑1 (MQ=255)
                                                                                                                                                cgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGca                                                                                          <  1:401812/251‑1 (MQ=255)
                                                                                                                                                cgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGc                                                                                           <  1:196448/250‑1 (MQ=255)
                                                                                                                                                 gtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTTCCCCCTTTTTGTTCCGCTCCCGTAGCGCCGCCGCAGGAAAATCAGGGCTCCCCCCATGACCCGGCCGACGATAGGGGAAACAAATAGCGAGGGGAAAATTTTgggcaactctcgccggctttaaacac                                                                                          >  1:64369/1‑227 (MQ=255)
                                                                                                                                                 gtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGcat                                                                                         <  1:489522/251‑1 (MQ=255)
                                                                                                                                                 gtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGcat                                                                                         <  2:230021/251‑1 (MQ=255)
                                                                                                                                                         cGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTCGATCTTCCGGGATGCGGTGGATACGGTCCCGCTTTTTCTTACCGCCCCACAAGCGGCGCCGCCGGAATATTAGGCCTCCCCCCATGACCCGGCCGACGATAGGGGGAACAAACAGCGAGGCGAAAATTTTGGGCACCTCACGCAGGGTTAACGcagccaacccc                                                                                 >  2:383133/1‑246 (MQ=255)
                                                                                                                                                                    gccacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAGTTTGCCTTTTTTCTTTTTTAGATCTTTCGGAATGCGGTGAATACGGGCACGCTTTTTCGTCCCCCCCCAGTAGCGGCGCAGCAGGAATAAAAGGCCTCCCGCAATGACCCGGCCGACGGTAGGGGAAAAAATCAGCGAGGAGAAAATTTTGGTCACCTCACGCAGGTTTAAAGCATCCATCACCGATGAGCCggg                                                                      >  1:483047/1‑250 (MQ=255)
                                                                                                                                                                     ccacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGTTTTGCGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTGTAACGCATCCATCGCCGATGAGCCGGtt                                                                     <  2:503553/251‑1 (MQ=255)
                                                                                                                                                                              cAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGGGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAAAATTAGGCCTCCCCCAATGACCAGGCCGACGAAAGGGGAAACAAACAGCGGGGAGAAAAATTTGGGCACTTCACGCAGGGTTAACGCATCCATCACCGCTGAACCGGTTTACAAcgcg                                                            >  1:173621/1‑251 (MQ=255)
                                                                                                                                                                                 cGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTg                                                         <  1:410543/251‑1 (MQ=255)
                                                                                                                                                                                 cGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAGGTTGCCTTTTTTCTTTTTTCGATCTTTCGGGATGCGGTGAATAAGGGCACGCTTTTTTGTCCCCCCCCCTTAACGGCGCAGCAGGAATATCAGGCCTCCCGCAATGAACAGGGCGCCGGGAGGGGGAACAAAAAGCGGGGGGAAAAATTTTGTGACTTCTCGCAGGTTTAAAACATCCATCACCGCTGCGCCGGTTataacggcggtg                                                         >  2:168359/1‑240 (MQ=255)
                                                                                                                                                                                    aaaCAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGtt                                                      <  1:481596/251‑1 (MQ=255)
                                                                                                                                                                                                             gAATGGCGGTTTTCGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCCGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCCACGATAGGGGAAAAAAACAGCGAGGAGAAAATTTTGGTCCACTTACGCAGGTTTAAACCCTCCATCACCGATTAGCCGGTTAACAGCGCGGTGGTTAAAACGACGcccgaaatcgccccac                             >  2:498237/1‑237 (MQ=255)
                                                                                                                                                                                                             gAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCaa                             <  1:168359/251‑1 (MQ=255)
                                                                                                                                                                                                              aaTGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAAt                            <  2:533661/251‑1 (MQ=255)
                                                                                                                                                                                                                          aCGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCg                 <  1:453172/251‑1 (MQ=255)
                                                                                                                                                                                                                               ttGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAAc              <  2:184005/250‑1 (MQ=255)
                                                                                                                                                                                                                                tGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACt             <  1:93251/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTgg           >  1:244298/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tttttCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGGCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAg         <  2:517154/249‑1 (MQ=255)
                                                                                                                                                                                                                                           tttttGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTa  >  1:358629/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          
GCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTA  >  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: