Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,328	0	C	A	43.5%	80.1 / 37.8	62	intergenic (+337/‑76)	yahL/yahM	uncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (23/12);  new base A (0/27);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 1.70e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGA  >  NC_000913/345179‑345469                                                                                                                                                        |                                                                                                                                              aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCa                                                                                                                                              <  2:161771/151‑1 (MQ=255) aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCa                                                                                                                                              <  1:42321/151‑1 (MQ=255)  aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCaa                                                                                                                                             <  1:82614/151‑1 (MQ=255)   aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCaaa                                                                                                                                            <  1:76865/151‑1 (MQ=255)      aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAAc                                                                                                                                         <  2:502312/151‑1 (MQ=255)        gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCg                                                                                                                                       >  1:437033/1‑151 (MQ=255)        gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGCGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAAcc                                                                                                                                        <  1:217754/150‑1 (MQ=255)             cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGc                                                                                                                                      <  1:691411/147‑1 (MQ=255)             cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCAtt                                                                                                                                  >  1:205501/1‑151 (MQ=255)                gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGc                                                                                                                               >  2:172144/1‑151 (MQ=255)                gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGc                                                                                                                               >  1:22536/1‑151 (MQ=255)                          tgtgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa                                                                                                                      >  1:440432/1‑150 (MQ=255)                           gtgAATGTGTTGGTTGACATTCATAAAAAAAAAATCAAAATCCCATCATGTTTGTGGGGGTTTTTTTCTTATTCTTACCCACCGCACCCAAGCCTAACAAAAAA‑TATGC‑TGTAGGTCTCATATAAAAAACCGCCATTAGCTCATCAGGaa                                                                                                                     <  2:159521/150‑1 (MQ=255)                              aaTGTGTTGGTTGCCTTTCATAAAAAAAAAATCATAATCCCATCATTTTTGTGTGGGGTTTTTTCTGTATCTTACGCTCCGCACTCAACCCTGAAAAAAAA‑TATGC‑TGTAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAag                                                                                                                    <  1:248725/148‑1 (MQ=255)                                    ttggttgACATTCATATAAAAAAAACCATAATCCCATCTTGTTTGTGGGGGTTTTTTTCTTTATCTTACGCACCGCACTCAGGCCTGACAAAAAA‑TATGC‑TGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacg                                                                                                            <  1:296885/150‑1 (MQ=255)                                      gtttgACATTCATAAAAAAAAAATCATAATCCCATCATGTTTGTGTGGGGTTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAA‑TATGC‑TGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                             <  1:11421/145‑1 (MQ=255)                                        ttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                       >  1:139643/1‑151 (MQ=255)                                        ttgACATTAATATAAAAAAAACCTAAATCCCATCATTTTTGTGGGGGTTTTTTTCTTTATCTTACCCACCCCACACAAGCCTAAAAAAAAA‑TATGC‑TGAAAGCCTCAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACaa                                                                                                        <  2:34230/150‑1 (MQ=255)                                             aTTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                   >  1:458046/1‑150 (MQ=255)                                                 atatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                                              >  2:506429/1‑151 (MQ=255)                                                   aaaaaaaaaaTCAAAATCCCATCTTGTTTGGGGGGGGTTTTTTCTGTATTTTCCGCATCGCCCTCAAGCCGGAAAAAAAA‑TATGC‑TGAAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttg                                                                                               <  1:271546/145‑1 (MQ=255)                                                    tGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                                            >  1:244135/1‑150 (MQ=255)                                                      aaaaaaaaCAAAATCCCATATTGTTTGGGGGGGTTTTTTTCTTTCTCTACCCAACCGCCCCCAACCCAAAAAAAAAA‑TATGC‑GGTAGGGCTCATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                                            <  2:359487/148‑1 (MQ=255)                                                      aaaaaaaTCATAATCCCATCATTTTTGTGTGGGTTTTTTTCTGTATCTTCCGCACCGCACTAAAGCCGGAAAAAAAA‑TATGC‑TGTAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                         <  1:458052/151‑1 (MQ=255)                                                      aaaaaaaTCAAAATTCCATCTTTTTTGTGTGTGGTTTTTTCTTTATCTTACGCACCCCACACAACCCTAACAAAAAA‑TATGC‑TGTAAGACTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                                           <  2:491525/149‑1 (MQ=255)                                                          aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                      >  1:284804/1‑150 (MQ=255)                                                              cATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACggg                                                                                 >  1:164069/1‑151 (MQ=255)                                                               aaaaTCCCATAATTTTTGGGGGGGTTTTTTTTTTTTTCTCCCCCACCGCACTCAAGCCTAAAAAAAAA‑TATGG‑TGAAAGCCTAATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg                                                                                   <  2:186138/146‑1 (MQ=255)                                                                 atgcctatcatGTTTGTGGGGGTTTTTTTCTTTTTCTTCCCCCTCGCACTAAACCCTGA‑AAAAAA‑TATTCTTGAAAGACCGACTTCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                                <  1:528433/143‑1 (MQ=255)                                                                   ttCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                            >  2:289285/1‑151 (MQ=255)                                                                     ccaccttttttGGGTGGGGTCTTTTCGGTATATTACACACCCCCCCCAAGCTGGACAGAAAA‑AATGC‑TGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGa                                                                           <  1:335336/142‑1 (MQ=255)                                                                      catcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                                         >  2:230326/1‑151 (MQ=255)                                                                           tttttGTGGGGGTTTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAA‑TATGC‑TGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                      <  1:34308/147‑1 (MQ=255)                                                                           tGTTGGTGGGGGGTCTTTTCTGTTTCTTCCGCACCGCACACAACCCGAAAAAAAAA‑GATGC‑TGAAGGGCAAATATCAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGcc                                                                       <  2:356646/148‑1 (MQ=255)                                                                             tttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGa                                                                   >  1:588725/1‑150 (MQ=255)                                                                             tttGTGGGGGTTTTTTTCTGTATCTTACCCATCGCACTAAACCCTGACAAAAAA‑TATGT‑TGTAAGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                     <  1:432498/148‑1 (MQ=255)                                                                               tgtgtgGGGTTTTTTCTGTTTCTTACACACCGCACACAACCCTGACAAAAAA‑TATGC‑TGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                                 <  2:237222/150‑1 (MQ=255)                                                                                    gggtttttttCTGTTTCTAACCCCCCGCCCTAAACCGTGAAAAAAAA‑TATGC‑TGAAGCGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                           <  2:54305/145‑1 (MQ=255)                                                                                    ggggTTTTTTTTTTATTTTACCCATCCCACTAAAGCCTGAAAAAAAA‑TATGT‑TGTAGGGCAAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                                 <  1:541779/145‑1 (MQ=255)                                                                                     gttttttttttttttCTTCCGCATCTCACCAAAGCCTAAAAAAAAA‑AATGT‑TATAGTTCAAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcc                                                          <  2:681185/137‑1 (MQ=255)                                                                                        tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                       >  2:279676/1‑151 (MQ=255)                                                                                          ttttCTATATCTCACGCACCCCAC‑CAAACCAGAAAAAAAATTATGA‑GGTAGGCCTCATTTAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                <  2:407204/140‑1 (MQ=255)                                                                                            ttCTTTATCTTCCGCATCGCACTCAAGCCTAACAAAAAA‑TATGC‑TGTAAGGCTAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                    <  1:515870/150‑1 (MQ=255)                                                                                             tCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTc                                                  >  1:263479/1‑151 (MQ=255)                                                                                                   tCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                                            >  2:64820/1‑151 (MQ=255)                                                                                                             cGCACTCAAGCCAGAAAAAAAA‑TATGC‑TGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGaa                                    <  1:556776/149‑1 (MQ=255)                                                                                                             cGCACTAAAGCCAGACAAAAAA‑TATGC‑TGTAAGTCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt                                   <  1:603422/150‑1 (MQ=255)                                                                                                                aCTCAACCTTGAAAAAAAA‑TATGC‑TGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                 <  1:224850/149‑1 (MQ=255)                                                                                                                 ctaaaccctaaaaaaaaa‑TATGC‑TGAAGGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa                               <  2:288840/136‑1 (MQ=255)                                                                                                                 cccAAGCCTGACAAAAAA‑TAGGC‑TGTAAGCCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa                               <  1:718359/148‑1 (MQ=255)                                                                                                                  tCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                  >  1:360867/1‑146 (MQ=255)                                                                                                                    aaGCCTGAAAAAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                              <  1:29423/148‑1 (MQ=255)                                                                                                                     accgtaaCAAAAAA‑TATGT‑TTTAGGGTTAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                              <  2:690350/141‑1 (MQ=255)                                                                                                                            cAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                   >  2:29877/1‑151 (MQ=255)                                                                                                                            cAGAAAA‑TATGC‑TGTAAGCCTCATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatt                     <  2:602859/149‑1 (MQ=255)                                                                                                                               aaaa‑TATGT‑TGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcat                <  1:132066/151‑1 (MQ=255)                                                                                                                               aaaa‑TAGGT‑TGAAAGGCCCATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatca                 <  1:277790/150‑1 (MQ=255)                                                                                                                                  a‑tatGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATccc             >  2:635370/1‑151 (MQ=255)                                                                                                                                             agggcctatatAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg    <  2:716883/144‑1 (MQ=255)                                                                                                                                             aaGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg    >  1:478648/1‑151 (MQ=255)                                                                                                                                                gCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGa  >  1:89309/1‑150 (MQ=255)                                                                                                                                                  tCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtta                       >  1:94242/1‑127 (MQ=255)                                                                                                                                                        |                                                                                                                                              AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAA‑TATGC‑TGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGA  >  NC_000913/345179‑345469

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