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

AAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTG  >  NC_000913/345182‑345475                                                                                                                                                    |                                                                                                                                                    aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaaa                                                                                                                                                  >  2:378457/1‑150 (MQ=255)  aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaaaaa                                                                                                                                                <  1:615783/151‑1 (MQ=255)  aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaaaaa                                                                                                                                                <  2:425976/151‑1 (MQ=255)  aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaaaa                                                                                                                                                 >  2:851633/1‑150 (MQ=255)  aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaaaa                                                                                                                                                 >  2:193138/1‑150 (MQ=255)          cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑GTCATATCAAAAACCGc                                                                                                                                            <  2:331012/147‑1 (MQ=255)            tGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCAtt                                                                                                                                        >  2:328622/1‑149 (MQ=255)            ggCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCAtt                                                                                                                                        >  2:13287/2‑149 (MQ=255)                gATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCtc                                                                                                                                   >  2:875863/1‑150 (MQ=255)                 aTGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCtcat                                                                                                                                 >  1:59053/1‑151 (MQ=255)                  tGTTTTGTGATTGTGTTGTTTGACATTAATATGAAAAAAACCATAATCCCATCATGTTTGTGGGGGTTTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGACAAAAAATATGTTGAAAGG‑CTCATATAAAAAACCGCCATTAGCtcat                                                                                                                                 <  2:64171/150‑1 (MQ=255)                    ttttGTGAAGGTGTTGGTTGACATTCATAAAAAAAAAATCATAATCCCATCATTTTTGTGGGGGTTTTTTTCTTTATCTTACCCATCGCACCCAACCCTAAAAAAAAATATGCGGAAAGG‑CTAATATAAAAAACCGCCATTAGCtcatca                                                                                                                               <  1:743742/150‑1 (MQ=255)                         tgAATGTGTTGTTTAACATTCATATAAAAAAAAACATAATCCAATCTTGTTTGTGGGGGTTTTTTTCTGTTTCTTCCCCACCGCACCCAAGCTTAAAAAAAAATATGCTGAAAGG‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAag                                                                                                                          <  1:397813/150‑1 (MQ=255)                           atttttttggttgACTTTAATATGAAAAAAATCATAATCCCATCATGTTTGGGGGGGTTTTTTTCTGTTTCTTCCGCACCCCACTCAAGCCTGACAAAAAATATGCTGAAAGC‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAaga                                                                                                                         <  2:343074/143‑1 (MQ=255)                           attttgtTGGTTGACTTTAATATAAAAAAAATCAAAATCCCATAATGTTTGTGGGGGTTTTTTTTTTTATCTTCCCCACCCCCCCCAAGCCTGAAAAAAAATAGCTTGAAGGG‑CTCATAAAAAAAACCGCCATTAGCTCATCAGGaa                                                                                                                           <  1:388181/143‑1 (MQ=255)                             tgtgTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                       >  1:224150/1‑149 (MQ=255)                                  tggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                                >  2:451459/1‑151 (MQ=255)                                        attttCATAAGAAAAAAAAATTACTCCCATCAGGTTGGTGTGGGTTTTTTTCTTTATCTTCCCCATCGCTCCCAGGCATAACATAAAAAATCCT‑TAAGGCCCCATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCa                                                                                                           <  2:490110/147‑1 (MQ=255)                                        aaaTTCATATAAAAAAAATAATAATCCCATCTTTTTTGGGTGGGTTTTTTTCTGTATCTCACGAATCGCACTAAAGCCTGAAAAAAAATCGGCTGAAAGG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAg                                                                                                          <  1:485868/149‑1 (MQ=255)                                                atGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                                                  >  1:541096/1‑151 (MQ=255)                                                   aaaaaaaTCATAATCCCATCATTTTTGGGGGGGTTTTTTTCTGTATCTTCCGCATCCCACTCAAGCCTGAAAAAAAATATGCTGTAAGC‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                               <  2:721716/151‑1 (MQ=255)                                                   aaaaaaaCAAAAATCCCATCATGTTGTTTTGGGTTCTTTTCTTTATCTCCCCCACCCCCCTCAACCCAGACAAAAATTTTGTTAAAAGG‑CTAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                               <  2:711481/151‑1 (MQ=255)                                                     aaaaacaaaaaTCCCTTCTTGTTTGGGGGGGTTTTTTTTTTTATCTTACGCTTCGCACGGAGCCTTGAAAAAAAATTCTTTGAAAGC‑CCAATACAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                             <  1:51245/142‑1 (MQ=255)                                                      aaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                            >  1:425631/1‑151 (MQ=255)                                                         aTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc                                                                                          >  2:248832/1‑150 (MQ=255)                                                              aaTTCCATCATGTTTGTGTGGGTTCTTTTTTGTATCTTACGCACCCCACTCACGCCTGAAAAAAAATATGCTGAAAGG‑CTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg                                                                                         <  1:118736/146‑1 (MQ=255)                                                              aaTCCCTTCATGTTTGTGGGGGGTTTTTTCTTTACCTACCGCATCGCACTAAGGCATAAAAAAAAATTGGCTGTAGGG‑CTCATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                                      <  2:431860/149‑1 (MQ=255)                                                                      catGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                            >  1:507774/1‑151 (MQ=255)                                                                      catGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                            >  1:519939/1‑151 (MQ=255)                                                                           ttGTGGGGGTTTTTTTCTTTATCTTCCGCACCCCCCACAACCCTAACAAAAAATATGCTGAAAGC‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                        <  1:57981/150‑1 (MQ=255)                                                                              tgtgGGGTCTTTTCTTTATCTTGCGCATCGCACGCAGGCTTGACAAAAAAATTGCTGTAAGG‑CTGATACAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                                     <  2:388295/150‑1 (MQ=255)                                                                                 gggtttttttCTGTATCTTACGCACCGCACTCAACCCTGACAAAAAATATGCTGAAAGC‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCgg                                                                   <  2:425009/143‑1 (MQ=255)                                                                                 ggggTCTTTTCTGTATCTTGCCCCTCCCACTCAGGCGTAAAAAAAAATATGGTGGAAGG‑CTAAAAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCg                                                                    <  1:722453/148‑1 (MQ=255)                                                                                   gtttttttCTTTATCTTACGCACCCCACTCAAGCCTGAAAAAAAATATGCTGTAAGT‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                               <  2:748982/147‑1 (MQ=255)                                                                                     tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                             >  2:232637/1‑151 (MQ=255)                                                                                           cTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTc                                                        >  1:334488/1‑150 (MQ=255)                                                                                                tCTTACGCATCGCCCCCAAGCCTGAAAAAAAATAGGCTGTAAGG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa                                                    <  2:844513/149‑1 (MQ=255)                                                                                                   tACGCACCTCACACAAGCCGGAAAAAAAAGATGCTTTAAGG‑CTCATACAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActct                                                <  1:745242/150‑1 (MQ=255)                                                                                                       caccccacccagcccTGAAAAAAAATATGCTGAAAGC‑CCAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                               <  1:231334/134‑1 (MQ=255)                                                                                                       cATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                            >  2:550722/1‑150 (MQ=255)                                                                                                        aTCCCACACAAGCCTGACAAAAAATATGTTGAAAGG‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                           <  2:224191/150‑1 (MQ=255)                                                                                                         tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt                                         >  2:425824/1‑151 (MQ=255)                                                                                                         tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt                                         >  1:741192/1‑151 (MQ=255)                                                                                                           gCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                        >  2:639045/1‑150 (MQ=255)                                                                                                                 aagccgaaaaaaaaaTATGCGGTAAGG‑CTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                                    <  1:36378/137‑1 (MQ=255)                                                                                                                 aacccTGAAAAAAAATATGCTGTAAGG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa                                     <  2:59062/144‑1 (MQ=255)                                                                                                                  agccaaaCAAAAAATAGGGTGAAAGG‑CTAATAAAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                                    <  1:461161/141‑1 (MQ=255)                                                                                                                  aGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc                                >  1:149044/1‑151 (MQ=255)                                                                                                                        aaaaaaaaTATGCTGAAAGG‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                            <  1:178288/145‑1 (MQ=255)                                                                                                                          aaaaaaTATTCTAAAAGG‑CACATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                        <  2:463368/149‑1 (MQ=255)                                                                                                                          aaaaaaTATGCTGTAAGG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                        <  1:452894/149‑1 (MQ=255)                                                                                                                            aaaaTATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcat                      >  1:119298/1‑151 (MQ=255)                                                                                                                                 atGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg                 >  2:604578/1‑151 (MQ=255)                                                                                                                                 atGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg                 >  2:617705/1‑151 (MQ=255)                                                                                                                                      gTAAGG‑CTTATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaaa             <  2:254734/150‑1 (MQ=255)                                                                                                                                       tAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt           >  1:566513/1‑151 (MQ=255)                                                                                                                                        aaGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg          >  1:510346/1‑151 (MQ=255)                                                                                                                                        aaGG‑CAAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg            <  2:492467/149‑1 (MQ=255)                                                                                                                                             cTCATAACAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGaa       <  1:232598/150‑1 (MQ=255)                                                                                                                                                atatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACAtt   >  2:452326/1‑151 (MQ=255)                                                                                                                                                 tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTg  >  2:311254/1‑151 (MQ=255)                                                                                                                                                 tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTg  >  1:798803/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTG  >  NC_000913/345182‑345475

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