Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,786,207	C→T	6.0%	L302L (CTG→CTA)	ppsA ←	phosphoenolpyruvate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,786,207	0	C	T	6.0%	208.8 / 2.7	67	L302L (CTG→CTA)	ppsA	phosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base T (2/2);  total (32/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 = 8.58e-01

GTATAACGCTCCATGACCTGACCGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCG  >  NC_000913/1786061‑1786355                                                                                                                                                   |                                                                                                                                                     gTATAACGCTCCATGACCTGACCGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAg                                                                                                                                                    >  4:212016/1‑149 (MQ=255)        gCTCCATGACCTGACCGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGc                                                                                                                                             <  1:46907/149‑1 (MQ=255)            cATGACCTGACCGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGCCCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCActt                                                                                                                                         >  7:327709/1‑149 (MQ=255)                cactgaccGCGTGAGCGCACGGTTTCCGGACGCGCCGGCACACTGCACAGTTTACCGGTGTGGCCATCTTTCTCCCACTCAATATCCAGCGGGCGACCGTATTGTGTCTCAATTTGTCCGGCCTGTTTTGCCAGTTCCTGCACttcttc                                                                                                                                     <  3:76212/147‑1 (MQ=255)                acctgaccGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCTATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACttcttc                                                                                                                                     <  7:230530/149‑1 (MQ=255)                   tgaccGCGTGCGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGTCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGtt                                                                                                                                  <  6:245819/149‑1 (MQ=255)                     accGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGtt                                                                                                                                  <  6:164319/147‑1 (MQ=255)                     accGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTg                                                                                                                                 <  4:115384/148‑1 (MQ=255)                     accGCGTGAGCGCACGGTTTCCGGACGCGCCTGAACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTg                                                                                                                                 <  3:316147/148‑1 (MQ=255)                       cgcgTGAGCGCACGGTTTCCGGCCGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGt                                                                                                                               >  2:11913/1‑148 (MQ=255)                         cgTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAg                                                                                                                            >  4:87800/1‑149 (MQ=255)                             aGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCga                                                                                                                         <  4:130413/148‑1 (MQ=255)                               cgcACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGagaa                                                                                                                      <  8:375591/149‑1 (MQ=255)                                gcACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCgag                                                                                                                        >  2:62841/1‑146 (MQ=255)                                gcACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCgag                                                                                                                        >  7:43511/1‑146 (MQ=255)                                 cACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGagaaga                                                                                                                    <  4:40673/149‑1 (MQ=255)                                  aCGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGAt                                                                                                                   <  3:87800/149‑1 (MQ=255)                                  aCGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGAt                                                                                                                   <  1:2380/149‑1 (MQ=255)                                  aCGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGAt                                                                                                                   >  3:64877/1‑149 (MQ=255)                                      tttCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCa                                                                                                               <  4:356344/149‑1 (MQ=255)                                      gttCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCa                                                                                                               <  5:287944/148‑1 (MQ=255)                                          cgtaCGCGCCTGCACAATGAACAGTTTACCGTTGTGGCCATCTTTCTCCTACTAAATATCCATCTGGCTACCGTAGTGTTTATCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCt                                                                                                           <  6:126338/146‑1 (MQ=255)                                                 gcCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTg                                                                                                    >  4:241236/1‑149 (MQ=255)                                                  ccTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGc                                                                                                   <  4:140804/149‑1 (MQ=255)                                                  ccTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGc                                                                                                   <  7:184737/149‑1 (MQ=255)                                                    tGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCgg                                                                                                 >  3:85450/1‑149 (MQ=255)                                                      cacaATGAACAGTTTACCGGTGTGGCCATCTATCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCgg                                                                                                 <  4:64877/147‑1 (MQ=255)                                                       acaATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTAc                                                                                              <  3:241236/149‑1 (MQ=255)                                                            gAACAGTTTACCGGTGTTGCAATCTTTCGCACACTCAATATCCATCGGGCGAACGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCATTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACg                                                                                             <  8:86504/145‑1 (MQ=255)                                                            gAACAGTTTACCGGTGTGGCCATATTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCATGCGGTACGtag                                                                                          >  7:86504/1‑146 (MQ=255)                                                             aaCAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCTAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTc                                                                                        >  3:318203/1‑149 (MQ=255)                                                             aaCAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCTAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTc                                                                                        >  3:323610/1‑149 (MQ=255)                                                             aaCAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTc                                                                                        >  4:100911/1‑149 (MQ=255)                                                                  tttACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGAttt                                                                                   <  1:62841/149‑1 (MQ=255)                                                                  tttACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGAttt                                                                                   >  3:6073/1‑149 (MQ=255)                                                                     aCCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTa                                                                                 <  2:109154/148‑1 (MQ=255)                                                                            tgGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCt                                                                          <  3:91194/148‑1 (MQ=255)                                                                              gCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGc                                                                       >  6:155050/1‑149 (MQ=255)                                                                              gCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTGCTCAATTTGTACGGCCTGTTTTGCCAGTTACTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGc                                                                       >  1:132607/1‑149 (MQ=255)                                                                                cATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCg                                                                     >  3:121769/1‑149 (MQ=255)                                                                                  tCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGcc                                                                      >  6:24408/1‑146 (MQ=255)                                                                                  tCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGcc                                                                      <  5:24408/146‑1 (MQ=255)                                                                                        gCCCACTAAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTAAGGCCTGATTTGCCACTTCCAGCACTTCTTCGTTTTTCAGCGATAATATGTCACGCTGATCCTGAGGTAAGTCTTCGATTTTAACCTGCTTGCCGTCCTCATg                                                             >  8:215665/1‑149 (MQ=255)                                                                                                aaTATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGgcg                                                     >  2:233182/1‑149 (MQ=255)                                                                                                aaTATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGgc                                                      <  5:73994/148‑1 (MQ=255)                                                                                                 atatCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCTAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGgcgc                                                    <  4:318203/149‑1 (MQ=255)                                                                                                 atatCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCTAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGgcgc                                                    <  4:323610/149‑1 (MQ=255)                                                                                                 atatCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGgcgc                                                    >  7:99838/1‑149 (MQ=255)                                                                                                      cATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTaaa                                               >  3:254464/1‑149 (MQ=255)                                                                                                               aCCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGCCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGa                                      >  4:180178/1‑149 (MQ=255)                                                                                                                 cGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGAtt                                    >  8:26497/1‑149 (MQ=255)                                                                                                                 cGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGAtt                                    >  4:234341/1‑149 (MQ=255)                                                                                                                     gtgtTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGAttttt                                 <  8:327709/148‑1 (MQ=255)                                                                                                                      tgtTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGc                                                                                                                           <  8:409254/57‑1 (MQ=255)                                                                                                                      tgtTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGc                                                                                                                           >  7:409254/1‑57 (MQ=255)                                                                                                                            tcAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGA‑TTTTTTCGAccc                         >  7:14440/1‑148 (MQ=255)                                                                                                                              aaTTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGAcc                          <  4:85450/146‑1 (MQ=255)                                                                                                                              aaTTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCa                       >  1:395762/1‑149 (MQ=255)                                                                                                                                     aCGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGc                                                                  <  3:74504/99‑1 (MQ=255)                                                                                                                                     aCGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGc                                                                  >  4:74504/1‑99 (MQ=255)                                                                                                                                      cGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTgcgg               >  5:408146/1‑149 (MQ=255)                                                                                                                                      cGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGCCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTgcg                <  3:180178/148‑1 (MQ=255)                                                                                                                                       ggCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTgcggc              <  2:395762/149‑1 (MQ=255)                                                                                                                                         ccTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTgcggcg             >  3:60971/1‑148 (MQ=255)                                                                                                                                          cTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCa           >  4:119430/1‑149 (MQ=255)                                                                                                                                                tGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATa      <  7:26497/148‑1 (MQ=255)                                                                                                                                                tGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATa      <  3:234341/148‑1 (MQ=255)                                                                                                                                                   cAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCg  >  4:180159/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GTATAACGCTCCATGACCTGACCGCGTGAGCGCACGGTTTCCGGACGCGCCTGCACAATGAACAGTTTACCGGTGTGGCCATCTTTCGCCCACTCAATATCCATCGGGCGACCGTAGTGTTTCTCAATTTGTACGGCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCG  >  NC_000913/1786061‑1786355

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