Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	296,502	A→T	5.9%	L501Q (CTG→CAG)	gshA ←	glutamate‑‑cysteine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	296,502	0	A	T	5.9%	203.1 / 2.0	67	L501Q (CTG→CAG)	gshA	glutamate‑‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (28/35);  new base T (2/2);  total (30/37)
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.29e-01

GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCTCGTCATGCTCGGTC  >  NC_002947/296359‑296648                                                                                                                                                |                                                                                                                                                   gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCg                                                                                                                                               >  6:242314/1‑149 (MQ=255)  ccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCgg                                                                                                                                              <  5:46195/149‑1 (MQ=255)  ccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCg                                                                                                                                               >  2:206116/1‑148 (MQ=255)   cAGGATCGAGGCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGa                                                                                                                                             >  2:242585/1‑149 (MQ=255)     ggATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACt                                                                                                                                           <  6:308238/149‑1 (MQ=255)       aTCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGc                                                                                                                                         <  1:206116/149‑1 (MQ=255)            ctctATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgc                                                                                                                                    <  2:96398/149‑1 (MQ=255)            cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgc                                                                                                                                    >  8:136165/1‑149 (MQ=255)               gATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcga                                                                                                                                  <  2:156815/148‑1 (MQ=255)                        ctggcTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTcgcg                                                                                                                        >  6:219998/1‑149 (MQ=255)                           gctgctGGCGCGTGGTTTGCTCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCgg                                                                                                                       <  6:347055/147‑1 (MQ=255)                           gctgctGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGc                                                                                                                     >  2:242409/1‑149 (MQ=255)                               ctgGCGCGTGGTTTGCGCGGGGTGGGGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCa                                                                                                                   <  5:84332/147‑1 (MQ=255)                                tgGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTc                                                                                                                           <  4:174141/138‑1 (MQ=255)                                tgGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTc                                                                                                                           >  3:174141/1‑138 (MQ=255)                                    gcgTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGtctc                                                                                                            <  7:230014/149‑1 (MQ=255)                                    gcgTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGtctc                                                                                                            >  3:10699/1‑149 (MQ=255)                                             gcgcGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAAtgctt                                                                                                    <  8:165106/148‑1 (MQ=255)                                                  gggTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCTGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTgctgct                                                                                              >  2:334220/1‑149 (MQ=255)                                                       gtggtgCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCg                                                                                         <  6:244491/149‑1 (MQ=255)                                                           tgCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGc                                                                                     <  8:257434/149‑1 (MQ=255)                                                            gCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGc                                                                                     >  8:285354/1‑148 (MQ=255)                                                              ggTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCt                                                                                    >  4:250883/1‑147 (MQ=255)                                                                ttCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGc                                                                                <  6:170798/149‑1 (MQ=255)                                                                  catcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCgg                                                                              >  4:100177/1‑149 (MQ=255)                                                                  catcaACTGCTCATTGCCAGGATGCTGGCCTTGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCg                                                                               <  4:196364/148‑1 (MQ=255)                                                                   atcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGt                                                                             >  6:342800/1‑149 (MQ=255)                                                                    tcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTtgctgc                                                                                               >  3:31618/1‑130 (MQ=255)                                                                    tcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTtgctgc                                                                                               <  4:31618/130‑1 (MQ=255)                                                                    tcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGCTTTGCTCCTGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGCCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGtt                                                                            >  3:193883/1‑149 (MQ=255)                                                                      aaCTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCTGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGTGGTTGc                                                                          <  1:334220/149‑1 (MQ=255)                                                                      aaCTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGc                                                                          >  3:291793/1‑149 (MQ=255)                                                                      aaCTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGc                                                                          <  5:219998/149‑1 (MQ=255)                                                                         tgctgATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTcg                                                                       >  3:118217/1‑149 (MQ=255)                                                                         tgcggATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGACCAGGTCGAAGTCACCGACTTCGTTTTGCTCAAGCTCTGAATGACTGGCGAGCGACTCGCGGGCCAGAGTCGCGAATGCTTGCTGCTTTTCTGTGCTCCGCGGTTGCTct                                                                       >  3:146824/1‑148 (MQ=255)                                                                           ctgATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTcgcg                                                                     <  5:342800/149‑1 (MQ=255)                                                                           ctgATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTcgc                                                                      <  1:242409/148‑1 (MQ=255)                                                                           ctgATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCCCAGCGGTTGCTcgc                                                                      <  7:285354/148‑1 (MQ=255)                                                                            tgATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTcgcg                                                                     >  4:302341/1‑148 (MQ=255)                                                                              aTTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTTGTTGCTCAGCGGTTGCTCGCGGa                                                                   >  5:1433/1‑148 (MQ=255)                                                                               ttGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGaaa                                                                 <  7:136165/149‑1 (MQ=255)                                                                                tGCCAGGATGTTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAg                                                                >  5:278263/1‑149 (MQ=255)                                                                                  ccAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGtc                                                              <  2:328633/149‑1 (MQ=255)                                                                                   cAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCgg                                                                    <  1:274847/142‑1 (MQ=255)                                                                                   cAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCgg                                                                    >  2:274847/1‑142 (MQ=255)                                                                                    aGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTTGTTGCTCAGCGGTTGCTCGCGGAAAGtctc                                                            <  6:1433/149‑1 (MQ=255)                                                                                     ggATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCg                                                           >  4:151731/1‑149 (MQ=255)                                                                                      gATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCgg                                                          <  3:151731/149‑1 (MQ=255)                                                                                       aTGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGc                                                         <  3:250883/149‑1 (MQ=255)                                                                                        tGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCg                                                        >  6:305231/1‑149 (MQ=255)                                                                                          cTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTg                                                      <  7:10846/149‑1 (MQ=255)                                                                                              ccTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGcaca                                                  <  1:242585/149‑1 (MQ=255)                                                                                               cTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCacac                                                 >  2:141066/1‑149 (MQ=255)                                                                                                                aCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGACTGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCgag                                <  5:242314/149‑1 (MQ=255)                                                                                                                    gTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACt                            >  3:143451/1‑149 (MQ=255)                                                                                                                    gTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACt                            >  8:15762/1‑149 (MQ=255)                                                                                                                    gTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGAGGAGTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCTGCGTGCACACGGCTCTGGCGCAGCGAGCAcg                            >  8:313523/1‑148 (MQ=255)                                                                                                                     acGAAGTCGGTGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTg                           >  2:295800/2‑149 (MQ=255)                                                                                                                        aaGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGAc                        <  4:50005/149‑1 (MQ=255)                                                                                                                           tCACCGACTTCGGTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGaa                     <  2:93227/149‑1 (MQ=255)                                                                                                                             aCCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAg                    <  1:141066/148‑1 (MQ=255)                                                                                                                             aCCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAg                    <  3:100177/148‑1 (MQ=255)                                                                                                                                      gTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCTCGTCa           <  7:15762/148‑1 (MQ=255)                                                                                                                                           gCTCCTGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGCCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCTCGTCATGCTc      <  4:193883/148‑1 (MQ=255)                                                                                                                                            cTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCTCGTCATGCTCgg    <  1:295800/149‑1 (MQ=255)                                                                                                                                              ccAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGAc                        >  4:37835/1‑127 (MQ=255)                                                                                                                                               cAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCTCGTCATGCTCGGTc  <  6:278263/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCTCGTCATGCTCGGTC  >  NC_002947/296359‑296648

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: ‑