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

AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345179‑345477                                                                                                                                                       |                                                                                                                                                      aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATc                                                                                                                                                       <  1:295399/150‑1 (MQ=255)   aaaaTTATCTCATGCAGGTGTTTTGTGATTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGACTGACAGAAAATATGCTGTAAGGCT‑CATATCaaa                                                                                                                                                    <  1:213139/151‑1 (MQ=255)   aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCaa                                                                                                                                                     >  1:578070/1‑150 (MQ=255)    aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCaaaa                                                                                                                                                   <  1:332659/151‑1 (MQ=255)      aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAAc                                                                                                                                                 >  2:268694/1‑151 (MQ=255)             cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCa                                                                                                                                            <  2:454948/149‑1 (MQ=255)              aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCAtt                                                                                                                                          >  2:1056958/1‑150 (MQ=255)                 cAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCt                                                                                                                                      >  1:991351/1‑151 (MQ=255)                   gATGTTTTGTGAATGTGTTGGTTGACATTCATATAAAAAAAATCATAATCCCATCATTTTTGTGTGGGGTTTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATAAAAAACCGCCATTAGCtc                                                                                                                                     <  1:113808/150‑1 (MQ=255)                     tGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCtcatc                                                                                                                                  >  1:938783/1‑151 (MQ=255)                        tttGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAgg                                                                                                                               >  1:122914/1‑151 (MQ=255)                         ttGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAgg                                                                                                                               >  2:223486/1‑150 (MQ=255)                            tgAATGTGTTGGTTGACATTCATATAAAAAAAATCATAATTCCATCATTTTTGTGTGGGGTTTTTTTTGTATCTTACGCACCGCACTCAAGCCTGACAAAAAATATGTTGAAAGTCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAag                                                                                                                            <  2:753740/150‑1 (MQ=255)                              aaTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:550485/1‑151 (MQ=255)                               aTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:1120215/1‑150 (MQ=255)                                  tgttgtTTGACATTCATATAAAAAAACTCATATTCCCATATTGTTTGTGTGGTGTTTTTTCTGTTTCTTACCCACCCCACCCACCCCTGACAAAAAAAATGCTAAAAGGCT‑AATATAAAAAGCCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         <  2:647689/147‑1 (MQ=255)                                   gttggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacg                                                                                                                    >  2:88806/1‑151 (MQ=255)                                      ggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                                  >  2:419761/1‑150 (MQ=255)                                          gACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                              >  1:203683/1‑150 (MQ=255)                                             aTTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                           >  2:794236/1‑150 (MQ=255)                                             aTTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                           >  1:435106/1‑150 (MQ=255)                                                  tatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATC‑TAAGCATCGCACTCCAACCCGGCAGAAAATATGCTTGAAGGCT‑CATTTTCAAAACCGCCCTTAACCCATCAGGAAGAGCAGACGACCAACAGGAttgtt                                                                                                     >  1:543668/1‑150 (MQ=255)                                                      aaaaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                                 >  2:81966/1‑151 (MQ=255)                                                         aaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                              >  1:811513/1‑151 (MQ=255)                                                           aaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc                                                                                            >  1:395242/1‑151 (MQ=255)                                                               aTAATCCCATCTTTTTTGTGTGGGTTTTTTTCTGTATCTTACGCATCGCCCTCAAGCCTGAAAAAAAATATGTTGTAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgg                                                                                          <  1:121963/149‑1 (MQ=255)                                                                 aaTCCCATCATTTTTGTGTGGGTTTTTTTCTGTATCTTACGCATCCCACCAAAGCCGAACAAAAAATATGCTGAAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGt                                                                                       <  2:714319/150‑1 (MQ=255)                                                                  atcccATCGTTTTTGTGTGGGTTTTTTTCTTTATCTTGCACCTCCCACCAACGCCTGACAAAAAATATGCTTAAACGCT‑AATAACAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGt                                                                                       <  1:794364/146‑1 (MQ=255)                                                                          atGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                             >  2:350679/1‑151 (MQ=255)                                                                          atGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                             >  2:462704/1‑151 (MQ=255)                                                                           tttttGGGGGGGGTTTTTTCTTTATCTTACCCATCGCACCCAAGCCGGAAAAAAAATAGGCTGAAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                              <  2:224216/147‑1 (MQ=255)                                                                            gTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                            >  2:22617/1‑150 (MQ=255)                                                                            gTTTGGGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                            >  1:380877/1‑150 (MQ=255)                                                                              ttgttgggggTTTTTTTTTTATCTTACGCATCTCACTCAACCTTGACAAAAAATAGGCTGTAAGCCC‑TATTTAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                          <  2:275393/144‑1 (MQ=255)                                                                                  gtgGGTTCTTTTCTGTATCTTACGCACCGCACTCAAGCCTAACAAAAAATATGCTGAAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                                       <  2:671898/149‑1 (MQ=255)                                                                                    ggggTTTTTTCTGTATCTTCCGCATCGCACTAAAGCCTGACAAAAAATATGCTGTAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                        <  2:757793/146‑1 (MQ=255)                                                                                      ggTTTTTTCTTTATCTTACGCACCGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                                 <  2:726592/151‑1 (MQ=255)                                                                                          ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt                                                             >  2:149663/1‑151 (MQ=255)                                                                                            ttCTGTCTCTCACGCATCGCACTCAAGCCTAAAAAAAAATATTCTGTAAGTCT‑AATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                            <  1:64416/150‑1 (MQ=255)                                                                                               tGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCa                                                         >  2:367728/1‑150 (MQ=255)                                                                                                    cTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                                                   >  2:1077386/1‑151 (MQ=255)                                                                                                     ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActct                                                  >  1:923476/1‑151 (MQ=255)                                                                                                     ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCAGTCCAACTTTTCAAAActct                                                  >  2:329832/1‑151 (MQ=255)                                                                                                       aCGCATCCCACTCAGGCTGAACAAAAAATATGTTGTAAGCTT‑CAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                                                   <  1:999019/148‑1 (MQ=255)                                                                                                        cgcaccgcacTCAACCCGAAAAAAAAATATGCTGAAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActct                                                  <  2:189486/148‑1 (MQ=255)                                                                                                        cgcaccgcacTAAACCCTGAAAAAAAATATGCTGAAAGGCT‑AATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                                                   <  1:269070/147‑1 (MQ=255)                                                                                                        cccaccGCACACAAGCCTGACAAAAAATATGCTGTAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                 <  1:374367/144‑1 (MQ=255)                                                                                                        cGCACCCCACTCAAGCCTAAAAAAAAATATGCTGTAAGCCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                                                    <  1:911037/146‑1 (MQ=255)                                                                                                           aTCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             >  1:854580/1‑150 (MQ=255)                                                                                                           aTCCCACTCAACCCTAAAAAAAAATATGCTGAAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             <  1:874710/150‑1 (MQ=255)                                                                                                             cccACTCAACCCGAACAAAAAATTTGCTGTAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          <  2:986866/149‑1 (MQ=255)                                                                                                             cGCCCTCAAGCCGGAAAAAAAATAGGTTGAAAGGCA‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttc                                                <  2:288025/145‑1 (MQ=255)                                                                                                             cGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          >  1:712849/1‑151 (MQ=255)                                                                                                                  aaaaGCCGAACAGAAAATATGATTGAAGCTA‑AATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa                                       <  1:815584/147‑1 (MQ=255)                                                                                                                    aacccTGAAAAAAAATATGCTGAAAGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa                                       <  2:513526/144‑1 (MQ=255)                                                                                                                    aaGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc                                   >  2:71258/1‑151 (MQ=255)                                                                                                                      gCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCt                                 >  1:753130/1‑151 (MQ=255)                                                                                                                          gACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatt                             >  1:747644/1‑151 (MQ=255)                                                                                                                            aaaaaaaTATG‑TGGCAGGCTACATTTAAAACACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtta                               <  2:291587/144‑1 (MQ=255)                                                                                                                             aaaaaaTATGCTGTAAGGCT‑AATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                           <  1:101640/148‑1 (MQ=255)                                                                                                                             aaaaaaTATGCTGAAAGGCT‑AATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatta                            <  1:448531/147‑1 (MQ=255)                                                                                                                                  atatGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATccc                     >  2:786143/1‑151 (MQ=255)                                                                                                                                   tatGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                    >  2:685780/1‑151 (MQ=255)                                                                                                                                      gctaaaaGCCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGc                  <  2:938639/145‑1 (MQ=255)                                                                                                                                         gTAAGCCT‑ATTATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa                 <  1:620098/148‑1 (MQ=255)                                                                                                                                          tAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg              >  2:255948/1‑150 (MQ=255)                                                                                                                                           aaGGTT‑AATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa                <  1:1057114/147‑1 (MQ=255)                                                                                                                                           aaGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaaa               <  1:667114/148‑1 (MQ=255)                                                                                                                                           aaGGCT‑CATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa                <  1:493884/147‑1 (MQ=255)                                                                                                                                           aaGGCT‑AATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa                <  2:654139/147‑1 (MQ=255)                                                                                                                                            aGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg           >  2:664322/1‑151 (MQ=255)                                                                                                                                            aGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg           >  1:792742/1‑151 (MQ=255)                                                                                                                                              ccT‑CATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGaa         <  2:756040/150‑1 (MQ=255)                                                                                                                                                   ataaCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAAc        <  1:994943/148‑1 (MQ=255)                                                                                                                                                       cAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  1:252979/1‑150 (MQ=255)                                                                                                                                                       cAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  2:370172/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCT‑CATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345179‑345477

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