Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,328	0	C	A	42.4%	89.5 / 31.2	59	intergenic (+337/‑76)	yahL/yahM	uncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (30/4);  new base A (0/25);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 7.84e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTG  >  NC_000913/345179‑345467                                                                                                                                                       |                                                                                                                                            aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATc                                                                                                                                             >  2:903191/1‑150 (MQ=255) aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATc                                                                                                                                             >  1:457231/1‑150 (MQ=255)  aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTACATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCaa                                                                                                                                           <  2:700882/151‑1 (MQ=255)       tGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAAcc                                                                                                                                      >  1:133418/1‑151 (MQ=255)          tctcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGcc                                                                                                                                   >  1:1076492/1‑151 (MQ=255)                gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCAtt                                                                                                                                >  1:758146/1‑148 (MQ=255)                          tgtaattgttttgtttAACATTCATATAAAAAAAATCAAAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTCCGCATCGCACTCAAGCCTGAAAAAAAATATGCTGTAAGGCTC‑ATATAAAAAACCGCCATTAGCTCATCAgg                                                                                                                     <  2:283534/143‑1 (MQ=255)                              attgtgTTGGTTGACATTAATATGAAAAAAACCATAATCCCATCATGTTTGTGTGGGTTCTTTTCTTTATCTTCCCAACCCCACTCAGGCCTGAAAAAAAATATGCTGAAAGCCTA‑ATAAAAAAAACCGCCATTAGCTCATCAGGa                                                                                                                    <  1:757062/144‑1 (MQ=255)                              aaTGTGTTGTTTAACATTCATATAAAAAAAATAATAATCCCATCATTTTTGTGTGGGGTTTTTTCTGTATCTTACGCACCCCACTCAACCCTGAAAAAAAATAGGCTGTAAGGCTA‑ATACAAAAAACCGCCATTAGCTCATCAGGa                                                                                                                    <  1:159198/146‑1 (MQ=255)                                  tgtTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                           >  1:849858/1‑151 (MQ=255)                                     tggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                        >  2:61052/1‑151 (MQ=255)                                     tggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                        >  1:428249/1‑151 (MQ=255)                                     tggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                        >  2:703099/1‑151 (MQ=255)                                        ttgACATTCATATAAAAAAAATCAAAATCCCATCATTTTTGTGTGGGTTTTTTTCTGTATCTTACGCACCCCACTAAAGCCTGAAAAAAAATATGCTGAAAGGCTC‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACa                                                                                                       <  1:1001611/149‑1 (MQ=255)                                         tgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                    >  2:17480/1‑151 (MQ=255)                                          gACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                    >  2:637718/1‑150 (MQ=255)                                             aTTAATATAAAAAAAATCATAAGTCCATCATGTTTGGGTGGGGTCTTTTCTGTATCTTCCCCACCGCACTCAAGCCTGACAAAAAATATGCAGTAAGGCTA‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAg                                                                                                  <  1:98325/149‑1 (MQ=255)                                              ttCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt                                                                                               >  1:172353/1‑151 (MQ=255)                                               tCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                              >  1:236394/1‑151 (MQ=255)                                               tCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                              >  2:1027285/1‑151 (MQ=255)                                                 atatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                                            >  1:157082/1‑151 (MQ=255)                                                 atatAAAAAAAATCATAATCCCATCATTTTTGTGTGGGTTCTTTTCTGTATCTTACGCACCGCACTCAAGCCTGACAAAAAATATGCTGAAAGGCTC‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                              <  2:92084/149‑1 (MQ=255)                                                     gAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                                         >  1:166919/1‑150 (MQ=255)                                                      aaaaaaaTAAAAATCCCATCATGTTTGTGTGGGTTCTTTTCTGTATCTTCCGCATCGCACACAAGCCTGACAAAAAATATGCTGAAAGGCTC‑ATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                       <  2:391599/151‑1 (MQ=255)                                                       aaaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                      >  1:1156994/1‑151 (MQ=255)                                                        aaaaaTGATAATCCCAGAATTTTTGGGTGGGTTCTTTTCTGTTTGTACCGATCCCCACTCGGGCGTGAAAGAAATTAGGCT‑TAAGGCCCGATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                     <  1:136671/151‑1 (MQ=255)                                                         aaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                    >  2:1027577/1‑151 (MQ=255)                                                              cATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACggg                                                                               >  1:563764/1‑151 (MQ=255)                                                               atatTTCCACCATGTTTGGGGGGGGTTTTTTTTGTCTCTTACGAATCGCACACAACTTTGACAAAAAATATGCTGTAAGGCTT‑GTAAAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACggg                                                                               <  2:804866/150‑1 (MQ=255)                                                                 aaTTCCATCATTTTGGTGTGGGTTTTTTTCTTTATCTTACCCACCCCCCCAAAGCCTGACAAAAAATTTGTTGTAAGGCTA‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                              <  1:628495/149‑1 (MQ=255)                                                                   ttCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                          >  1:876061/1‑151 (MQ=255)                                                                       atcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc                                                                      >  2:919278/1‑151 (MQ=255)                                                                         catGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                    >  2:816524/1‑151 (MQ=255)                                                                               tgtgtgGGGTTTTTTCTGTATCTTCCGCATCGCCCTCAGGCCTGAAAAAAAATATGCTGTAAGGCTC‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                               <  2:299993/150‑1 (MQ=255)                                                                                   tgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                          >  2:315530/1‑151 (MQ=255)                                                                                    ggttttttttCTGTATTTTCCCCATCCCACTGAAGCCTTACAGAAAAGATGGTGTAAGGCCC‑TGAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                               <  2:1066922/139‑1 (MQ=255)                                                                                    ggggTTTTTTTTTTATCTTACCCATCCCCCCCAAGCCTGAAAAAAAATATGCTGTAAGGCTA‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                             <  1:1109257/147‑1 (MQ=255)                                                                                       gTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                       >  2:506610/1‑150 (MQ=255)                                                                                            ttCTGTTTCTTACGCATCACACCAACGCTTAACAAAAAATATCCGTTAGGGCTA‑ATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                    <  2:822365/148‑1 (MQ=255)                                                                                            ttCTGTATCTTCCCAACCGCACTAAAGCCTGAAAAAAAATATGCTGAAAGCCTA‑ATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                  <  2:356556/150‑1 (MQ=255)                                                                                                gTATCTTCCGCATCGCACTCAAGCCTAAAAAAAAATATGCTGTAAGGCTA‑ATATAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTCCGAGGCCTCGATGGCGGTCCAACTTTTc                                                <  1:223038/148‑1 (MQ=255)                                                                                                gTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaa                                              >  1:618161/1‑150 (MQ=255)                                                                                                     ttACGCATCGCACTAAAGCCTGAAAAAAAATATGTTGTAAGGCTC‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                             <  1:551081/146‑1 (MQ=255)                                                                                                     ttACGCATCGCACCCAAGCCAGAAAAAAAATATGCTGAAGGGCTC‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa                                            <  1:424214/147‑1 (MQ=255)                                                                                                          cATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                    >  1:1134609/1‑150 (MQ=255)                                                                                                             cGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGCTA‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                    <  2:731594/147‑1 (MQ=255)                                                                                                                    aggcataaaaaaaatTATGTTGAAAGGCAC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACATCCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATc                           <  2:524146/138‑1 (MQ=255)                                                                                                                    aaGCCTAAAAAAAAATATGTTATAGGCCTC‑ATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                          <  1:707836/150‑1 (MQ=255)                                                                                                                      gccagaaaaaaaaTATGCTGTAAGGCTC‑ATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                            <  2:532371/137‑1 (MQ=255)                                                                                                                      ccctaaaaaaaaaTATGCTGAAGGCCTA‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCt                       <  1:709086/142‑1 (MQ=255)                                                                                                                       ccTGAAAAAAAATATGTTGTAAGGCTC‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                      <  2:428255/151‑1 (MQ=255)                                                                                                                         tGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                    >  1:979468/1‑151 (MQ=255)                                                                                                                             aaaaaaTATGCTGAAAGGCTA‑ATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                    <  2:1071459/145‑1 (MQ=255)                                                                                                                             aGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                >  1:712466/1‑151 (MQ=255)                                                                                                                                    atGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg         >  1:433800/1‑151 (MQ=255)                                                                                                                                      gttgaagggCTA‑AAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa       <  1:393526/144‑1 (MQ=255)                                                                                                                                        tGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa      >  2:754470/1‑150 (MQ=255)                                                                                                                                        tGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa      >  2:294723/1‑150 (MQ=255)                                                                                                                                           aaGGCTC‑ATATCAAAAACCGCCTTTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg  <  1:113787/151‑1 (MQ=255)                                                                                                                                           aaGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg  >  1:284572/1‑151 (MQ=255)                                                                                                                                             ggCTA‑ATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg  <  2:788484/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                            AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTC‑ATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTG  >  NC_000913/345179‑345467

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