Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,286,965	0	A	G	25.0%	76.2 / 12.3	61	intergenic (+439/+101)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (25/17);  new base G (6/8);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGC  >  NC_000913/1286823‑1287105                                                                                                                                               |                                                                                                                                             ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                              >  1:930188/1‑143 (MQ=255)    cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                                           >  2:936832/1‑143 (MQ=255)    cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGc                                                                                                                                           >  1:1109072/1‑143 (MQ=255)         ttGGGCAGAATATTTGATTTCGGATTCGTGTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                                      <  1:329025/143‑1 (MQ=255)             gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGCGATGATCAGCGTAGTAGGATTACTCGGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGCGCCgc                                                                                                                                  >  2:604958/1‑142 (MQ=14)              cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTGTCGCTGACGCTCGAGCCGaa                                                                                                                                 >  1:764096/1‑143 (MQ=255)              cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                                 >  1:812942/1‑143 (MQ=255)              cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTATCTCGCTGCCGCTCGCGCCGaa                                                                                                                                 >  1:371379/1‑143 (MQ=255)                gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCgagccgagcc                                                                                                                               >  1:124545/1‑143 (MQ=255)                gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                               >  1:983750/1‑143 (MQ=255)                      ttGATTGCGGATGCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                         <  1:808759/143‑1 (MQ=255)                        gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGa                                                                                                                       <  2:812942/143‑1 (MQ=9)                             cGGATTCGTTTGAGAGTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                       <  2:68759/122‑1 (MQ=9)                             cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                       >  1:68759/1‑122 (MQ=255)                             cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                       >  1:391955/1‑122 (MQ=255)                             cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                       <  2:391955/122‑1 (MQ=255)                                     tttAAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGCAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGACCCTTGGTCGACGCTCCCCATCCt                                                                                                          >  2:1189651/1‑143 (MQ=11)                                        gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGt                                                                                                                          >  2:588051/1‑124 (MQ=255)                                        gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGt                                                                                                                          <  1:588051/124‑1 (MQ=255)                                        gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTGcc                                                                                                       >  1:207258/1‑143 (MQ=255)                                          gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTGGTCAGAGCTGCTGATCCTTcccc                                                                                                     >  2:521260/1‑143 (MQ=11)                                                     cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGa                                                                                          >  1:81933/1‑143 (MQ=17)                                                      ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAgca                                                                                         >  1:914912/1‑141 (MQ=38)                                                       gttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                        <  1:1150312/142‑1 (MQ=38)                                                        ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGTCAGAata                                                                                       >  1:1120392/1‑143 (MQ=255)                                                           aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGTGCCGCACCTTAGTCGAAGCTTCTCATACTTCCCCGCCTGGGCCGAATAttt                                                                                    >  1:13912/1‑143 (MQ=11)                                                           aaaGTGATGGTGGTGGGGGAAGGATTACGCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTCCCGTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                    <  1:311746/143‑1 (MQ=255)                                                              gcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                                 <  1:104491/93‑1 (MQ=255)                                                              gcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                                 >  2:104491/3‑95 (MQ=37)                                                               tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGCACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                >  1:400005/1‑143 (MQ=255)                                                                     tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                          <  1:49705/143‑1 (MQ=255)                                                                       gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGCATATTTGATTGCGGATTc                                                                        >  1:489574/1‑143 (MQ=255)                                                                          ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGCAGCTTCTCATCCTTCCCCGCATGGGCAGAATACTGAATTGCGGCCTCGCt                                                                     >  2:89328/1‑143 (MQ=11)                                                                           gggAGGGGTTCCTCAGCGCTGCGCGCGTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                    <  2:11641/143‑1 (MQ=11)                                                                               aGGATTACTCAGGGCTGCTCGCTTCGCCCTCCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                <  2:723263/143‑1 (MQ=255)                                                                                 gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTCATTGCGGGTGCGCTTGGGAGt                                                              >  1:1053415/1‑143 (MQ=14)                                                                                      cTCCGCGCTTCCCCCTCCGCCCGCCCGGCCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                         <  2:129404/143‑1 (MQ=17)                                                                                       tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGCGCCGAACCTTGGTCGCAGCTTCTCATCCTTCCTCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                        >  1:227423/1‑143 (MQ=255)                                                                                       tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGCGGATTAGCAGCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                        >  1:993878/1‑143 (MQ=17)                                                                                       tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                        >  1:1013016/1‑143 (MQ=255)                                                                                         aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                    <  1:499283/113‑1 (MQ=40)                                                                                         aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                    >  2:499283/1‑113 (MQ=39)                                                                                             cTGCGCGCTTCGCCCTTCGGGGCGTCGGCTGCGGCCGCGCTCTCCCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                  <  2:1109072/143‑1 (MQ=17)                                                                                                cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGCTTGGGAGCTCAGGGACCTTTGaa                                               >  2:709977/1‑143 (MQ=255)                                                                                                 gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGaaa                                              >  1:854924/1‑143 (MQ=255)                                                                                                  cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGCCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAAATTAGACGGCGGAGTCGTCGGGGAGTTCAGGGCCCTTTGagac                                             >  2:1144441/1‑140 (MQ=11)                                                                                                    cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCGCGCATTGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                                           <  1:437963/143‑1 (MQ=255)                                                                                                     tgcGCCGTGCGGGTCGGTTCGTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCGTCGTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                          <  2:662044/141‑1 (MQ=14)                                                                                                      tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGCAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGGCAGTTCGGGGCCTTTTGAAAGGGac                                         >  2:765305/1‑142 (MQ=11)                                                                                                             ccGGGTCGCTGCCTGCGGCACCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                  <  1:344320/142‑1 (MQ=255)                                                                                                              cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTGCTCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                 <  2:1212901/143‑1 (MQ=255)                                                                                                                ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAACGTGATGGTGGTgggg                               >  1:324223/1‑143 (MQ=255)                                                                                                                   cGTTGCCTGCGCCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                            <  1:667936/143‑1 (MQ=255)                                                                                                                   cGTTGCCGGCGGCAATGCTCTCGCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                            <  2:1144615/143‑1 (MQ=255)                                                                                                                     ctGCCTGCGGCAGCGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                          <  2:1137715/142‑1 (MQ=255)                                                                                                                     cggCCGGCGTCAACGCTCTCTCGCTGACGCTCCATCCGAGCCTTAGTCGAATCTTCTCATCCTGCCCCGTATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                          <  2:653244/141‑1 (MQ=17)                                                                                                                      tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTTGTGGGGGAAGGa                         <  1:1362335/143‑1 (MQ=255)                                                                                                                       gCCTGCTGCATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGGATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                        <  1:211913/143‑1 (MQ=255)                                                                                                                       gCCTGCGGCCACGCCCTCCCCCTCGCGCGTGCGCTCAACCTGAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                        <  2:82928/143‑1 (MQ=14)                                                                                                                       gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTCGGTCGAAGCCTCCCATCCTTCCCCGCCTGGGCACAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGACTTTCTAAAGTGCCGCTGGTGGGGGAAGGAt                        >  2:548813/1‑143 (MQ=14)                                                                                                                         cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGTGACGTTTGAAAGTGATGGTGGTGGGGGGAGGACTa                      >  1:443734/1‑143 (MQ=255)                                                                                                                           gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACt                    >  1:343727/1‑143 (MQ=255)                                                                                                                                             tGACCCTCGAGCCGAACCTTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc  <  1:820040/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             TTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGC  >  NC_000913/1286823‑1287105

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