Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,035,938	C→T	5.5%	V81V (GTG→GTA)	PP_0894 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,035,938	0	C	T	5.5%	220.6 / 2.2	72	V81V (GTG→GTA)	PP_0894	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/38);  new base T (2/2);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

CCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGTTGATGTC  >  NC_002947/1035790‑1036084                                                                                                                                                     |                                                                                                                                                    ccTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGc                                                                                                                                                     <  1:20914/149‑1 (MQ=255)     cTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACg                                                                                                                                                  <  7:183473/148‑1 (MQ=255)        cACCCAGTCATTGAGGTTGTACTAGTTCGTGACCCGGGCAATCTTGCCACAGTGAAAATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGt                                                                                                                                                <  2:381810/147‑1 (MQ=255)             aGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTcc                                                                                                                                         <  1:17423/149‑1 (MQ=255)              gTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTccg                                                                                                                                        >  8:186933/1‑149 (MQ=255)                 aTTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTccgccg                                                                                                                                     >  2:34053/1‑149 (MQ=255)                       gTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCcgcgc                                                                                                                               >  1:278701/1‑149 (MQ=255)                           tagtagTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCt                                                                                                                           <  8:89013/149‑1 (MQ=255)                           tagtagTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCt                                                                                                                           <  4:82700/149‑1 (MQ=255)                           tagtagTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCt                                                                                                                           <  4:192623/149‑1 (MQ=255)                             gtagtTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGTACGGTGAACTCCGCCGCCGCGCGGCTGc                                                                                                                         >  3:44505/1‑149 (MQ=255)                             gtagtTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGTACGGTGAACTCCGCCGCCGCGCGGCTGc                                                                                                                         >  3:39660/1‑149 (MQ=255)                                 ttGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑Tc                                                                                                                    >  3:198495/1‑149 (MQ=255)                                   tgTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCg                                                                                                                   <  8:289129/147‑1 (MQ=255)                                     tGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCg                                                                                                                >  5:12829/1‑149 (MQ=255)                                       aCCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGtt                                                                                                              <  3:20537/149‑1 (MQ=255)                                                 tCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATcac                                                                                                     >  3:172091/1‑148 (MQ=255)                                                     gCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATcaccac                                                                                                  <  7:88024/147‑1 (MQ=255)                                                     gCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATcaccac                                                                                                  <  6:393102/147‑1 (MQ=255)                                                      ccACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGAt                                                                                               <  4:172091/149‑1 (MQ=255)                                                          aGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCG                                                                                           <  4:313404‑M1/149‑1 (MQ=255)                                                           gTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACCATGTCGG                                                                                          >  7:142070‑M1/1‑149 (MQ=255)                                                                 ataGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCAGTCGGCAACCCTTCGTCATCGGCAAGCTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGgc                                                                                    >  1:79281/1‑149 (MQ=255)                                                                  taGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGgcg                                                                                   >  5:332479/1‑149 (MQ=255)                                                                  taGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGgcg                                                                                   <  6:326378/149‑1 (MQ=255)                                                                       gaaCGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTcgc                                                                              >  7:365426/1‑149 (MQ=255)                                                                           gCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCgg                                                                            >  2:94881/1‑147 (MQ=255)                                                                             cccTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTagc                                                                        <  7:186933/149‑1 (MQ=255)                                                                             cccTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTag                                                                         >  6:147777/1‑148 (MQ=255)                                                                             cccTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTag                                                                         >  6:335050/1‑148 (MQ=255)                                                                              ccTGCCGGCAGCACATGGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGttcgccgcgttcat                                                                       >  7:212499‑M1/1‑135 (MQ=255)                                                                                tGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTagcagc                                                                     >  6:325100/1‑149 (MQ=255)                                                                                    gccagcaCATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGtt                                                                 <  5:335050/147‑1 (MQ=255)                                                                                     gcagcaCATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTc                                                                <  6:259482/149‑1 (MQ=255)                                                                                        gcaCATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATc                                                             <  2:229627/149‑1 (MQ=255)                                                                                         cacaTAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCt                                                            >  5:206586/1‑149 (MQ=255)                                                                                            aTAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATc                                                             >  5:45687/1‑145 (MQ=255)                                                                                            aTAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATc                                                             <  6:45687/145‑1 (MQ=255)                                                                                             tAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTc                                                        <  6:332479/149‑1 (MQ=255)                                                                                                         tgtgctgcCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACgcgg                                             <  4:344743/146‑1 (MQ=255)                                                                                                                tCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCGTCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAAcct                                      >  4:165320‑M2/1‑145 (MQ=255)                                                                                                                tCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCGTCACCACGATGTTGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAAcctg                                     <  2:110159‑M2/149‑5 (MQ=255)                                                                                                                  ggCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAAcctggc                                   <  7:267709‑M2/149‑7 (MQ=255)                                                                                                                   gCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGttcgccgcgttcatggaaaagatgaaccgctgctaccgcgagcgcctggcc                                  <  6:12829‑M1/149‑52 (MQ=255)                                                                                                                     aaCCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGttcgccgcgttcatggaaaagatgaaccgctgctaccgcgagcgcctggccga                                <  1:34053‑M1/149‑54 (MQ=255)                                                                                                                     aaCCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGttcgccgcgttcatggaaaagatgaaccgctgctaccgcgagcgcctggccg                                 <  8:212499‑M1/148‑53 (MQ=255)                                                                                                                       cccTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGttcgccgcgttcatggaaaagatgaaccgctgctaccgcgagcgcctggccgaca                              <  8:142070‑M1/149‑56 (MQ=255)                                                                                                                       cccTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGc                              <  4:198495/149‑1 (MQ=255)                                                                                                                         cTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGcc                             >  3:63773/1‑148 (MQ=255)                                                                                                                         cTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGcc                             >  4:26987/1‑148 (MQ=255)                                                                                                                            cGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCAt                          <  2:278701/148‑1 (MQ=255)                                                                                                                              tCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTg                       <  1:94881/149‑1 (MQ=255)                                                                                                                                ggtggCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGgc                     <  3:19999/146‑1 (MQ=255)                                                                                                                                aTCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACgcggcg                                           >  4:56571/1‑127 (MQ=255)                                                                                                                                aTCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAAcctggcc                                  >  4:292801‑M2/1‑129 (MQ=255)                                                                                                                                aTCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAAcctggcc                                  <  3:292801‑M2/136‑8 (MQ=255)                                                                                                                                aTCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGgc                     <  8:267709/149‑1 (MQ=255)                                                                                                                                 tCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGgcg                    >  4:138239/1‑149 (MQ=255)                                                                                                                                   ggCCAGGTATTGGCCATTCACGGTGAACTCCGCCGCCTCGCGGCTGCCA‑TCGGCGTTCTGCATCACCAAGAAGTCGGCCAGGCCCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCCGCGAACCTGCCCATGCCCAACTGGCGCt                  >  3:306402/1‑149 (MQ=255)                                                                                                                                   ggCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGgcgc                   <  8:365426/148‑1 (MQ=255)                                                                                                                                    gCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTc                 >  8:134554/1‑149 (MQ=255)                                                                                                                                    gCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTc                 >  4:343049/1‑149 (MQ=255)                                                                                                                                       aGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGc                              <  3:325432/133‑1 (MQ=255)                                                                                                                                       aGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGc                              >  4:325432/1‑133 (MQ=255)                                                                                                                                          tATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTg           <  7:134554/149‑1 (MQ=255)                                                                                                                                          tATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTg           <  3:26987/149‑1 (MQ=255)                                                                                                                                            ttGGCCATGTACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGt         <  4:44505/149‑1 (MQ=255)                                                                                                                                            ttGGCCATGTACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGt         <  4:39660/149‑1 (MQ=255)                                                                                                                                             tGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGGGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGt         <  2:79281/148‑1 (MQ=255)                                                                                                                                               gCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCCGGCGCTCGCGGTAGCAGCGGGTCATCGTTTCCATGAACGCGGCGAACCTGGCCTTGCCCCTCTGGCGCGCGCCCTGGGTgc      >  2:122852/1‑148 (MQ=255)                                                                                                                                                ccATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGTTGAt     <  3:138239/149‑1 (MQ=255)                                                                                                                                                  aTGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGTTGATGt   >  4:305814/1‑149 (MQ=255)                                                                                                                                                    gCACGGTGAACTCCGCCGCCGCGCGGCTGCCATTCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGTTGATGTc  >  1:47350/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGGCCATGCACGGTGAACTCCGCCGCCGCGCGGCTGCCA‑TCGGCGTTCTGCATCACCACGATGTCGGCCAGGCGCTCGCGGTAGCAGCGGTTCATCTTTTCCATGAACGCGGCGAACCTGGCCTTGCCCATCTGGCGCTCGCCCTGGTTGATGTC  >  NC_002947/1035790‑1036084

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