Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,254	A→T	6.7%	intergenic (+32/‑115)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,254	0	A	T	6.7%	254.9 / 6.4	85	intergenic (+32/‑115)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base T (4/3);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTT  >  NC_000913/781106‑781402                                                                                                                                                     |                                                                                                                                                     ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGTCCCACCAGTTTTAACATAAATCTCTGTTGTTAAGCTTGat                                                                                                                                                     >  2:1457398/1‑149 (MQ=14) ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                                     <  1:484150/150‑1 (MQ=255) ggaggGGGTCGGCGGTGCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                                     <  1:484112/150‑1 (MQ=255)       ggTCGGCGGTTCGATCCCGTCACCACCCACC‑CCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAc                                                                                                                                              <  2:594311/150‑1 (MQ=255)        gTCGGCGGTTCGATCCCGGCCACACCCCCCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAAc                                                                                                                                              <  2:1177967/150‑1 (MQ=255)         tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGTTCGAATACTGCACGACCCACCAATTTTAACATCAAAAACAGATGTTAAGCGTGAAGGATAACg                                                                                                                                             >  1:525172/1‑150 (MQ=11)            gcggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAACTCAGATGAGAGGCGTGAAGGATGACGTTg                                                                                                                                          >  2:1690727/1‑150 (MQ=14)               gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcgc                                                                                                                                       <  1:224013/150‑1 (MQ=255)                ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGcc                                                                                                                                      >  1:847605/1‑150 (MQ=255)                ttCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGcc                                                                                                                                      <  1:975874/150‑1 (MQ=255)                    aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTGAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCCGTTTTAACATCAAACTCAGATGTTAAGCGGGAAGGATAACGTTGCGCCAGCa                                                                                                                                  >  1:1309464/1‑150 (MQ=255)                    aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCa                                                                                                                                  >  1:1259304/1‑150 (MQ=255)                       ccGTCATCACCCACCACTCGGGTCGTTAGCTCATTTGGTAGAGCTGTTGAATTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACg                                                                                                                               <  1:1452028/150‑1 (MQ=255)                          tcatcaCCCACCAC‑CGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGccc                                                                                                                           >  2:1002698/1‑150 (MQ=12)                              cacccaccACGCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTACTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTCACATCCCCCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCCGCAACGGCCCGTa                                                                                                                        <  2:1610152/150‑1 (MQ=14)                                  caccacTCGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGAATTTTAATCAAATTGTCGCAGGTTAGAATACAGCACCACCCACCAGGTTTAACATCAAAAACAGATGTTAAGCGTGAAGGATAACCTTGCGCCAGCAACGGCCCGTAGGgc                                                                                                                    >  1:204647/1‑150 (MQ=14)                                      acTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgagc                                                                                                                >  2:977532/1‑150 (MQ=255)                                             tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                         <  2:155549/150‑1 (MQ=255)                                             tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                         >  2:132000/1‑150 (MQ=255)                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                              >  1:1627372/1‑150 (MQ=255)                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATAGGTCGCACGTACGAATACAGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCGACGGCCCGGAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                              >  1:52096/1‑150 (MQ=14)                                                         ggggTAGCGCAGTTGACTTTTAAACATTTGGCCGCAGGTCCTAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTATGGGTTACGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                                             <  1:1002698/148‑1 (MQ=14)                                                           ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACATCATAAACAGAAGATAACCGTAAAGAAAAACGTTCCTAAAACAACGGACCATAGGCCGAGCGAAGCAATTCATCCTGCaaca                                                                                           >  2:1014955/1‑147 (MQ=2)                                                                agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAGGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                                      >  1:698016/1‑150 (MQ=255)                                                                agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                                      >  1:1340995/1‑150 (MQ=255)                                                                      ttaattttAATCATTGGGTCGCTGTTTCTAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATTTTAAGCGTTATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCt                                                                                <  1:74059/146‑1 (MQ=14)                                                                      tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATAAAACTCAGATGTTAATCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCt                                                                                >  2:917009/1‑150 (MQ=255)                                                                        aCTTTTAATCACTTGGTCGCAGGTTCGAATCCTGCACGACCCACCCGTTTTAACATCAAACTCAGATGTTTAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTa                                                                              <  2:794895/150‑1 (MQ=255)                                                                          ttttAATGTTTTGGTCGTTGGTTGGATTCCTGCACGACCCACCTGTTTTAACATCAAACTCTGATGTTATGGGTGATGGATAAGGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                            <  1:1168536/150‑1 (MQ=17)                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGGGTGAAGGAAAACGTTGCGCAAGAAACGGCCCGTAGGGCGAGCGAAGCGAGTCATACTTCACGAACCACCAATCATaaa                                                                            >  2:1556841/1‑150 (MQ=14)                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCCCCAATCATaaa                                                                            >  2:1036618/1‑150 (MQ=255)                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCATaaa                                                                            >  2:1469347/1‑150 (MQ=255)                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCAACCTGCACGACCCACCAATCTTaaa                                                                            >  2:1469257/1‑150 (MQ=255)                                                                               aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                       >  2:860699/1‑150 (MQ=255)                                                                               aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                       >  2:1077131/1‑150 (MQ=255)                                                                               aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGAAGATAAGCGGGAAGGATAACGTTCCGCCAGCAACGGCCCGTAGGCCTAGCGAAGCGAGTCATCCTGCACGACCCGCCAATCTTaaattttg                                                                       >  2:1077241/1‑145 (MQ=17)                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                              >  1:270937/1‑126 (MQ=25)                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                              <  2:270937/126‑1 (MQ=255)                                                                                    tgggacGCAGGGTCGAAGCCTGCACGCCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcccc                                                                  <  1:49891/145‑1 (MQ=255)                                                                                      ggTCGCTGGTTCTAATCCTGCACGACCCACATGTTTTAACATAAAACTCAGATTTTAAGCGTGAATGTTAAGGTTGCGCCTGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                <  1:602933/150‑1 (MQ=18)                                                                                      ggTCGCTGGTTCGAATCCTGCACGACCCACCATTTTTAACATCAAACTCAGATGTTATGCGTGATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                <  1:638489/150‑1 (MQ=255)                                                                                      ggTCGCTGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCTTATGTTAAGCGTGATGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                <  1:602927/150‑1 (MQ=255)                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                              >  2:252003/1‑150 (MQ=255)                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAAAATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGAGAGTCATCCTGCACGACCCACCAATCTTAAAGATTTGCCCCGatt                                                              >  2:321947/1‑148 (MQ=255)                                                                                           cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTAAAAGATTGGCCCCGAGTaaa                                                           >  2:1699220/1‑150 (MQ=255)                                                                                              gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                                        <  1:977532/150‑1 (MQ=255)                                                                                              gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                                        >  1:1439913/1‑150 (MQ=255)                                                                                               ttCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATc                                                       >  1:396882/1‑150 (MQ=255)                                                                                                tCGAATCCTGCACGACCCACCAGTTTTAACATCCAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                                      <  2:46697/150‑1 (MQ=255)                                                                                                tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCt                                                      >  1:783183/1‑150 (MQ=255)                                                                                                      ccTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                                <  1:252003/150‑1 (MQ=255)                                                                                                         gCACGACCCACCAGTTTTAACATCTAAGTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                             <  2:1676249/150‑1 (MQ=255)                                                                                                              aCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAccc                                        <  2:900942/150‑1 (MQ=255)                                                                                                                   ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGACTAAAAAACTTTCAGGTAACACCCGGATg                                   >  2:1083765/1‑150 (MQ=255)                                                                                                                    caGTTTTAACATCAATCTCCTATTTTATGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTTGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                                  <  1:425568/150‑1 (MQ=25)                                                                                                                    caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTAAAAGATTTGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                                  >  2:294787/1‑150 (MQ=255)                                                                                                                      gTTTTAACATCTAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                <  2:653112/150‑1 (MQ=255)                                                                                                                      gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                <  2:1220159/150‑1 (MQ=255)                                                                                                                      gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                <  2:350571/150‑1 (MQ=255)                                                                                                                       ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTAAAATATAGCCCCCGAGTAGATATATTACAGGTAAAGCCCGAAAGGGTc                               >  1:347357/1‑150 (MQ=17)                                                                                                                       ttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTAAAAGAGAGGCCCCGAGTCAAAATCTTTCAGCTCACACCCGGAAGGGTc                               >  1:347352/1‑150 (MQ=18)                                                                                                                        tttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCg                              <  1:132000/150‑1 (MQ=255)                                                                                                                         ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGGAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGAGTCGt                             >  1:1605758/1‑150 (MQ=255)                                                                                                                              aTCAAACTGAGATGTTACGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                        <  2:1701/150‑1 (MQ=255)                                                                                                                              aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                        <  2:847605/150‑1 (MQ=255)                                                                                                                              aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                        >  2:1552775/1‑150 (MQ=255)                                                                                                                              aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCt                        >  2:1552707/1‑150 (MQ=255)                                                                                                                                 gaagacAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTCAAGTTTGGCCCCGAGTAAAAATCTTTAAGGTAACACCCGAGTGTGTCGTGATCTCAg                     >  2:19083/6‑150 (MQ=17)                                                                                                                                 gaagacAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATATTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCCTATGGCTCGTTAGATCAg                     >  2:383846/6‑150 (MQ=18)                                                                                                                                  gagtttgATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGAAACACCCGTATGGGTCGTTAGCTCAGt                    >  1:1323949/7‑150 (MQ=255)                                                                                                                                  gagtttgATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGAAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt                    >  1:1629484/7‑150 (MQ=255)                                                                                                                                     tCAGATGTTAAGCGTGAATGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGATCGAATCGAGTCATCCTGCACGACCCACCAATCATAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg                 <  2:996043/150‑1 (MQ=255)                                                                                                                                     tCAGATGTTAAGCGTGAAGGATAACGTTGCTCCAGCAACGGCCCGTAGGGCGAGCGAATCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg                 <  2:582276/150‑1 (MQ=255)                                                                                                                                     tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTgg                 >  2:1504885/1‑150 (MQ=255)                                                                                                                                      cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGACCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGt                <  2:1420091/150‑1 (MQ=255)                                                                                                                                       aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa               >  2:1090005/1‑150 (MQ=255)                                                                                                                                       aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa               >  2:1089985/1‑150 (MQ=255)                                                                                                                                       aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa               <  1:614705/150‑1 (MQ=255)                                                                                                                                          tGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                            >  2:726242/1‑86 (MQ=12)                                                                                                                                          tGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                            <  1:726242/86‑1 (MQ=12)                                                                                                                                          tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTagag            >  1:9989/1‑150 (MQ=255)                                                                                                                                            ttAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa          >  1:1678494/1‑150 (MQ=255)                                                                                                                                            ttAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa          <  1:19083/150‑1 (MQ=255)                                                                                                                                            ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa          <  1:383846/150‑1 (MQ=255)                                                                                                                                                    gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt  <  1:1552775/150‑1 (MQ=255)                                                                                                                                                    gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt  <  2:1627372/150‑1 (MQ=255)                                                                                                                                                    gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt  <  1:1552707/150‑1 (MQ=255)                                                                                                                                                    gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt  >  1:1707783/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTT  >  NC_000913/781106‑781402

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