Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,844	0	A	C	47.5%	80.4 / 17.7	60	L728V (TTA→GTA)	arcB	aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (26/5);  new base C (0/28);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 6.86e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGG  >  NC_000913/3350703‑3350991                                                                                                                                                    |                                                                                                                                                     cTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCGGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTg                                                                                                                                             >  1:791962/1‑150 (MQ=255)     tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGAcccc                                                                                                                                        >  1:1135008/1‑151 (MQ=255)      gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCccg                                                                                                                                       <  1:99141/151‑1 (MQ=255)       ccacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCccgc                                                                                                                                      <  2:603967/151‑1 (MQ=255)               gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCtt                                                                                                                              >  2:1056818/1‑151 (MQ=255)                cTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGCCCCCGCCGCAc                                                                                                                                 <  1:655666/147‑1 (MQ=255)                cTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCttt                                                                                                                             >  1:168484/1‑151 (MQ=255)                 tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCttt                                                                                                                             >  1:316404/1‑150 (MQ=255)                       cccttctcctccgtGCCCCCACTCTTTTTTCACCTCTCCAATCCATTCCCCGACGTTC‑TCTCCCCAGCCC‑GAAGGGTCGGGAGAC‑TGATTTTGCTG‑CCCCGCTTTTTGAAG‑AGGG‑CGTAC‑CCCCACCGCCCCCGCCGCACCTTTAAtttt                                                                                                                       <  1:175953/139‑1 (MQ=255)                         cTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATg                                                                                                                    >  2:984717/1‑151 (MQ=255)                             tccgcgttaccGCCTCTTTTTTTACATCTCTTAAATCCATTCCCCGCCGTTA‑TCTCCCCGGGCC‑GGAAGCTCGGGCGAC‑TGATTTTGCTG‑CCCCCGTTTTTGCAA‑GTGC‑GGTAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGt                                                                                                                   <  2:92574/138‑1 (MQ=255)                                  cgtGCCCCCCCTCTTCTTTCATCTCTACAACCCATTCCCCGCCTTTA‑TCTTCCCAGGCC‑GAAAGGTCGGGAGAC‑TGAATTTGCTA‑CCCCAGTTGTTGCAG‑ATGG‑CGTAC‑CCCCTCTACCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                               <  2:1207589/147‑1 (MQ=255)                                    tGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                          >  2:430292/1‑150 (MQ=255)                                     gccGCCTCTCTTCTTTCTTCTC‑TCAAACCCTTCCCCCGCCTTATTCTCCCCAGGCC‑GGAGGTTCAGAGAAC‑TG‑ATTTGGTGTCCCCAGTTTTTGCAG‑TGGG‑CGTAC‑CCCCACGCCCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                               <  2:646436/142‑1 (MQ=255)                                        ccccctctTCTTTCACCTCTTCAACCCATCCCCCGACTTTA‑TCTCCCAGGGCC‑GGAAGGTCAGGAGAC‑TGAATTGGCGG‑ACCCGGTTTTTGCAG‑ATGG‑CGTAC‑CCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                                      <  2:479819/146‑1 (MQ=255)                                               ttttttCATCTCTCCAATCCATCCCCCGACTTTT‑TCTCCCCAGGCG‑GAAAGGTCAGGAGAT‑TAATTTTGCTG‑CCCCATCTGTGGCAG‑AGGG‑CGTAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                               <  2:508674/147‑1 (MQ=255)                                               ttcttcCTTCTCTCCAACCCATTCACCGACGTTA‑TCTCCCCGGGCG‑GAAGGGTCGGGAGAC‑TGATTTTGCGG‑CCCCAGTTGTTGCAG‑ATGG‑CGTCC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt                                                                                                 <  2:412512/148‑1 (MQ=255)                                               ttcttcCCTCTCTCCACTCCCTCCCCCGGCTTTA‑TCTCCCCGGGCC‑GGAAGGTCAGGACAC‑TGATTTTGCTC‑ACCCATTTTTTG‑AG‑GTGGCCTTAC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt                                                                                                 <  1:250866/148‑1 (MQ=255)                                                tatttCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCCGGAGAC‑TGATTTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑ACCCACTGAACGCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                             >  1:1063352/3‑151 (MQ=255)                                                 cttTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                            >  2:648779/1‑151 (MQ=255)                                                  tttCATCTCTTCAACCCATCCCCCGACTTTA‑TCTCCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑CCCCAGTTGTTGCAG‑ATGG‑CGAAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                            <  2:780428/150‑1 (MQ=255)                                                  ttccttctctcCAACCCATCCCCCGCCTTTA‑TCTCCCCAGGCC‑GGAAGGCCGGGAGAC‑TAAATTGGCGG‑CCCCAGCTGTTGAAG‑ATGG‑CGTAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                             <  2:405525/144‑1 (MQ=255)                                                  ttccATCTCTCAAACCCATTCCCCGCCGTTA‑TCTCCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑CCCCAGTTTTTGCAG‑AGGG‑CGTAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                            <  2:149969/147‑1 (MQ=255)                                                       tctcTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                      >  2:71860/1‑151 (MQ=255)                                                          ctacAACCCATCCACCGCTTTTT‑TCTCCCCGGGCC‑GAAGGGTCAGGAAAC‑GAATTTTGCT‑‑CCCCCTTTTTTTCAGTAGGG‑CTTAC‑CCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                    <  1:1103474/147‑1 (MQ=255)                                                           ttCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                  >  2:651298/1‑151 (MQ=255)                                                           taaactCCATCCCCCGACTTTA‑TCTCCCCAGGCA‑GGGAGGTGAGGAGAC‑TAATTTTGCAG‑CCCCAGTTGTTGCAG‑ATGG‑CCTAC‑CCCCTCTAACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                          <  1:1091499/138‑1 (MQ=255)                                                            tCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                 >  2:902847/1‑151 (MQ=255)                                                              aacccATTCCCCGACTTTA‑TCTCCCCAGGCC‑GGAAGGTCAGGAGAC‑TAATTTTGCTG‑CCCCAGTTGTTGCAG‑ATGG‑CGTAC‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                      <  2:206195/141‑1 (MQ=255)                                                                     tCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                        >  2:627483/1‑151 (MQ=255)                                                                      ccccgccgTTA‑TCTCCCCAGGCC‑GGAGGGTCAGGAGAC‑TGATTTTGCTG‑CCCCATTTTTTGCAG‑ATGG‑CTTAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                        <  1:65156/148‑1 (MQ=255)                                                                      ccccGACTTTT‑TCTTCCCAGGCG‑GGAGGGTCAGGAGAC‑TGAATTTGCTG‑CCCCAGCTGTTGCAG‑ATGG‑CGTAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                         <  2:297819/147‑1 (MQ=255)                                                                      ccccGACTTTT‑TCTCCCCAGGCC‑GGAAGGTCGGGAGAC‑TGAATTTGCTC‑ACCCGGTTTTTGCAG‑A‑GG‑CCTAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                  <  1:109514/138‑1 (MQ=255)                                                                      cACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                        >  2:109866/1‑150 (MQ=255)                                                                       cccgccgTTA‑TCCCCCCAGGCG‑GGAAGG‑CCGGAAACTTGATTTTCCTG‑CCCCGTCTGTTGAAG‑ATGG‑CGAAC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                                      <  1:679050/150‑1 (MQ=255)                                                                       aCCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                                      >  2:280256/1‑151 (MQ=255)                                                                       aCCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                                      >  2:411467/1‑151 (MQ=255)                                                                         cGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                    >  1:1188801/1‑151 (MQ=255)                                                                         cGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                    >  2:627387/1‑151 (MQ=255)                                                                              ttA‑TCTCCCCGGGCC‑GAAGGGTCAGGGGAC‑GAATTTTGCGC‑CCCCAGTTGTTCCAG‑ATGG‑CGTAC‑CCCCCCGACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                   <  2:620942/147‑1 (MQ=255)                                                                                ttttctCCCCGGCCGGGAGGGTCGGGAACC‑TAAATTTGCGC‑CCCCATTTGTGGAAG‑AGGG‑CGACC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                <  1:700592/145‑1 (MQ=255)                                                                                      ccaaGGCC‑GGAAGGCCAGAAGAC‑TGATTTGGCTA‑CCCCAGCTTTTGCGG‑TGGG‑CGTAC‑CCCCCCTGGCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                          <  2:1024865/146‑1 (MQ=255)                                                                                          ggCA‑GGAGGGTCAGGAGGA‑TGAATTTGCTA‑ACCCATGTGTTAGGG‑AGCG‑TTAAG‑CCCCGCCGGCGCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa                                                     <  2:1115009/150‑1 (MQ=255)                                                                                            cc‑GGAAGGGCGGGAAAC‑TGATTTTCCTC‑ACCCGGTTGTTGCAG‑ATGG‑CATAA‑CCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                      <  1:603947/147‑1 (MQ=255)                                                                                                  aGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                              >  1:372525/1‑150 (MQ=255)                                                                                                    gTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                           >  2:357663/1‑151 (MQ=255)                                                                                                               tGATTTTGCTG‑CCCCAGTTGTTGCAG‑ATGG‑CGTAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                    <  2:924025/148‑1 (MQ=255)                                                                                                               tGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCa                                  >  2:1040774/1‑150 (MQ=255)                                                                                                                 aaTTTGCTG‑CCCCAGCTTTTGAAG‑ATGG‑CGTAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                               <  2:258796/151‑1 (MQ=255)                                                                                                                     tGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                           >  2:273511/1‑151 (MQ=255)                                                                                                                      tgcggccccAGTTGTTGAAG‑AGGG‑CGTAC‑CCCCCCGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAacac                             <  2:660170/143‑1 (MQ=255)                                                                                                                        tG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                         >  1:155765/1‑150 (MQ=255)                                                                                                                            cccAGCTTTTGCAG‑AGGG‑CGTAC‑CCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTa                        <  2:947842/148‑1 (MQ=255)                                                                                                                              cagtttttGCAG‑ATGG‑CTAAC‑CCCCCCTACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAccc                    <  1:257696/144‑1 (MQ=255)                                                                                                                                    ttGAAG‑AGGC‑CGTAC‑CCCCCCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAACACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc                 <  1:1056766/147‑1 (MQ=255)                                                                                                                                    tggCAG‑AGGG‑CGTAC‑CCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt              <  2:898421/148‑1 (MQ=255)                                                                                                                                      gCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCACACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGAt           >  2:853589/1‑151 (MQ=255)                                                                                                                                         g‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc          >  2:593607/1‑149 (MQ=255)                                                                                                                                             gg‑CGTAC‑CCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCtt     <  2:209852/151‑1 (MQ=255)                                                                                                                                                 gTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCgg  >  2:164009/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGG  >  NC_000913/3350703‑3350991

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