Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	175,694	0	T	G	41.7%	95.1 / 14.3	60	G196G (GGT→GGG)	clcA	H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base T (8/27);  new base G (25/0);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 4.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGATTTTTAATCATGAAGTTGCGTTGATTGACGTCGGTAAACTTTCTGAT  >  NC_000913/175549‑175844                                                                                                                                                  |                                                                                                                                                        aaGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTAtttt                                                                                                                                                    <  1:1001865/151‑1 (MQ=255)   gggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTg                                                                                                                                                   <  1:380144/150‑1 (MQ=255)    ggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGtt                                                                                                                                                 <  1:784511/151‑1 (MQ=255)    ggCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGGGCTTGATATTTTTCGCCTGAAAGGGGACGAAGCTCGCCATACGCTGCTGGCCACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATtttttt                                                                                                                                                 >  1:1028319/1‑151 (MQ=255)       cAACCGTGCAGATCGGCGGTAACATTTGTCGTATGGTGGTTGATATTTTTCGCCCGAAAGGGGACGCAAACCGCCATACGCTGCTGGCAACCGGTGCTGCTCCGGGGCTGGCCTCGGCCTTTAACGCGCCCCTGGCGGGGACTTTGttttt                                                                                                                                              >  1:767057/1‑149 (MQ=255)        aaCCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATtttttttttt                                                                                                                                             >  2:356366/1‑151 (MQ=255)           cGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATc                                                                                                                                          <  2:1030076/151‑1 (MQ=255)            gTGCAGATCGGCGGTAACATTTGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCCTACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGGTtttttttttt                                                                                                                                             >  2:821814/1‑147 (MQ=255)            gTGCAGATCGGCGGTAACATTGGCCGTATGGGGCTTGATATTTTCCGCCTGAAAGGGGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATtttttttttt                                                                                                                                             >  2:785187/1‑147 (MQ=255)              gCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGa                                                                                                                                        <  2:913527/150‑1 (MQ=255)               cAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGGGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGGTttttttttttttttaa                                                                                                                                       >  2:705781/1‑144 (MQ=255)                  aTCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAgaga                                                                                                                                   <  2:325462/151‑1 (MQ=255)                    cggcggTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTTTTTTTGAAaaga                                                                                                                                   >  1:444467/1‑149 (MQ=255)                            aCATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGAAGCGTcc                                                                                                                            >  2:337127/1‑148 (MQ=255)                             cATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTTTTTATTATCGGAAAGATGCGTCCg                                                                                                                           >  1:68355/1‑148 (MQ=255)                                     gTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTTTTATCGAAGAGATGCGTCCGCAGGTTCg                                                                                                                   >  2:876356/1‑148 (MQ=255)                                            gCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTTTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑Gtt                                                                                                         >  2:811183/1‑150 (MQ=255)                                               tGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTaa                                                                                                       <  1:38758/149‑1 (MQ=255)                                               tGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAAt                                                                                                      <  1:808808/150‑1 (MQ=255)                                                   atTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCg                                                                                                  <  1:15370/150‑1 (MQ=255)                                                     tttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGGTTTTGTTTATTATCGAAGAGATGCGTCCGCCGTTTTGCTATAA‑GTTAATTTCGGtt                                                                                               >  2:17273/1‑151 (MQ=255)                                                      ttCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGCAGAGATGCGTCCCCAGGTTTGCTTTTC‑GTTAATTTTGAtt                                                                                               >  2:644059/1‑150 (MQ=255)                                                        ccGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTaaa                                                                                            <  2:413861/151‑1 (MQ=255)                                                        ccGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTaa                                                                                             >  2:220762/1‑150 (MQ=255)                                                         cGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTaa                                                                                             >  2:1089283/1‑149 (MQ=255)                                                         cGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGGAGAGATGCGTCCGCAGTTTCGCTTTAC‑GTTAATTTCGGTTaa                                                                                             >  2:1002647/1‑149 (MQ=255)                                                              gAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTTTTTATTATCGAAGAGATGTGTCCCCGGTTTTGACATAC‑GTTATTTTTGCTTaaaaag                                                                                         >  1:419956/1‑145 (MQ=255)                                                              gAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCgg                                                                                        <  2:552809/149‑1 (MQ=255)                                                              gAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCgg                                                                                        <  1:86896/149‑1 (MQ=255)                                                              gAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGGAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTTGATTAAAGCggg                                                                                       >  1:748063/1‑149 (MQ=255)                                                               aaaGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGtat                                                                                     <  1:1101792/151‑1 (MQ=255)                                                               aaaGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAACgg                                                                                        >  1:205285/1‑148 (MQ=255)                                                                aaGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGtatt                                                                                    <  2:244405/151‑1 (MQ=255)                                                                     gACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTTTTTTTTTTCGAAGAGATGCGTCCGCGGTTTTGCTATAC‑GTTAAATTTGATTTAAGCGGgatttatt                                                                                >  2:201541/1‑150 (MQ=255)                                                                      aCGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTTGCTATAC‑GTTTATTTCGGTTAAAGCGGTATTTTTTgg                                                                              >  2:544821/1‑151 (MQ=255)                                                                       cGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGGAGAGATGCGTCCGCCGTTTCGCTCTAC‑GTTAATTTCGATTAAAGCGGgatttatt                                                                                >  1:291651/1‑148 (MQ=255)                                                                         aaGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTgg                                                                              >  1:142563/1‑148 (MQ=255)                                                                             tCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGGTTAAAGCGGGATTTATTGGTGTCAtt                                                                       >  1:1057653/1‑151 (MQ=255)                                                                              cGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTTGATTAAAGCGGGATTTTTTGGTGTCAtt                                                                       >  2:212648/1‑150 (MQ=255)                                                                                ccATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTTTTTTCGAAGAGATGCGTCCGCAGTTTCGCTATACGGTTAATTT‑GATTTAAGCGGTATTTATTTGTGTCAtt                                                                       >  1:636740/1‑148 (MQ=255)                                                                                ccATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGGTTATTATCTGAGAGGTGCGTCCGCAGTTTTGCTATAC‑GTTTATTTCGGTTAAAGCGGGATTTATTGGGGTAGTTAt                                                                     >  2:801042/1‑150 (MQ=255)                                                                                 cATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCCAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGt                                                                   >  2:990008/1‑151 (MQ=255)                                                                                  aTACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTc                                                                  <  2:508976/151‑1 (MQ=255)                                                                                      gctgctGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTTTTATCGAAGAGGTGCGTCCGCAGTTTCGCTATAC‑GGTAATTTCGGTTTAAACGGGATTTATTGGGGTCCTTATTTCGGcc                                                              >  1:706668/1‑151 (MQ=255)                                                                                         gctgGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCAtt                                                           <  1:75800/151‑1 (MQ=255)                                                                                            ggCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATTC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACAATTATg                                                        <  1:348310/151‑1 (MQ=255)                                                                                            ggCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATg                                                        <  1:89464/151‑1 (MQ=255)                                                                                              cAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTGTTGGAGAGGTGAGTCCGCCGTTTTGGCAGAG‑GTTGATTTTGATTAAAGCGGGATTTATTGGGGGAAGTTTGTTGACCCTTATg                                                        >  1:32544/1‑149 (MQ=255)                                                                                              cAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTa                                                      <  1:508255/151‑1 (MQ=255)                                                                                                aCCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTAcc                                                    <  1:640385/151‑1 (MQ=255)                                                                                                   ggTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGa                                                 <  2:472208/151‑1 (MQ=255)                                                                                                      gctgctgcGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTTTTATTGAAGAGATGCGTCCGCAGTTTCGCTATAC‑CTTAATTTCGATTAAAGCGGGATTTATTGGTGTCATTATGTCGACCATTATGTatccgcttt                                              >  2:563126/1‑143 (MQ=255)                                                                                                      gctgctgcGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGTAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGAttt                                              <  1:656359/151‑1 (MQ=255)                                                                                                              ggggCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGATTTTTAATCAt                                      <  2:562526/151‑1 (MQ=255)                                                                                                                    ggctgCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGATTTTTAATCATGAAGt                                 <  2:378301/150‑1 (MQ=255)                                                                                                                          ggCCTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGAAGAGATGCGTCCGCAGGTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGGTTTTTAATCATGAAGTTGCGtt                           >  2:520358/1‑150 (MQ=255)                                                                                                                             cTTTAACGCGCCGCTGGCGGGGATTTTGTTTATTATCGGAGAGATTCGTCCGCAGTTTCGCTATAC‑GTTTATTTCTTTTAAAGAGGGATTTTTTTGGGTCATTTTTTCGACCATTATGTACCGGATTTTTAATCATGAAGTTGCGttgat                        >  1:387750/1‑150 (MQ=255)                                                                                                                                            ggCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGATTTTTAATCATGAAGTTGCGTTGATTGACGTCGGTAAACt         <  2:544681/150‑1 (MQ=255)                                                                                                                                               gggTATTTTGTTTAATGTCGAAGAGATGCGTGCGCAGTCTCGCTATAC‑GTGAATTTCGATTAAAGCGGTATTTATTGGTGGCAGTATGTAGGCCATTAGGTAACGGATATGTAATCATGAAGTTGAGTTGACTGACGTCGGTAGACTTTCt     <  1:383461/151‑1 (MQ=255)                                                                                                                                                 gTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGATTTTTAATCATGAAGTTGCGTTGATTGACGTCGGTAAACTTTCTGa   <  1:127243/151‑1 (MQ=255)                                                                                                                                                  tATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGATTTTTAATCATGAAGTTGCGTTGATTGACGTCGGTAAACTTTCTGAt  >  1:702454/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                        AAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGAAGCTCGCCATACGCTGCTGGCAACCGGTGCTGCTGCGGGGCTGGCTGCGGCCTTTAACGCGCCGCTGGCGGGTATTTTGTTTATTATCGAAGAGATGCGTCCGCAGTTTCGCTATAC‑GTTAATTTCGATTAAAGCGGTATTTATTGGTGTCATTATGTCGACCATTATGTACCGGATTTTTAATCATGAAGTTGCGTTGATTGACGTCGGTAAACTTTCTGAT  >  NC_000913/175549‑175844

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