Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	48.3%	76.6 / 14.0	60	S725A (TCA→GCA)	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 (25/5);  new base C (0/28);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 8.35e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350707‑3351003                                                                                                                                                          |                                                                                                                                                       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:411467/1‑151 (MQ=255)                                                                   aCCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                                                  >  2:280256/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)                                                                                                                                                  ccccACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTcc           <  2:238864/150‑1 (MQ=255)                                                                                                                                                          aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc  >  1:761410/1‑151 (MQ=255)                                                                                                                                                          aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAACGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc  >  2:309568/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTA‑TCTTCCCAGGCC‑GGAAGGTCAGGAGAC‑TGAATTTGCTG‑ACCCAGCTGTTGCAG‑ATGG‑CGTAA‑CCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350707‑3351003

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