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

AATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345183‑345477                                                                                                                                                  |                                                                                                                                                       aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaaa                                                                                                                                                    >  1:1085821/1‑150 (MQ=255) aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaaa                                                                                                                                                    >  1:442242/1‑150 (MQ=255)        tcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCGTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCAt                                                                                                                                           >  2:318337/1‑151 (MQ=255)         cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGcc                                                                                                                                             <  2:623139/148‑1 (MQ=255)         cATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGc                                                                                                                                              <  1:957904/147‑1 (MQ=255)          aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCAtt                                                                                                                                          >  2:84514/1‑150 (MQ=255)           tGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCAtt                                                                                                                                          >  1:176321/1‑149 (MQ=255)            gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGc                                                                                                                                       >  2:335809/1‑151 (MQ=255)                        tgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGa                                                                                                                              >  1:202457/1‑148 (MQ=255)                           aTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:540047/1‑150 (MQ=255)                            ttggttgttggACTTTCATATAAAAAAAATAATAATCCCATCTTTTTTGGGTGGGTTTTTTTTTGTATCTCCCGCCCCCCACACACGCCTCAAAAAAAATATGCTGTAGGACTAAAATAAAAAC‑CCGCCATTAGCTCATCAGGa                                                                                                                              <  1:502587/141‑1 (MQ=255)                              ttttggttgACATTAATATAAAAAAAATCATAATTCCATCTTTTTTGTGGGGGGTCTTTTCTGTATCTTCCGCACCCCACCCAAGCCTGACAAAAAATATGTTGTAGGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAg                                                                                                                       <  2:1015848/147‑1 (MQ=255)                                ttgtttgACATTCATATAAAAAAAATCATAATCCCATCATGTTTGTGTGGGGTCTTTTCTTTATCTTACGCACCCCACTCAACCCTAAAAAAAAATATGCTGAAAGCCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                                     <  2:567580/149‑1 (MQ=255)                                  gtttgACATTAATATAAAAAAAATCAAATTCCCAGCTTGTTTGTGGGGGGTTTTTTTTTTATTTACCGCACCGCCCCCACCCCTAACAAAAAATATGTTGGAAGGCT‑ATTTCAAAAATCCGCCATTAGCTCATCAGGAaga                                                                                                                           <  1:714066/139‑1 (MQ=255)                                   gttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACa                                                                                                                 >  1:904031/1‑150 (MQ=255)                                    ttgACATTCATATAAAAAAAATCATAATCCCATCATTTTTGTGTGGGGTCTTTTCTGTATCTTCCGCACCGCACACAAGCCTGAAAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACaa                                                                                                                <  2:271005/150‑1 (MQ=255)                                       aCATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                              >  2:12072/1‑149 (MQ=255)                                       aCATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                              >  1:525358/1‑149 (MQ=255)                                        cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                           >  1:977058/1‑151 (MQ=255)                                        cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                           >  1:871406/1‑151 (MQ=255)                                          ttCATATGAAAAAAATCATAATCCCATCATTTTTGTGTGGGTTCTTTTCTGTATCTTCCGCATCGCCCTCAAGCCTGAAAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAg                                                                                                            <  2:1060524/148‑1 (MQ=255)                                              tataAAAAAAATCATAATCCCATCTTTTTTGTGGGGGGTTTTTTTTGTATCTTCCGCACCGCACTCAAGCCTGAAAGAAAATATGCTGAAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttg                                                                                                       <  2:926682/149‑1 (MQ=255)                                               ataaaaaaaaTCATAATTCCATCATTTTTGTGGGGGGTTTTTTCTGTATCTTACGCATCCCCCTCAAGCCTGACAAAAAATATGTTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                                        <  2:850553/144‑1 (MQ=255)                                               ataaaaaaaaTCATAATTCCATCATGTTTGTGGGGGGTTTTTTCTGTATCTTCCGCATCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                                                    <  2:94806/148‑1 (MQ=255)                                                 gAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                                                   >  1:566698/1‑150 (MQ=255)                                                 aaaaaaaaTCATAATCCCATCATGTTTGTGTGGGTTTTTTTTTGTATCTTCCGCATCACCCCAAAGCCTGACAAAAAATATGCTGTAAGCCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                                                      <  1:786361/146‑1 (MQ=255)                                                  aaaaaaaTAATAATTCCATCATGTTTGTGTGGGTTTTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                                                   <  1:1037755/149‑1 (MQ=255)                                                  aaaaaaaTAATAATCCCATCATGTTTGTGGGGGTTTTTTTCTGTATCTTACCCATCCCACTCAGGCCTGAAAAAAAATATGCGGAAGGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                                 <  1:44054/151‑1 (MQ=255)                                                   aaaaaaTAATATCTCCATCATTTTTGGGTGGGTTTTTTTATGTCTCTTCCCCACCGCCTTAATCCCTGACAAAAAATATTCTTAATGGCTGATTTAAAAAT‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                                <  2:569544/151‑1 (MQ=255)                                                     aaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                              >  2:42514/1‑151 (MQ=255)                                                      aaaTCATAATCCCATCATGTTTGTGTGGGGTTTTTTCTGTATCTTACGCATCGCACTCAAGCCTGAAAAAAAATATGCTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGTCAACCAGGATTGTTGTATGGTa                                                                                             <  2:71442/151‑1 (MQ=255)                                                      aaaGCAAAATCCCATCATTTTTGTGGGGGTTTTTTTCTGTTTCTAACGCACCGCACCCAAGCGTGAAAAAAAATATGGTGTAGGGCTCATATAAAAAG‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTa                                                                                             <  1:954360/151‑1 (MQ=255)                                                            tAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGt                                                                                       >  1:790345/1‑151 (MQ=255)                                                             aaTCCCATCTTTTTTGTGGGGGGTTTTTTCTGTATCTTCCGCACCCCACCAAACCCTAACAAAAAATTTGCTGAAGGCCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                                        <  2:10710/149‑1 (MQ=255)                                                               ttCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                                    >  1:3703/1‑151 (MQ=255)                                                               tcccattttttttGTGGGGGGTTTTTTTTGTATCTTACGCATCCCACTCAAGCCTGAAAAAAAATTTGCTGAAAGGCTAATAAAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTc                                                                                     <  2:696228/140‑1 (MQ=255)                                                                tCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                                    >  1:480035/1‑150 (MQ=255)                                                                 ccATCATGTTTGTGTGGGGTTTTTTCTGTATCTTACCCACCGCACCCAAGCCAGAAAAAAAATATGGTGAAAGGCACATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg                                                                                  <  1:1150859/151‑1 (MQ=255)                                                                    tcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGcc                                                                               >  2:1019513/1‑151 (MQ=255)                                                                      atGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                             >  2:461936/1‑151 (MQ=255)                                                                       tGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                            >  2:661415/1‑151 (MQ=255)                                                                        gTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGCGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                            >  1:1115007/1‑150 (MQ=255)                                                                        gTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                            >  1:729419/1‑150 (MQ=255)                                                                         tttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                          >  1:30840/1‑151 (MQ=255)                                                                            gtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                                       >  2:757947/1‑151 (MQ=255)                                                                            gtgtgGGCTCTTTTCTATTCCTTCAGCTCCGCACCCACGCCTGAGAGAAAGTTTGTCGTGGGGCTCATATAAAAAT‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACCAACAGCATTGTTGTATGGTACTGTGTTCGAGGCCTCGATgga                                                                       <  1:905357/151‑2 (MQ=255)                                                                                ggttttttttCTGTATTTTCCGCACCGCCCCCAAGCCTGACAAAAAATATGCTGAAGGGTTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                          <  1:1097853/138‑1 (MQ=255)                                                                                  gtttttttCTTTATCTTCCGCATCGCACAAAAGCCTAAAAAAAAATAGGCGGAAGGCCACAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                                 <  1:1013361/147‑1 (MQ=255)                                                                                    tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                               >  1:508132/1‑151 (MQ=255)                                                                                      ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt                                                             >  2:204873/1‑151 (MQ=255)                                                                                      ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAAAAAATATGCTGAAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                              <  2:174752/150‑1 (MQ=255)                                                                                      ttttCTGTATCTTACGCACCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                               <  2:294778/149‑1 (MQ=255)                                                                                       tttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCTTATCACAAA‑CCGCCATTAGCTCTCCAGAACGAGCAGACGACAACCAGGATTGTTGCATGGTACGGGGTTCGAGGCCTCGATGCCGGACCAACttt                                                            >  1:150288/1‑151 (MQ=255)                                                                                         tCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTc                                                          >  1:630/1‑151 (MQ=255)                                                                                           tGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCa                                                         >  2:693667/1‑150 (MQ=255)                                                                                            gTATCTTCCGCACCGCCCCCAACCCCAAAAGAAAATATGTTGTAAGGCTAAAATAAAAAC‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTc                                                          <  2:266413/148‑1 (MQ=255)                                                                                              aTCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAc                                                     >  1:832055/1‑151 (MQ=255)                                                                                              aTCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAc                                                     >  2:495347/1‑151 (MQ=255)                                                                                                   aCGCATCCCCCCCAAGCCTAACAAAAAATATCTTGAAGGCCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActct                                                  <  1:84493/149‑1 (MQ=255)                                                                                                    cgattcccacacaAGCTTGACAAAAAATATGTTGAAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActct                                                  <  2:776919/141‑1 (MQ=255)                                                                                                     gCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                              >  1:547214/1‑151 (MQ=255)                                                                                                       accGCACTCAATCCTGAAAGAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             <  1:331937/148‑1 (MQ=255)                                                                                                       aTCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             >  2:583658/1‑150 (MQ=255)                                                                                                         cGCACTCAAGCCAGAAAGAAAATATGCTTTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                              <  1:1082398/147‑1 (MQ=255)                                                                                                           cACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          >  2:85223/1‑149 (MQ=255)                                                                                                           cACTCAAGCCGGACAAAAAATATGTTGTAAGCCTCATATCAAATA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                 <  1:697901/142‑1 (MQ=255)                                                                                                           cACTAAACCCTGAAAAAAAATATGCTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                 <  2:248455/142‑1 (MQ=255)                                                                                                               cAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                                    >  2:306824/1‑151 (MQ=255)                                                                                                                aagccgaaaaaaaaaTAGGCTGAAAGGCTAAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                         <  2:481974/134‑1 (MQ=255)                                                                                                                aaGCCTGAAAAAAAATATGTTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                                      <  2:525650/148‑1 (MQ=255)                                                                                                                   ccTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGCTGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                                >  1:732816/1‑151 (MQ=255)                                                                                                                     tGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                              >  1:222045/1‑151 (MQ=255)                                                                                                                         aaaaaaTATGCTGTAAGCCTAATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatta                            <  2:572541/147‑1 (MQ=255)                                                                                                                         aaaaaaTAGGTTGAAGGACTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                           <  1:965355/148‑1 (MQ=255)                                                                                                                            aaaTATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcatc                       >  2:681008/1‑151 (MQ=255)                                                                                                                            aaaTATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcatc                       >  1:303627/1‑151 (MQ=255)                                                                                                                              atatGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATccc                     >  2:720331/1‑151 (MQ=255)                                                                                                                                  gCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGc                  >  2:372835/1‑150 (MQ=255)                                                                                                                                   cTGTAAGGCTCATTTAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                    <  1:227282/147‑1 (MQ=255)                                                                                                                                     atAAGGCTAATATAAAAAG‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg              <  1:693339/150‑1 (MQ=255)                                                                                                                                       aaGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt             <  2:1086243/150‑1 (MQ=255)                                                                                                                                            tCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGCTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGAGGAAc        >  1:183298/1‑150 (MQ=255)                                                                                                                                                  cAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  2:395528/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                       AATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345183‑345477

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