Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	46.9%	103.3 / 10.8	67	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 (28/6);  new base C (0/30);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.33e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAG  >  NC_000913/3350704‑3351000                                                                                                                                                       |                                                                                                                                                    ttttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAc                                                                                                                                                    >  2:298390/1‑151 (MQ=255)  tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGa                                                                                                                                                     <  1:285899/149‑1 (MQ=255)    tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                                 >  2:841641/1‑151 (MQ=255)       cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                              <  1:18201/151‑1 (MQ=255)                tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                                      >  1:137466/1‑150 (MQ=255)                      ctcTTCTACGTCGTGCCCCCACTCTTTTTTCATCTCTTCAACCCATCCCCCGCCTTTTTCTCCCCGGGCGGGAAGGTCAGG‑AGACTGAATTTGCTGACCCTGTTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAAtttt                                                                                                                                <  2:1043238/148‑1 (MQ=255)                         ttCTCCGTCGTGCCCCCCCTCTTCTTTCATCTCTTCAATCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCGGG‑AGACTAATTTTGCTGCCCCATTTTTTGCAGATGGGGAACCCCCACGGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                              <  2:31768/149‑1 (MQ=255)                         ttCTACGTCGTGACGCCCCTCTTCTTTCATCTCTCCAATCCATCCCCCGACGTTATCTCCCCGGGCCGGAAGGTCAGG‑AGACTAATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                              <  1:1024873/149‑1 (MQ=255)                           cTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcc                                                                                                                          >  1:719728/1‑151 (MQ=255)                            tACGTCGTGCCCCCACTCTTCTTTCATCTCTTAAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGG‑AGCCTGATTTTGCGGACCAAGTTTTTGCAGTGGGCGTACCCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                                <  2:1012550/144‑1 (MQ=255)                            acgttcgtacccccACTCTTTTTTAATCTCTTAACCCAATCCCCCAACTTTATCTTCCCCGGCGGGAAGGTCAGG‑AGCCTGAATTTGCTCCCCTGGCTGTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAAtttt                                                                                                                                <  2:1062948/132‑1 (MQ=255)                              cgtcgtGCCCCCCCTCTTCTTTCATCTCTCAAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAGGGTCAGG‑AGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGACCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                        <  1:990965/150‑1 (MQ=255)                              cgccgtGCCCCCCCTCTTCTTTCATCTCTTAAATCGTTACCCCGCCTTTATCTCCCCAGGCCGGAAGTTAAGG‑ATACTGATTTTGCTGCCCAAGCTGTTGCAGAGGGCGAACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGtc                                                                                                                           <  1:687050/144‑1 (MQ=255)                                 cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                    >  2:928750/1‑151 (MQ=255)                                   tgccccccctctTTTTCCATCTCTCAAACCCATCCCCCGCCTTTATCTCCCCAGGCCGAAAGGTCAGG‑AAACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                     <  1:402062/140‑1 (MQ=255)                                    gccccccttcttttTTCCCCCCTCCAATCCATTCCCCGACGTTATCTCCCCAGGCGGGAAGGTCAGG‑AGACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                        <  1:125333/136‑1 (MQ=255)                                    gACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                                 >  2:633692/1‑151 (MQ=255)                                          aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                           >  2:503674/1‑151 (MQ=255)                                              ttcttcCCTCTCTCCAACCCATCCCCCGCCTTTATCTCCCCAGGCCGGAGGGTCAGG‑AGCCTGATTGTGCTGCCCCGGTTTTTGCAGATGGCCAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                         <  2:51930/149‑1 (MQ=255)                                               tcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                      >  2:749447/1‑151 (MQ=255)                                                   tCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                                  >  2:487657/1‑151 (MQ=255)                                                   tCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                                  >  1:678964/1‑151 (MQ=255)                                                       ctctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                              >  2:361359/1‑151 (MQ=255)                                                         ctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc                                                                                            >  2:1034684/1‑151 (MQ=255)                                                             aacccttcccccGACTTTATCTCCCCGGGCGGGAAGGCCAGG‑AGACTAATTTTGCGGCCCCAGTTTTTGCAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGt                                                                                         <  2:313899/140‑1 (MQ=255)                                                             aacccttcccccGACTTTATCTCCCCAGGCGGGAAGGTGGGG‑AGACTAAATTTGCTCACCCGGTTTTTGCAGTGGGCATACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                        <  2:478053/141‑1 (MQ=255)                                                             aacccttcccccGACGTTATCTCCCCAGGCCGAAAGGTCAGG‑AGACTAATTTTGCTCCCCCAGCTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                               <  2:378613/134‑1 (MQ=255)                                                             aacccATCCCCCGACTTTATCTTCCCAGGCGGGAAGGTCAGG‑AGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCATACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                               <  1:298446/141‑1 (MQ=255)                                                                caatgcACCGACGTTATCGTCCCAGGCGGAGAGGGCAGG‑AGACGGAATTTGCGGGCCCAGGTTTTGCAGATGCCGTACCCCCACGGGCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                        <  1:533941/127‑1 (MQ=255)                                                                   ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                                  >  1:882680/1‑151 (MQ=255)                                                                    cccccGTCGTTTTCTTCCAAGCCCGGAAGG‑CAGGTGGACGGACTTTGCTAACCCGCTTGTTGTAGAGGGCGAAACCCCATCGTCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                 <  1:226251/148‑1 (MQ=255)                                                                     ccccGACTTTATCTCCCCAGGCGGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                                  <  1:974047/147‑1 (MQ=255)                                                                     ccccGACTTTATCTCCCCAGGCCGGAGGGTCAGG‑AGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCTAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                                  <  2:350039/147‑1 (MQ=255)                                                                     ccccGACGTTATCCCCCCAGGCCGGAAGGCCAGG‑AGACTAAATTTGCTCCCCCGGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                 <  1:239025/148‑1 (MQ=255)                                                                        cGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                             >  2:1078174/1‑151 (MQ=255)                                                                              tATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                       >  2:573865/1‑151 (MQ=255)                                                                              tATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                       >  2:234881/1‑151 (MQ=255)                                                                                 cTTCCCAGGCCGGACGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt                                                                    >  2:593561/1‑151 (MQ=255)                                                                                    cccAGGCCGGAAGGTCAGG‑AACCTAAATTTCCTCCCCCAGTTTTTGCAGATGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                                 <  2:613700/151‑1 (MQ=255)                                                                                    ccagggCGGGAAGGTCGGG‑AGACTGAATTTGCTGACCCAGTTGTTGAAGAGGGCGAACCCCCACGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGa                                                                  <  1:301891/146‑1 (MQ=255)                                                                                     ccAGGCCGAAGGGCAAGG‑AGCTGGATTTTGCTGCCCCATTTTTTGTAGATGGCTTACCCCCATTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACa                                                                <  2:1068470/151‑1 (MQ=255)                                                                                            ggAAGGTCAGG‑AGACTGAATTTCCTGACCCAGTTTTTGAAGGTGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                          <  1:762737/150‑1 (MQ=255)                                                                                            gaaaGGTCAGG‑AGACGAATTTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCg                                                         <  2:140061/149‑1 (MQ=255)                                                                                              agggTCGGG‑AGACTAATTTTGCTCCCCCAGCTGTTGCAGAGGGCTAACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                        <  1:18150/148‑1 (MQ=255)                                                                                              aaGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                       >  1:239378/1‑151 (MQ=255)                                                                                                 ggttagcaggcggattttGCTGACCCATCTGTTCGAGAGGGCGAACCCCCCCTGTCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcat                                                     <  1:1057251/136‑1 (MQ=255)                                                                                                    aGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCatcatc                                                 >  1:459323/1‑151 (MQ=255)                                                                                                      g‑agaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                               >  1:272758/1‑151 (MQ=255)                                                                                                         gaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                             >  2:555837/1‑151 (MQ=255)                                                                                                          aCTGATTTTGCTGCCCCAGTTGTGGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTtctc                                            <  1:593686/151‑1 (MQ=255)                                                                                                             gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                                          >  1:1046908/1‑150 (MQ=255)                                                                                                              attttCCTGACCCGGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCa                                           <  2:356069/146‑1 (MQ=255)                                                                                                                 ttGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                                     >  2:535626/1‑151 (MQ=255)                                                                                                                   gCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                                    <  2:213538/150‑1 (MQ=255)                                                                                                                   gCTGACCCAGCTTTTGCAGAGGGCGAACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                                     <  2:418942/149‑1 (MQ=255)                                                                                                                   gCTCACCCAGTTGTTGGAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                                    <  1:642230/150‑1 (MQ=255)                                                                                                                    ctcccccGGCTGTTGCAGATGGCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                                   <  2:838347/146‑1 (MQ=255)                                                                                                                     tgccccGGCTTTTGCAGATGGCTTAACCCCTCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                                   <  1:488550/146‑1 (MQ=255)                                                                                                                     tcccccAGTTTTTGCAGATGGCTAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                                     <  2:1118016/144‑1 (MQ=255)                                                                                                                       aCCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc                               >  1:801375/1‑151 (MQ=255)                                                                                                                          cAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCg                            >  2:644270/1‑151 (MQ=255)                                                                                                                           aGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                           <  1:344217/151‑1 (MQ=255)                                                                                                                             cTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                         >  1:256128/1‑151 (MQ=255)                                                                                                                                  gCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                      >  2:386942/1‑149 (MQ=255)                                                                                                                                    aGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg                      <  1:742548/147‑1 (MQ=255)                                                                                                                                     gATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa                  <  1:1010794/150‑1 (MQ=255)                                                                                                                                             aaCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc         >  1:630439/1‑151 (MQ=255)                                                                                                                                                    cggACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAg  <  1:472827/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAG  >  NC_000913/3350704‑3351000

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: ‑