Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,006,299	G→A	5.6%	W116* (TGG→TGA)	zapC →	FtsZ stabilizer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,006,299	0	G	A	5.6%	227.1 / 2.6	71	W116* (TGG→TGA)	zapC	FtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base A (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTG  >  NC_000913/1006151‑1006446                                                                                                                                                     |                                                                                                                                                    ttAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTgg                                                                                                                                                     >  6:42217/1‑149 (MQ=255)    aTTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGtt                                                                                                                                                   <  2:232154/148‑1 (MQ=255)        aTCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGc                                                                                                                                               <  6:73725/148‑1 (MQ=255)         tCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGtt                                                                                                                                                   <  3:92641/143‑1 (MQ=255)         tCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGtt                                                                                                                                                   >  4:92641/1‑143 (MQ=255)         tCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGAGATGCGGCCTGCGTTTGGTTGAGCGa                                                                                                                                             >  6:49957/1‑149 (MQ=255)           tAAATAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATa                                                                                                                                           >  5:143567/1‑149 (MQ=255)            aaaGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATAc                                                                                                                                          >  8:147497/1‑149 (MQ=255)             aaGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATa                                                                                                                                           <  1:161420/147‑1 (MQ=255)             aaGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACt                                                                                                                                         <  2:38653/149‑1 (MQ=255)               gAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCa                                                                                                                                       >  3:46590/1‑149 (MQ=255)                 aCAAAAAGCC‑AACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCAc                                                                                                                                      >  7:203631/1‑147 (MQ=255)                     aaaGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCCTGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGGTGAGCGATACTCACGAGCa                                                                                                                                 <  4:123355/149‑1 (MQ=255)                     aaaGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAg                                                                                                                                   >  2:46599/1‑147 (MQ=255)                               tGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTg                                                                                                                       >  2:62730/1‑149 (MQ=255)                                 gTGCTGAATGCACTGGTATCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGc                                                                                                                      <  2:121856/148‑1 (MQ=255)                                 gTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCt                                                                                                                     <  2:88149/149‑1 (MQ=255)                                  tGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTg                                                                                                                    >  4:100208/1‑149 (MQ=255)                                   gCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCAGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTgg                                                                                                                   >  5:124880/1‑149 (MQ=255)                                    cTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTgg                                                                                                                   <  6:143567/148‑1 (MQ=255)                                      gAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGACGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCg                                                                                                                >  6:43702/1‑149 (MQ=255)                                        aTGCACTGGGAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGtt                                                                                                              <  7:101821/149‑1 (MQ=255)                                            ccTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGAGGCGGCCTGCGGTTGGTTGAGCGATACTCACGCGCAGGTCAATTTGCTGGTCGTTGCAt                                                                                                          >  8:53602/2‑149 (MQ=255)                                              tGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCt                                                                                                        >  8:68909/1‑149 (MQ=255)                                                gTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCAGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTgg                                                                                                      >  7:98382/1‑149 (MQ=255)                                                   gCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGa                                                                                                   <  4:46590/149‑1 (MQ=255)                                                      aTTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCGCGAGCAGGTCAATTTGCTg                                                                                                                    <  1:229074/129‑1 (MQ=255)                                                      aTTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCGCGAGCAGGTCAATTTGCTg                                                                                                                    >  2:229074/1‑129 (MQ=255)                                                      aTTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGGGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGaaaa                                                                                                <  4:65934/149‑1 (MQ=255)                                                            tATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGtt                                                                                                              <  8:81414/129‑1 (MQ=255)                                                            tATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGtt                                                                                                              >  7:81414/1‑129 (MQ=255)                                                              tctcAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCAc                                                                                        >  8:231003/1‑149 (MQ=255)                                                                 cAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTAt                                                                                     >  2:102663/1‑149 (MQ=255)                                                                  aaaCCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTAt                                                                                     <  3:129007/148‑1 (MQ=255)                                                                  aaaCCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGACGCACTAt                                                                                     <  6:112642/148‑1 (MQ=255)                                                                   aaCCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGATTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGc                                                                                   >  4:113173/1‑149 (MQ=255)                                                                   aaCCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGATTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGc                                                                                   >  4:116856/1‑149 (MQ=255)                                                                     ccGCAAATGCAGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCTATACTCACGAGCAGGTCAATTTGGTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCct                                                                                 >  5:199068/1‑149 (MQ=255)                                                                         aaaTGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCctgctg                                                                             <  8:203631/149‑1 (MQ=255)                                                                          aaTGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTgg                                                                            >  6:135992/1‑149 (MQ=255)                                                                               cGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACaa                                                                       >  4:88893/1‑149 (MQ=255)                                                                               cGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACaa                                                                       >  8:57730/1‑149 (MQ=255)                                                                                gAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAAc                                                                      >  5:86056/1‑149 (MQ=255)                                                                                   aaGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGATTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCt                                                                   <  3:116856/149‑1 (MQ=255)                                                                                   aaGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGATTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCt                                                                   <  3:113173/149‑1 (MQ=255)                                                                                      cTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCAGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCg                                                                <  8:98382/149‑1 (MQ=255)                                                                                       tGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCgt                                                               >  5:64640/1‑149 (MQ=255)                                                                                          cATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCgttgt                                                            <  5:159677/149‑1 (MQ=255)                                                                                          cATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGGCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCgttgt                                                            <  5:42217/149‑1 (MQ=255)                                                                                           aTTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCgttgtt                                                           <  3:215251/149‑1 (MQ=255)                                                                                            ttttGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTa                                                          >  7:4304/1‑149 (MQ=255)                                                                                             tttGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTAt                                                         >  5:2265/1‑149 (MQ=255)                                                                                                   tCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCgg                                                    <  6:124880/148‑1 (MQ=255)                                                                                                     ccATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACACCCCTGCGTTGTTATTGCGGGTc                                                 >  4:203631/1‑149 (MQ=255)                                                                                                          gTGAAATGTGGGTGCCAATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTcgcgc                                             >  2:33565/1‑148 (MQ=255)                                                                                                           tGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGa                                           <  1:46599/149‑1 (MQ=255)                                                                                                           tGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGa                                           <  1:62730/149‑1 (MQ=255)                                                                                                           tGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGa                                           <  7:68909/149‑1 (MQ=255)                                                                                                            gAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGAt                                          <  6:64640/149‑1 (MQ=255)                                                                                                             aaaTGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATg                                         <  7:147497/149‑1 (MQ=255)                                                                                                                  tgGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGtt                                    >  2:200694/1‑149 (MQ=255)                                                                                                                  tgGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGt                                     <  6:86056/148‑1 (MQ=255)                                                                                                                             gCCGGGCGATGCGGCCTGCGTGTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCa                         >  3:65911/1‑149 (MQ=255)                                                                                                                                     aTGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATc                 <  7:57730/149‑1 (MQ=255)                                                                                                                                      tGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCa                >  6:58620/1‑149 (MQ=255)                                                                                                                                      tGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCa                <  5:43702/149‑1 (MQ=255)                                                                                                                                         ggCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATg              >  5:21393/1‑148 (MQ=255)                                                                                                                                         ggCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATg              >  5:24603/1‑148 (MQ=255)                                                                                                                                           ccTGCGTTTGGGTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCGCAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAAc           <  5:49957/149‑1 (MQ=255)                                                                                                                                            cTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTg                                   >  6:160742/1‑124 (MQ=255)                                                                                                                                            cTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTg                                   <  5:160742/124‑1 (MQ=255)                                                                                                                                               cGTTTGGTTCAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACa       <  2:19068/149‑1 (MQ=255)                                                                                                                                                     gTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTg  >  7:99497/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTGAAATGTGGGTGCCCATGCCGGGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTG  >  NC_000913/1006151‑1006446

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