Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,577,123	0	A	C	24.6%	124.9 / 13.1	61	V166G (GTA→GGA)	yddB	putative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (34/12);  new base C (0/15);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 2.46e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGC  >  NC_000913/1576986‑1577262                                                                                                                                          |                                                                                                                                            ttCCTGAAGGGCCTTGTTTAAATGCCCTCTTGTTATTCTCATCGATTTCGGATGTTAACCGGTCCGACCGGGAAGTGCGTTAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCATTAAcc                                                                                                                                           >  2:220297/1‑140 (MQ=255) ttCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcc                                                                                                                                           >  1:283569/1‑140 (MQ=255)  tCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATACCCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTcccc                                                                                                                                          <  1:360628/140‑1 (MQ=255)  tCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAccc                                                                                                                                          >  2:283633/1‑140 (MQ=255)   ccTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcccc                                                                                                                                         >  2:950037/1‑140 (MQ=255)   ccTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcccc                                                                                                                                         >  2:37311/1‑140 (MQ=255)      gAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGcc                                                                                                                                      >  2:720319/1‑140 (MQ=255)      gAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGcc                                                                                                                                      >  2:69522/1‑140 (MQ=255)        agagCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCAt                                                                                                                                    >  2:672683/1‑140 (MQ=255)            ccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAccc                                                                                                                                          <  1:674656/130‑1 (MQ=255)            ccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAccc                                                                                                                                          >  2:154633/1‑130 (MQ=255)            ccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAccc                                                                                                                                          <  1:154633/130‑1 (MQ=255)            ccTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAccc                                                                                                                                          >  2:674656/1‑130 (MQ=255)             cTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAg                                                                                                                               >  2:376449/1‑140 (MQ=255)              ttGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGc                                                                                                                              >  1:322530/1‑140 (MQ=255)                gATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGa                                                                                                                            <  1:69522/140‑1 (MQ=255)                      aTGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt                                                                                                                      >  2:86008/1‑140 (MQ=255)                           ctctTGTTATTCCCATCGATATGCGATGTTAACCAGCCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAc                                                                                                                 <  2:730610/140‑1 (MQ=255)                                tttATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGca                                                                                                            <  2:283569/139‑1 (MQ=255)                                     tctcATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa                                                                                                       <  2:58208/140‑1 (MQ=255)                                       tcATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACt                                                                                                     <  1:376449/140‑1 (MQ=255)                                        cATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTg                                                                                                    >  1:653180/1‑140 (MQ=255)                                        cATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTg                                                                                                    >  1:223578/1‑140 (MQ=255)                                         aTCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGt                                                                                                   >  1:638455/1‑140 (MQ=255)                                           cGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCa                                                                                                 >  2:903981/1‑140 (MQ=255)                                             atatGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATa                                                                                               <  2:322530/140‑1 (MQ=255)                                             atatGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATa                                                                                               <  2:653180/140‑1 (MQ=255)                                              tatGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa                                                                                              >  1:890707/1‑140 (MQ=255)                                                 gCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGc                                                                                           >  1:299783/1‑140 (MQ=255)                                                 gCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGc                                                                                           >  1:306321/1‑140 (MQ=255)                                                  cgttGTTAACCAGTCCGACCGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCg                                                                                          <  2:638455/137‑1 (MQ=255)                                                  cGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCg                                                                                          <  1:672683/140‑1 (MQ=255)                                                     tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTca                                                                                       >  2:794134/1‑140 (MQ=255)                                                       ttAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcaca                                                                                     <  2:299783/140‑1 (MQ=255)                                                               tCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg                                                                             >  2:597621/1‑140 (MQ=255)                                                               tCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg                                                                             >  2:295695/1‑140 (MQ=255)                                                                 cGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa                                                                           >  2:688655/1‑140 (MQ=255)                                                                   aaCGCGTAGGGCGATAACCCAATTTCACTTTACTAAAAGCAGGGTTTAGGCTTTTAATCTTTCCATAAATTCCCCCGCCATTAAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAg                                                                         <  2:938494/140‑1 (MQ=255)                                                                     cgcgTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt                                                                       <  2:722874/140‑1 (MQ=255)                                                                     cgcgTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt                                                                       <  1:220297/140‑1 (MQ=255)                                                                     cgcgTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCt                                                                       <  1:543315/140‑1 (MQ=255)                                                                        gTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGAc                                                                    <  1:37311/140‑1 (MQ=255)                                                                          aGTGCGATACCCAAATTTCCCCTTGCTACCATCAGCGTTGAGGCGTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACAt                                                                  <  2:223578/140‑1 (MQ=255)                                                                             gCGATACCCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa                                                               <  2:890707/140‑1 (MQ=255)                                                                                aTAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGat                                                            >  2:377444/1‑140 (MQ=255)                                                                                  aaCCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                          <  2:752873/140‑1 (MQ=255)                                                                                  aaCCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                          <  1:281709/140‑1 (MQ=255)                                                                                  aaCCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa                                                          <  1:597621/140‑1 (MQ=255)                                                                                   ccccAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATAAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat                                                         <  2:144646/139‑1 (MQ=255)                                                                                    cccAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataata                                                        >  2:330456/1‑140 (MQ=255)                                                                                     ccAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAc                                                       <  1:794134/140‑1 (MQ=255)                                                                                               cTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAATACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACa                                             >  1:772967/1‑140 (MQ=255)                                                                                               cTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACa                                             <  1:283633/140‑1 (MQ=255)                                                                                                ttGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACAt                                            <  1:377444/140‑1 (MQ=255)                                                                                                 tGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATc                                           >  1:836453/1‑140 (MQ=255)                                                                                                   cTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATccc                                         >  2:633890/1‑140 (MQ=255)                                                                                                             cGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTg                               >  2:35155/1‑140 (MQ=255)                                                                                                                  aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGc                          >  2:852093/1‑140 (MQ=255)                                                                                                                   aGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCa                         <  1:503358/140‑1 (MQ=255)                                                                                                                      ttttGATCTTTGCATCATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAATACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCacta                      <  2:772967/139‑1 (MQ=255)                                                                                                                         tGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTg                   >  1:861809/1‑140 (MQ=255)                                                                                                                             ttttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCAt               <  1:950037/139‑1 (MQ=255)                                                                                                                                       aTTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCAtt     >  2:794126/1‑140 (MQ=255)                                                                                                                                       aTTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCAtt     >  2:806737/1‑140 (MQ=255)                                                                                                                                        ttACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTc    >  1:882189/1‑140 (MQ=255)                                                                                                                                          aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGc  >  2:802242/1‑140 (MQ=255)                                                                                                                                          aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGc  >  1:257121/1‑140 (MQ=255)                                                                                                                                          aCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGc  >  1:654279/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                            TTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCCCTGTGACATCCCGCTAATATTGGTTGCACTACTGGCATCGGACTCATTCGC  >  NC_000913/1576986‑1577262

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