Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,349,856	T→C	E84E (GAA→GAG)	yciW ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,349,856	0	T	C	94.1%	234.6 / 1.6	68	E84E (GAA→GAG)	yciW	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (27/37);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GGGCATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTG  >  NC_000913/1349711‑1350003                                                                                                                                                  |                                                                                                                                                    gggCATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCAt                                                                                                                                                 >  1:977012/1‑150 (MQ=255)     atatTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCg                                                                                                                                            >  2:1348614/1‑151 (MQ=255)      tatTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCgg                                                                                                                                           >  2:1364583/1‑151 (MQ=255)         tGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTg                                                                                                                                         <  2:1319317/150‑1 (MQ=255)          gCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTg                                                                                                                                         >  1:61244/1‑149 (MQ=255)           cGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGa                                                                                                                                        >  2:1006639/1‑149 (MQ=255)           cGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGCTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGa                                                                                                                                        >  1:505838/1‑149 (MQ=255)             tGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAg                                                                                                                                      >  2:1284948/1‑149 (MQ=255)              gATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGt                                                                                                                                    <  2:1322630/150‑1 (MQ=255)              gATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTc                                                                                                                                   <  2:1274942/151‑1 (MQ=255)                       cAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACgg                                                                                                                          >  1:501435/1‑151 (MQ=255)                         aCCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGTATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGt                                                                                                                        <  2:316567/151‑1 (MQ=255)                          ccGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTc                                                                                                                       <  1:581969/151‑1 (MQ=255)                          ccGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTc                                                                                                                       <  1:990235/151‑1 (MQ=255)                          ccGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTc                                                                                                                       <  2:795997/151‑1 (MQ=255)                               gcttccTTCCCGTGGTGGAACCAGGCCCGGGGAGCGGCGTAAGTCGTGCTGACTAGTAGTTACATAACCGCTTGACGCCGCAACCCGGAGAAATCGTCAGGGCCGACGGCAAGCGCGCATAAGCGTTGAAGGTCTGTGATCGGGTcgcg                                                                                                                    <  1:804801/144‑1 (MQ=255)                               gatTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGt                                                                                                                  <  2:490226/151‑1 (MQ=255)                                   aCTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGc                                                                                                              <  1:1063476/151‑1 (MQ=255)                                    cTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCt                                                                                                             >  2:1301626/1‑151 (MQ=255)                                      ttCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTgg                                                                                                           >  2:1099721/1‑151 (MQ=255)                                      ttCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTgg                                                                                                           >  2:1347893/1‑151 (MQ=255)                                        cccTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTc                                                                                                         <  1:860446/151‑1 (MQ=255)                                               ggAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAAc                                                                                                  >  2:191082/1‑151 (MQ=255)                                               ggAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAAc                                                                                                  <  1:605412/151‑1 (MQ=255)                                                 aaTCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACaa                                                                                                <  1:1343307/151‑1 (MQ=255)                                                  aTCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACaaa                                                                                               <  2:609533/151‑1 (MQ=255)                                                  aTCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACaa                                                                                                >  2:186288/1‑150 (MQ=255)                                                      ggCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATcc                                                                                            >  1:1103784/1‑150 (MQ=255)                                                      ggCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCa                                                                                           <  2:381402/151‑1 (MQ=255)                                                        cccGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACt                                                                                         >  1:710196/1‑151 (MQ=255)                                                        cccGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACt                                                                                         >  2:404210/1‑151 (MQ=255)                                                         ccGGGGAGCGGCGTACGTCGCGCCGAATAGGAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAATGGGTCACGGTAAGCTGGACAGGGAACAAATCCACt                                                                                         >  1:142528/1‑150 (MQ=255)                                                                  ggcgTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAacac                                                                               <  2:753873/151‑1 (MQ=255)                                                                  ggcgTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAacac                                                                               <  2:977062/151‑1 (MQ=255)                                                                  ggcgTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAacac                                                                               <  1:846742/151‑1 (MQ=255)                                                                   gcgTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACAcg                                                                              >  2:301246/1‑151 (MQ=255)                                                                   gcgTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACAcg                                                                              <  2:1162600/151‑1 (MQ=255)                                                                        aGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAg                                                                         <  1:733396/151‑1 (MQ=255)                                                                           cgcgCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGt                                                                      <  1:898423/151‑1 (MQ=255)                                                                           cgcgCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAggg                                                                       <  2:1388041/150‑1 (MQ=255)                                                                                   aTAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCAt                                                              <  2:622414/151‑1 (MQ=255)                                                                                     agtagtTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATg                                                             <  1:525097/150‑1 (MQ=255)                                                                                      gtagtTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGaa                                                           <  2:875710/151‑1 (MQ=255)                                                                                        agtTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGaa                                                           >  1:391755/1‑149 (MQ=255)                                                                                           tACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAaagga                                                       <  1:950406/150‑1 (MQ=255)                                                                                             cATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAaaggaagg                                                    >  2:1200411/1‑151 (MQ=255)                                                                                                    gCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAg                                             <  1:916697/151‑1 (MQ=255)                                                                                                      tGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCg                                           >  2:126998/1‑151 (MQ=255)                                                                                                      tGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCg                                           >  2:1010653/1‑151 (MQ=255)                                                                                                       ggACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGc                                            <  2:63011/149‑1 (MQ=255)                                                                                                       ggACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCg                                           <  1:634828/150‑1 (MQ=255)                                                                                                         aCGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGttt                                        >  1:164658/1‑151 (MQ=255)                                                                                                               cAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCttttt                                  >  2:541828/1‑151 (MQ=255)                                                                                                                aaaCCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTc                                 >  2:1401434/1‑151 (MQ=255)                                                                                                                 aaCCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTc                                 >  1:382377/1‑150 (MQ=255)                                                                                                                 aaCCTGAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGGGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACATGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCATTTCttttt                                  <  1:542399/149‑1 (MQ=255)                                                                                                                      gAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAgg                             <  2:428261/149‑1 (MQ=255)                                                                                                                      gAGCAACCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAgg                             <  2:1228542/149‑1 (MQ=255)                                                                                                                         cAACCGTCAGGGCGGTCGACAAGCGCGCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAGGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGGAGGGGCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCAc                        <  1:176002/151‑1 (MQ=255)                                                                                                                           aCCGTCAGGGCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGt                      <  1:949599/151‑1 (MQ=255)                                                                                                                                    gCCGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGaa             <  1:840304/151‑1 (MQ=255)                                                                                                                                      cGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGGACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGaa             >  2:616275/1‑149 (MQ=255)                                                                                                                                      cGTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAg           <  1:1355670/151‑1 (MQ=255)                                                                                                                                       gTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAgg          >  1:487316/1‑151 (MQ=255)                                                                                                                                       gTCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAgg          <  2:1177393/151‑1 (MQ=255)                                                                                                                                        tCGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCAGGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGACAGGTAGGGGCACAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGaaa            >  1:1229675/1‑148 (MQ=255)                                                                                                                                         cGACAAGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGa        <  1:234466/151‑1 (MQ=255)                                                                                                                                             aaGCGCTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGc    <  1:1339938/151‑1 (MQ=255)                                                                                                                                               gcgcTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTg  <  2:1182213/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GGGCATATTGCGTGATTTGTGCCAACCGATGATTACTTTCCCTGGTGGAATCAGGCCCGGGGAGCGGCGTAAGTCGCGCCGAATAGTAGTTACATAACCGCTGGACGCCGCAAACCTGAGCAACCGTCAGGGCCGTCGACAAGCGTTCATAAGCGGTGAAGGTCTGTGAACGGGTCACGGTAAGCTGGTCAGGGAACAAATCCACTGCAAGAACACGCGCAGGGTCCAGCCATGAAAGGAAGGGGGCCAGCGTTTCTTTTTCAAGGATCACGTCGAGTAGAAAGGGATCGCTG  >  NC_000913/1349711‑1350003

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