Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	848,505	C→G	6.1%	R22P (CGG→CCG)	PP_0730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	848,505	0	C	G	6.1%	202.4 / 3.3	66	R22P (CGG→CCG)	PP_0730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

GCTGAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCC  >  NC_002947/848357‑848653                                                                                                                                                     |                                                                                                                                                     gCTGAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGcc                                                                                                                                                      >  7:65444/1‑149 (MQ=255) gCTGAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGcc                                                                                                                                                      >  5:191948/1‑149 (MQ=255)    gAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGtt                                                                                                                                                   >  7:302686/1‑149 (MQ=255)          gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCt                                                                                                                                              <  2:343327/148‑1 (MQ=255)          gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCt                                                                                                                                              <  2:145421/148‑1 (MQ=255)          gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTg                                                                                                                                             <  4:4533/149‑1 (MQ=255)           ccgccTTCCAGATTTTAGATACGCAAATAAGCATCCAAATATCGCAGGTAAGGCTCGTGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            <  7:307107/149‑1 (MQ=255)           ccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            <  5:145989/149‑1 (MQ=255)            cgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGcc                                                                                                                                           <  2:143509/149‑1 (MQ=255)            cgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAACGGCTTGGCTGcc                                                                                                                                           <  6:308414/149‑1 (MQ=255)               cTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            >  8:259632/1‑145 (MQ=255)               cTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            <  7:259632/145‑1 (MQ=255)                ttCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCa                                                                                                                                        >  3:47520/1‑148 (MQ=255)                  ccAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCaa                                                                                                                                       <  8:302686/147‑1 (MQ=255)                      aTTATCGATACGAAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTTATCGAGATCTTCAGGAACTAGGCGTCGCTCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCGAAAAAGc                                                                                                                                   >  4:183885/1‑147 (MQ=255)                            gATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGc                                                                                                                           >  1:155939/1‑149 (MQ=255)                                cGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCt                                                                                                                       >  1:4912/1‑149 (MQ=255)                                    aaTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACg                                                                                                                   <  3:110770/149‑1 (MQ=255)                                    aaTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACg                                                                                                                   <  2:155939/149‑1 (MQ=255)                                     aTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGAGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACgg                                                                                                                  >  5:285281/1‑149 (MQ=255)                                         gCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCt                                                                                                              >  6:177479/1‑149 (MQ=255)                                           aTCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGt                                                                                                            >  7:264804/1‑149 (MQ=255)                                              cacaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATGGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcga                                                                                                         <  8:65444/149‑1 (MQ=255)                                               acaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGt                                                                                                                                                    <  7:51006/105‑1 (MQ=255)                                               acaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGt                                                                                                                                                    >  8:51006/1‑105 (MQ=255)                                               acaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcgac                                                                                                        >  4:196621/1‑149 (MQ=255)                                                         ggTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCTTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACg                                                                                              <  5:104910/149‑1 (MQ=255)                                                               gCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTc                                                                                        <  6:191948/149‑1 (MQ=255)                                                                cTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCa                                                                                       <  8:311314/149‑1 (MQ=255)                                                                 tCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAt                                                                                      >  1:193756/1‑149 (MQ=255)                                                                   gagaCGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGCCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                     >  3:167175/1‑148 (MQ=255)                                                                    agaCGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTg                                                                                    >  5:251568/1‑148 (MQ=255)                                                                          gAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTc                                                                             >  7:87316/1‑149 (MQ=255)                                                                           aGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCg                                                                            <  8:262120/149‑1 (MQ=255)                                                                            ggggtcgtcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCgg                                                                           <  5:48471/146‑1 (MQ=255)                                                                            ggTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTct                                                                            >  5:310034/1‑147 (MQ=255)                                                                                   tcCTTGATCGAGAGCGTCAGGAACGAGGCGTCGCGCTGATGAATGCAGGTGGCGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCa                                                                    <  6:371997/149‑1 (MQ=255)                                                                                   tcCTTGATCAAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTAc                                                                      >  8:274482/1‑147 (MQ=255)                                                                                     cTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAgg                                                                  <  3:196621/149‑1 (MQ=255)                                                                                        taTCGAGATCGTCAGGAACGATGCGTCGCGCTGATGAATGAAGGTGACGTGGCATCCAAGACGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTTTCGACGAAAACCACTTGTTTCATTGGAACCTCGGATTACCAggg                                                                 <  3:183885/146‑1 (MQ=255)                                                                                        gATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTCTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGAt                                                               >  5:32640/1‑149 (MQ=255)                                                                                         aTCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCGCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACCAGTTTCATTGGAACCTCGGCTTACCAGGGAt                                                               >  8:111067/1‑148 (MQ=255)                                                                                          tCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATg                                                              >  4:63718/1‑148 (MQ=255)                                                                                          tCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGa                                                             <  4:47520/149‑1 (MQ=255)                                                                                              gaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGt                                                         <  2:4912/149‑1 (MQ=255)                                                                                              gaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGt                                                         <  1:150033/149‑1 (MQ=255)                                                                                               aTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCGAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCATGGATGAAAGt                                                         <  6:35912/148‑1 (MQ=255)                                                                                                   tCAGGAACGAGGCGGCGCGCTGATGAATGAAGGTGACGGGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATa                                                    <  8:120545/149‑1 (MQ=255)                                                                                                        aaCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCt                                               <  7:274482/149‑1 (MQ=255)                                                                                                            aGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGc                                           >  5:137128/1‑149 (MQ=255)                                                                                                               cggcgctCTGATGAATGAAGGTGCAGTGGCAGCCAAGCCGTTTGGCAGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCg                                        <  2:104333/146‑1 (MQ=255)                                                                                                                gTCGCGCTGATGAATGAAGGTGACGTGGCAGCCACGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCg                                        >  6:20745/1‑148 (MQ=255)                                                                                                                   gcgcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGGCCGaa                                      >  2:133865/1‑147 (MQ=255)                                                                                                                   gcgcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACa                                    <  1:131010/149‑1 (MQ=255)                                                                                                                   gcgcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACa                                    <  2:193756/149‑1 (MQ=255)                                                                                                                             atgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACGTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          <  5:177479/149‑1 (MQ=255)                                                                                                                             atgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          <  8:87316/149‑1 (MQ=255)                                                                                                                              tgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaag                         >  2:363023/1‑149 (MQ=255)                                                                                                                                  ggTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaag                     <  1:363023/149‑1 (MQ=255)                                                                                                                                  ggTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaa                      <  3:116912/148‑1 (MQ=255)                                                                                                                                  ggTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagc                        >  4:116912/1‑146 (MQ=255)                                                                                                                                    tGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGACCGAACAGCATTGCGAAGCAAGCg                   >  3:204742/1‑149 (MQ=255)                                                                                                                                        gTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCg               >  6:156551/1‑149 (MQ=255)                                                                                                                                        gTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCg               >  8:139279/1‑149 (MQ=255)                                                                                                                                          ggCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTc             >  4:83507/1‑149 (MQ=255)                                                                                                                                           gCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTg                                                 >  2:254605/1‑112 (MQ=255)                                                                                                                                           gCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTg                                                 <  1:254605/112‑1 (MQ=255)                                                                                                                                                cAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGt       >  6:329394/1‑149 (MQ=255)                                                                                                                                                     cGTTTGGCTGACTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGcc  <  8:264804/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCTGAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCC  >  NC_002947/848357‑848653

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