Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	42.0%	137.1 / 17.2	81	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 (37/10);  new base C (0/34);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 3.82e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACA  >  NC_000913/3350705‑3350998                                                                                                                                                       |                                                                                                                                                  tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTTATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGa                                                                                                                                                   <  1:946075/149‑1 (MQ=255) tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAc                                                                                                                                                  <  1:1476881/150‑1 (MQ=255)   tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAccc                                                                                                                                                >  2:274615/1‑150 (MQ=255)    gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                               <  2:1382852/150‑1 (MQ=255)    gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                               <  2:1204468/150‑1 (MQ=255)      cacccacGCTTTCAGCACTTCTCCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                            <  2:1207423/151‑1 (MQ=255)      cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                            >  1:1007500/1‑151 (MQ=255)      cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                            >  1:414897/1‑151 (MQ=255)             gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCtt                                                                                                                                     >  2:436757/1‑151 (MQ=255)               tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                                    >  1:951/1‑150 (MQ=255)               tttCACCACTTCTCCGTGGTGCCCCCCCTCTTTTTTCTTCTCTTCAACCCATCCCCCGACGTTTTCTCCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGCCCTAGTTGTTGCAGTTGACGCACCCCCTCTGACCCCGCCGCACCtt                                                                                                                                     <  2:1496333/149‑1 (MQ=255)               ttccAGCACTCCTCCGTCGTGCCCCCACCCTTCTTCCATCTCTCAAACCCATCCCCCGACGTTACCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTAATTTTGCTGACCCAGCTGTTGAAGATGGCGTACCCCCCCGGACCCCGCCGCACCt                                                                                                                                      <  2:797117/145‑1 (MQ=255)                ttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTa                                                                                                                                   >  2:348041/1‑150 (MQ=255)                   aGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAttt                                                                                                                               >  2:1455354/1‑151 (MQ=255)                     cACTTCTCCGTCGTGCCCCCACTCTTCTTTCATCTCTTCAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGG‑TCA‑GGAGACTGATTTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAAtttt                                                                                                                              <  1:21049/150‑1 (MQ=255)                             cgtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                                     >  2:63666/1‑151 (MQ=255)                              gtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                                    >  1:775912/1‑151 (MQ=255)                               tcgtGACCCCCCTCTTCTTTCATCTTTTCAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGG‑TCA‑GGAGACTAAATTTGCTGCCCCAGCTGTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                                     <  1:100418/149‑1 (MQ=255)                                cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                  >  1:323900/1‑151 (MQ=255)                                 gtGACGCCCCTCTTCTTTCATCTCTTCAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAGGG‑TCA‑GGAGACTAAATTTGCTGCCCCAGTTTTTGCAGAGGGCGTAACCCCCCTGCCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                      <  1:1344182/146‑1 (MQ=255)                                  tgccGCCCCTCTTCTTCCTTCTCTTAACCCCTCCCCCCGACTTTTTCTCCCCAGG‑CGGAAGGAGCG‑GGATCCTGATTTGGCTGCCCCTGTTTTTGCAGAGGCCGTACCCCCATCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                                 <  2:844211/147‑1 (MQ=255)                                       cctctctTCTTCCACCTCTCAAACCCAGCCACCGACTTTATCTTCCCAGGCCGGAGGG‑TCA‑GGAAATTTAATTTGCTGCCCTAGTTGTTGCTGAGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                                    <  2:628371/139‑1 (MQ=255)                                       cccctttttttccctctcTTAAACCCATCCCCCGACTTTATCTCCCCAGGCGGGAGGG‑TCA‑GGAGACTAATTTTCCTCCCCCAGTTTTTGCAGAGGGCTAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                                            <  2:1243614/136‑1 (MQ=255)                                       cccccctttttCCATCTCTCCAACCCACCCCCCGACTTTTTCTCCCAGGGCAGGAAGG‑CAA‑GGAGACTGATTTTCCTGCCCCACTTGTTGCGGTTGCCCTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                              <  2:626408/143‑1 (MQ=255)                                         aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                         >  1:1187552/1‑151 (MQ=255)                                             ttctttctTCTCTTCAACCCATCCCCCGTTTTTATCCCCCGAGGCGGGAAGG‑TCA‑GGAGTCTAATTTTCCTGACCCAGTTTTTGAGGAGGGCATACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                       <  2:843146/149‑1 (MQ=255)                                             ttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                                     >  1:1403438/1‑151 (MQ=255)                                              tcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                    >  1:447138/1‑151 (MQ=255)                                               cttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                   >  2:867165/1‑151 (MQ=255)                                                ttccttctcTACAACCCATCCCCCGCTGTTATCTTCCCAGGCGGGAAGG‑TCA‑GGAGACTGATTTTGCTGCCCCATCTTTTGCATATGGCTTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                   <  2:347640/145‑1 (MQ=255)                                                 ttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                                 >  1:308511/1‑151 (MQ=255)                                                     tctcTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                             >  1:796090/1‑151 (MQ=255)                                                         ttcaccccATTCCCCGCCGTTATCTCCCCAGCCCGAAGGG‑CCG‑GGAGTCTGATTTGGCGGCCCCAGTTGTTGAAGATGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc                                                                                          <  2:617123/144‑1 (MQ=255)                                                         ttCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                         >  1:337366/1‑151 (MQ=255)                                                         ttCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                         >  1:1378619/1‑151 (MQ=255)                                                             aTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                      >  2:137458/1‑150 (MQ=255)                                                              tCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCATCAATGCCTTTTTTATCCTGCGCCGTc                                                                                      >  2:996655/1‑149 (MQ=255)                                                               caaTTCCCCGCCTTTATCTCCCCAGGCCGGAGGA‑TCA‑GGAGACTGATTTTGCTGCCCCAGTTTTTGCGGATCGCAAAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                         <  2:743287/143‑1 (MQ=255)                                                                cATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                      >  1:1222477/1‑147 (MQ=255)                                                                   tCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                               >  2:546059/1‑151 (MQ=255)                                                                   tCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                               >  1:1487264/1‑151 (MQ=255)                                                                      ccgcctttATCTTCCCAGGCGGGAGGG‑TCG‑GGAGCCTAAATTTCCTGCCCTAGTTTTTGAAGTTGGCTTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                                             <  2:518991/144‑1 (MQ=255)                                                                      ccgcctttATCTCCCCAGGCGGGAGGG‑TCG‑GGAGACTGATTTTGCTGCCCCGGTTGTGGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                                  <  1:848501/139‑1 (MQ=255)                                                                        gACGTTGTCCTCCGGGGCGGGAAGG‑TGG‑GGAGCCTGAATTTTCTGGCCCAGCTGTTGCAGACGGCATAACCCCACAGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                          <  1:1388023/151‑1 (MQ=255)                                                                          cGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                          >  2:1056288/1‑149 (MQ=255)                                                                            tttatcccccAGGCCGAAAGG‑TGA‑GGAGACGGATTTTCCTGCCCCAGTTTTTACAGAGGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                         <  2:565030/125‑1 (MQ=255)                                                                               tCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGc                                                                   >  2:1471916/1‑151 (MQ=255)                                                                                 ttCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                                 >  1:159900/1‑151 (MQ=255)                                                                                   cccGGGCCGGAAGG‑TCG‑GGAGATTGAATTTGCTGCCCCAGTTGTTGCAGTTGGCGTAACCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt                                                                  <  2:1287093/148‑1 (MQ=255)                                                                                   cccAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                               >  2:674110/1‑151 (MQ=255)                                                                                        gCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCTGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGc                                                          >  2:1523513/1‑151 (MQ=255)                                                                                           ggAGGG‑TCG‑GGGGACTGAATTTGCTGCCCCGGGTGTTGAAGAGGGCGAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAg                                                           <  1:725757/147‑1 (MQ=255)                                                                                            gAAGG‑GCAGGGAAACTAATTTTGCTGCCCCGGCTTTTGAAGAGGGCTTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                        <  2:183428/150‑1 (MQ=255)                                                                                             aaGT‑TCA‑GGAGACTAATTTTGCTGCCCCAGCTTTTGCAGAGGGCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                     <  2:97708/151‑1 (MQ=255)                                                                                             aaGG‑TCA‑GGAGACTGATTTTGCTGCCCCAGCTGTTGAAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                      <  1:367942/150‑1 (MQ=255)                                                                                              aGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                     >  1:492694/1‑150 (MQ=255)                                                                                               gg‑TCA‑GGAGACTAATTTTCCTGCCCCATCTGTTGCAGAGGGCGAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                     <  2:580403/149‑1 (MQ=255)                                                                                                  tCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                                  >  1:1026776/1‑150 (MQ=255)                                                                                                   cA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                                                >  2:138889/1‑151 (MQ=255)                                                                                                      ggAGATTAATTTTGCTGCCCCAGTTGTTCAAGAGGGCTAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatca                                                 <  1:348055/148‑1 (MQ=255)                                                                                                      ggAGACTTATTTTGCTGACCCAGCTGTTGCAAAGGGCAAACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatca                                                 <  2:958881/148‑1 (MQ=255)                                                                                                              aaTTTGCTGCCCCCGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                                      <  1:753650/151‑1 (MQ=255)                                                                                                              aaTTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                                      >  2:1479510/1‑151 (MQ=255)                                                                                                                    cTGACCCAGCTGTTGGAGATGGCGTAACCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                                  <  1:122963/149‑1 (MQ=255)                                                                                                                     tGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                                >  2:744566/1‑150 (MQ=255)                                                                                                                       aCCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc                             >  2:351672/1‑151 (MQ=255)                                                                                                                              tGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                         <  1:213682/148‑1 (MQ=255)                                                                                                                                 tacAGAGGGCGAAACCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt                     <  2:298653/147‑1 (MQ=255)                                                                                                                                   cAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGAt                  <  1:196520/150‑1 (MQ=255)                                                                                                                                    agagGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa                <  1:975897/151‑1 (MQ=255)                                                                                                                                     gttGGCGTCCCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAg               <  2:997451/149‑1 (MQ=255)                                                                                                                                      aTGGCGTACCCCCCCTACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc                 <  1:878756/148‑1 (MQ=255)                                                                                                                                      aTGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGc              >  1:262173/1‑151 (MQ=255)                                                                                                                                       gggCGTACCCCCTCGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt             <  1:401672/150‑1 (MQ=255)                                                                                                                                        gccgaacccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCtt            <  2:254722/144‑1 (MQ=255)                                                                                                                                        gccGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAg               <  1:421781/146‑1 (MQ=255)                                                                                                                                         gcgaacccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc           <  2:498927/145‑1 (MQ=255)                                                                                                                                         gCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc           >  1:134533/1‑151 (MQ=255)                                                                                                                                          cttacccccACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTc           <  1:1144090/145‑1 (MQ=255)                                                                                                                                          cGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg          <  1:1458641/151‑1 (MQ=255)                                                                                                                                          cGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg          >  2:811654/1‑151 (MQ=255)                                                                                                                                                  cACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTaaa  <  1:276200/151‑3 (MQ=255)                                                                                                                                                       |                                                                                                                                                  TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGG‑TCA‑GGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACA  >  NC_000913/3350705‑3350998

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