Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,313,276	C→G	6.4%	intergenic (‑595/+119)	PP_16SG ← / ← msrQ	16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,313,276	0	C	G	6.4%	288.0 / 7.6	90	intergenic (‑595/+119)	PP_16SG/msrQ	16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base G (2/4);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCG  >  NC_002947/5313128‑5313423                                                                                                                                                     |                                                                                                                                                    aTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCgccgcc                                                                                                                                                     >  2:392091/1‑149 (MQ=255)    tCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGTCTGTGAGATCGAGCGccgcgcgg                                                                                                                                                  <  5:306696/149‑2 (MQ=255)    tCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGAGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGccgcccgc                                                                                                                                                  >  3:271742/1‑149 (MQ=255)        tGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                              >  4:57988/1‑149 (MQ=255)        tGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                              >  7:51626/1‑149 (MQ=255)          gAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCgcggcg                                                                                                                                            >  8:262777/1‑149 (MQ=255)           aaGTGCGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAACTGATATATAAGGGGAAATACCGTATCGGGCCTGCCAACTGTTCGCCACAACAGTATCAGCGCCTATGAGATAGAGCGCCGCCAGCcccgcg                                                                                                                                            >  8:341778/1‑148 (MQ=14)             gTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                         <  3:57259/149‑1 (MQ=255)                gggCGCATTATAAGTGGATTCTAGATGGCGTCAACCGTTAATTTCAAGAAAGTGAAATATAAGGTGAACTACCGTATCCGGGCTGACAAGTGTTCGCCACAACAGTTTCATCGCCTATGAGATAGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                       <  6:421908/148‑1 (MQ=14)                  gcgcATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                         >  3:98407/1‑144 (MQ=255)                  gcgcATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                         <  4:98407/144‑1 (MQ=255)                      aTTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                  >  6:93956/1‑147 (MQ=255)                      aTTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                  >  6:14411/1‑147 (MQ=255)                       ttATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                  >  6:107258/1‑146 (MQ=255)                       ttATAAGCGGATTCGCGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGCGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCATAg                                                                                                                               <  7:300478/149‑1 (MQ=14)                           aaGGGGATTCGAGAGGTCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTTTTCGCCACAACAGTTTCAGCGCCTATGAAATCGAGCACCGCCCCCGCGCCGCTCGAACTCACAGACGc                                                                                                                           >  8:176225/1‑149 (MQ=17)                             ggggATTCGAGAGTGAGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCa                                                                                                                         <  6:93730/149‑1 (MQ=255)                             ggggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGc                                                                                                                           >  5:142693/1‑147 (MQ=255)                             ggggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCa                                                                                                                         <  4:271742/149‑1 (MQ=255)                                   tCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCGcgggcg                                                                                                                                               >  6:259247/1‑121 (MQ=255)                                   tCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCGcgggcg                                                                                                                                               <  5:259247/121‑1 (MQ=255)                                     gagaGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAAc                                                                                                                  >  7:293911/1‑148 (MQ=255)                                     gagaGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAAc                                                                                                                  >  8:120564/1‑148 (MQ=255)                                        agGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCAACCGCGCGGAGCTCGATCTCACAGACGCCAGAGAAACAAAg                                                                                                              >  1:135423/1‑149 (MQ=255)                                            cGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                          >  6:50997/1‑149 (MQ=255)                                                aaCCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGcc                                                                                                       <  5:50997/148‑1 (MQ=255)                                                   cTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCctct                                                                                                   >  5:240886/1‑149 (MQ=255)                                                    tttAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCtctc                                                                                                  >  6:267467/1‑149 (MQ=255)                                                      tAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCtctc                                                                                                  <  8:51626/147‑1 (MQ=255)                                                       aaTTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCca                                                                                                <  5:267467/148‑1 (MQ=255)                                                          ttCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGc                                                                                            >  4:210808/1‑149 (MQ=255)                                                             aaGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCAt                                                                                         >  5:351535/1‑149 (MQ=255)                                                              aGAAAGTGAAATATAAGGGGAAATACCGTATCGGGTCTGCCAAGTGTTCGCCACAACAGTTTCATCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCCGAGAAACAAAGGCGGGCCCTCTCCACAGCCAtt                                                                                        <  5:388254/149‑1 (MQ=255)                                                              aGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCAtt                                                                                        <  6:142693/149‑1 (MQ=255)                                                                    tGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCctct                                                                                                   <  1:177969/132‑1 (MQ=255)                                                                    tGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCctct                                                                                                   >  2:177969/1‑132 (MQ=255)                                                                      aaaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                          >  8:20052/1‑91 (MQ=255)                                                                      aaaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtct                                                                                 >  5:412208/1‑148 (MQ=255)                                                                       aaTATAAGTGGAAATACCGTATCGGGGCTGCCAAGTGTTCCCCACAACATTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCa                                                                               <  7:341778/149‑1 (MQ=255)                                                                           taAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATGTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGa                                                                            <  4:79354/148‑1 (MQ=255)                                                                            aaGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGa                                                                            >  1:304630/1‑147 (MQ=255)                                                                              ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcac                                                                                               <  4:404179/126‑1 (MQ=255)                                                                              ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcac                                                                                               >  3:404179/1‑126 (MQ=255)                                                                              ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAgg                                                                        >  4:233114/1‑149 (MQ=255)                                                                                 gAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGccc                                                                     >  7:326231/1‑149 (MQ=255)                                                                                   aaTACCTTATCGTGGCTGCCAAGTGGTCGCCACAACAGGTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGATCGATCTAACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCtt                                                                   <  7:176225/149‑1 (MQ=39)                                                                                    aTACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCGATTGCTATCTCAAGACAGGCCCTTg                                                                  >  4:249909/1‑149 (MQ=255)                                                                                     tACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTgg                                                                 <  5:14411/149‑1 (MQ=255)                                                                                     tACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTgg                                                                 <  5:93956/149‑1 (MQ=255)                                                                                         gTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGc                                                             <  6:420373/149‑1 (MQ=255)                                                                                           aTCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccc                                                           >  7:4687/1‑149 (MQ=255)                                                                                            tCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccc                                                           <  8:4687/148‑1 (MQ=255)                                                                                              ggggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccac                                                        <  6:412208/149‑1 (MQ=255)                                                                                               gggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccacc                                                       >  1:188145/1‑149 (MQ=255)                                                                                               gggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccacc                                                       >  4:154272/1‑149 (MQ=255)                                                                                               gggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccacc                                                       <  2:304630/149‑1 (MQ=255)                                                                                                  cTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccacccca                                                    <  3:210808/149‑1 (MQ=255)                                                                                                       aaGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACtctc                                               <  7:262777/149‑1 (MQ=255)                                                                                                       aaGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACtctc                                               <  7:120564/149‑1 (MQ=255)                                                                                                         gtgtTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCa                                              >  7:295741/1‑148 (MQ=255)                                                                                                          tgtTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCaaa                                            <  1:392091/149‑1 (MQ=255)                                                                                                               gCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGcc                                       >  6:101255/1‑149 (MQ=255)                                                                                                                 cacaACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCaa                                                                              <  1:31500/108‑1 (MQ=255)                                                                                                                 cacaACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCaa                                                                              >  2:31500/1‑108 (MQ=255)                                                                                                                  acaacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTg                                                                  <  1:114919/119‑1 (MQ=255)                                                                                                                  acaacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTg                                                                  >  2:114919/1‑119 (MQ=255)                                                                                                                    aaaaGTTTCAGCGCCTATTAGATCGAGAGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTATCAACCGCCAGTGCTATCTCAAGACAGGCCCGTGGCCGCCCCAGACCACTCTCAAACCGCCCTAcc                                  >  7:149534/4‑149 (MQ=38)                                                                                                                      caGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCCAACCGCCCTACCCt                                >  5:423266/1‑149 (MQ=255)                                                                                                                      caGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCt                                <  6:240886/149‑1 (MQ=255)                                                                                                                       aGTTTCAGCGCCTATGAGATCGTGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTg                               <  2:135423/149‑1 (MQ=255)                                                                                                                       aGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTg                               <  3:57988/149‑1 (MQ=255)                                                                                                                        gTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGc                              >  5:214812/1‑149 (MQ=255)                                                                                                                             aGCTCCTATGAGTTCGAGCGCCGCCCTCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGTGCCCTCTCCACATCCATTTCTATATCAAGACAGGCACTTGGCAGCACCAGCCCACTCTCAAACCGCACTACCCTGCTACCt                         <  8:149534/149‑1 (MQ=255)                                                                                                                               cgcCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTgg                       <  2:188145/149‑1 (MQ=255)                                                                                                                                gcCTATGAGATCGAGCTCCGCCCGCGCGGCGCTCGATATCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTAAACAGCCATTGCTATCTCAAGTCAGGCCCTTCGCAGCCCCACCCCACTCGCAAACCGCCCAACCCTGCTACCAgg                       >  4:122455/1‑148 (MQ=255)                                                                                                                                 ccTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGc                      >  5:166556/1‑148 (MQ=255)                                                                                                                                    gtGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCac                                                              >  4:371770/2‑104 (MQ=255)                                                                                                                                    gtGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAg                                                              <  3:371770/104‑1 (MQ=255)                                                                                                                                       agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGGCCCACCCCACACTCCAACCGCCCTACCCAGCTACCTGGCAACCATg               >  7:41580/1‑149 (MQ=39)                                                                                                                                        gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg              <  6:351535/149‑1 (MQ=255)                                                                                                                                        gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg              <  5:107258/149‑1 (MQ=255)                                                                                                                                        gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg              <  8:293911/149‑1 (MQ=255)                                                                                                                                         aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg              <  8:326231/148‑1 (MQ=255)                                                                                                                                            gAGCTCCGCCCGCGCGGCGCTCGATCTCACAGACGCAACAGAAACAAAGGCGGGCACTCTACACAGCCATTGCTATCTCAAGACAGCCCATTGGCAGCCCCACCCCACTATCAAACCCACCAACAAATCTACCTGGCAACATAGGGAAg          >  8:325824/1‑149 (MQ=255)                                                                                                                                            gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCa                                              >  1:354048/1‑113 (MQ=255)                                                                                                                                            gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCa                                              <  2:354048/113‑1 (MQ=255)                                                                                                                                             aGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCATCCCCACCCCACTCTCAAACCGCCCTa                                    <  1:219014/122‑1 (MQ=255)                                                                                                                                             aGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTa                                    >  2:219014/1‑122 (MQ=255)                                                                                                                                             aGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGc                      >  3:202649/1‑136 (MQ=255)                                                                                                                                             aGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGc                      <  4:202649/136‑1 (MQ=255)                                                                                                                                              gcgcCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCa                     <  8:110914/136‑1 (MQ=255)                                                                                                                                              gcgcCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCa                     >  7:110914/1‑136 (MQ=255)                                                                                                                                               cgccgcCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgcgc       <  3:154272/149‑1 (MQ=255)                                                                                                                                                gccgccCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGcc      >  8:20596/1‑149 (MQ=255)                                                                                                                                                 ccgcgcggGCGGCGCTCGATCTCATAGGCGCCAGAGAAACAAAGGCGGGCCCTGTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCt                         <  8:78358/129‑1 (MQ=39)                                                                                                                                                 ccgcgcggGCGGCGCTCGATCTCATAGGCGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCt                         >  7:78358/1‑129 (MQ=255)                                                                                                                                                     ccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCGGGCAACCTTGGGAAGCGCCcgcg  >  5:397412/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCG  >  NC_002947/5313128‑5313423

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