Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,460,732	G→C	6.8%	H24Q (CAC→CAG)	PP_2153 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,460,732	0	G	C	6.8%	184.5 / 3.8	59	H24Q (CAC→CAG)	PP_2153	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base C (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

CGCCTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTATTCGC  >  NC_002947/2460584‑2460880                                                                                                                                                     |                                                                                                                                                     cgccTTCGGAAATGACCACCGGCAGCGGCTCGCCGCCTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAgg                                                                                                                                                      >  4:172705/1‑149 (MQ=255)  gccTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCACTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGACCTTGTCGCTTAACTGGCGCAACAGgt                                                                                                                                                     >  3:274891/1‑149 (MQ=255)     ttCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGtt                                                                                                                                                  >  5:182629/1‑149 (MQ=255)        ggAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGc                                                                                                                                                <  4:1682/148‑1 (MQ=255)        ggAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGc                                                                                                                                                >  6:210299/1‑148 (MQ=255)        ggAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGc                                                                                                                                                >  3:281975/1‑148 (MQ=255)           aaTGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACt                                                                                                                                            <  3:378490/149‑1 (MQ=255)           aaTGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACACCTGTTGCgcca                                                                                                                                            <  5:242134/149‑4 (MQ=255)            aTGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGTCGCAACACCTGTTGCgccag                                                                                                                                           <  6:45891/149‑5 (MQ=255)                  accGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTc                                                                                                                                     <  5:326760/149‑1 (MQ=255)                   ccGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCg                                                                                                                                    <  7:323457/149‑1 (MQ=255)                      gCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGc                                                                                                                                  >  4:111497/1‑148 (MQ=255)                       cAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGc                                                                                                                                                <  8:129179/133‑1 (MQ=255)                       cAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGc                                                                                                                                                >  7:129179/1‑133 (MQ=255)                             cTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTg                                                                                                                            >  4:92972/1‑147 (MQ=255)                                gCCGACTTCGCCATGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGGCAGGTGGAGc                                                                                                                        <  8:50753/148‑1 (MQ=255)                                    aCTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCt                                                                                                                   <  6:244308/149‑1 (MQ=255)                                    aCTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCt                                                                                                                   <  3:131982/149‑1 (MQ=255)                                      ttCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCt                                                                                                                   >  7:161532/1‑147 (MQ=255)                                      ttCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCt                                                                                                                   >  7:258643/1‑147 (MQ=255)                                        cGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGa                                                                                                                 <  7:35128/147‑1 (MQ=255)                                          ccAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAg                                                                                                             >  4:359614/1‑149 (MQ=255)                                            aaGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGt                                                                                                           >  7:237259/1‑149 (MQ=255)                                                       gCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAgcggc                                                                                                >  1:87087/1‑149 (MQ=255)                                                        cccAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTGGCTTACCTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGAAGCTCGCTGAGCAGGTTGAACAgcggc                                                                                                >  4:242913/1‑148 (MQ=255)                                                         ccAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAgcggcg                                                                                               >  8:52604/1‑148 (MQ=255)                                                         ccAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGa                                                                                              <  3:111497/149‑1 (MQ=255)                                                                 cgcgCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCtt                                                                                      >  7:303833/1‑149 (MQ=255)                                                                              tCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGACCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTgg                                                                         <  4:274891/149‑1 (MQ=255)                                                                                    cgcTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAg                                                                   <  3:172705/149‑1 (MQ=255)                                                                                       ttGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGc                                                                <  8:303833/149‑1 (MQ=255)                                                                                         gTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCt                                                                           >  6:45863/1‑136 (MQ=255)                                                                                         gTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCt                                                                           <  5:45863/136‑1 (MQ=255)                                                                                          ttCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGAGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTc                                                             <  5:391662/149‑1 (MQ=255)                                                                                            cTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCg                                                            <  8:237259/148‑1 (MQ=255)                                                                                               tgCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCcg                                                        <  2:87087/149‑1 (MQ=255)                                                                                               tgCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCACTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTAGCAACACATCTGGTTCAAGGGCTTCGGCcg                                                        >  1:332707/1‑149 (MQ=255)                                                                                                gCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCcgc                                                       >  8:236975/1‑149 (MQ=255)                                                                                                   cGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAAAACATCAGGTTAAAGGGCTTCCGCCGCGCt                                                    >  8:168208/1‑149 (MQ=255)                                                                                                    gCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCt                                                    >  6:278693/1‑148 (MQ=255)                                                                                                    gCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCt                                                    >  1:33645/1‑148 (MQ=255)                                                                                                         aaGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAgggg                                              >  7:90993/1‑149 (MQ=255)                                                                                                              ggcgAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTa                                         >  4:130862/1‑149 (MQ=255)                                                                                                              ggcgAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTa                                         >  4:376905/1‑149 (MQ=255)                                                                                                                    ggTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGaa                                   <  3:92972/149‑1 (MQ=255)                                                                                                                    ggTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGaa                                   <  6:182629/149‑1 (MQ=255)                                                                                                                      tGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAggggg                                             <  6:138386/137‑1 (MQ=255)                                                                                                                      tGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAggggg                                             >  5:138386/1‑137 (MQ=255)                                                                                                                         ggTCCTTGTCGCTTAACTGGCGAAACAGGTGTTGCGACTCGAAGTCCGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGAAAAACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGAGc                               >  5:327126/1‑148 (MQ=255)                                                                                                                         ggTCCTTGTCACTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTAGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcg                              <  2:332707/149‑1 (MQ=255)                                                                                                                                 tCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGcgatccgatc                      <  8:371365/149‑1 (MQ=255)                                                                                                                                        aCTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATg               <  3:359614/149‑1 (MQ=255)                                                                                                                                         cTGGCGCACCAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGc              <  4:281975/149‑1 (MQ=255)                                                                                                                                          tGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGTGCTTCGGCCGCGCTAAGGCGGCTTATTTGAAGTGCGATCCGATCTTCGATGCg             <  5:278693/149‑1 (MQ=255)                                                                                                                                          tGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCg             >  1:102690/1‑149 (MQ=255)                                                                                                                                           gtcgcAACACCTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCgg            >  4:353126/3‑149 (MQ=255)                                                                                                                                            tcgcAACACCTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCACCACATCTGGTTCAAGGGCTTAGGCCGCGCTAAGTGGGCTTATTTGAAGTGCGATCAGATCTTCGATGCgg            >  1:192714/2‑148 (MQ=255)                                                                                                                                            gcgcAACAGGTGTTGCGACTCGAAGGAGTCCAGTTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGgt           <  7:168208/149‑1 (MQ=255)                                                                                                                                             cgcAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGgta          <  7:52604/149‑1 (MQ=255)                                                                                                                                              gcAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGgtag         <  4:397857/149‑1 (MQ=255)                                                                                                                                               cAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGgtagt        <  5:210299/149‑1 (MQ=255)                                                                                                                                                     gtgtTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTATTcgc  <  3:130862/149‑1 (MQ=255)                                                                                                                                                     gtgtTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTATTcgc  >  4:86732/1‑149 (MQ=255)                                                                                                                                                     gtgtTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTATTcgc  >  8:379312/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGCCTTCGGAAATGACCACCGGCAGCGGCTCGCCGACTTCGCCAAGCAGGGTTTGCCCAACCACCGCGCTGAGCAAGTCGAGGCGCTTGTTCTGTGCCCGCAGGAAGGCGGCGAGGGTGCGGTCCTTGTCGCTTAACTGGCGCAACAGGTGTTGCGACTCGAAGTCGGCCAGGTGCAGCTCGCTGAGCAGGTTGAACAGCGGCGAATCATCTTGCAACACATCTGGTTCAAGGGCTTCGGCCGCGCTAAGGGGGCTTATTTGAAGTGCGATCCGATCTTCGATGCGGTAGTATTCGC  >  NC_002947/2460584‑2460880

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