Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,264,632	C→G	7.7%	G263A (GGC→GCC)	PP_1106 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,264,632	0	C	G	7.7%	209.7 / 5.6	65	G263A (GGC→GCC)	PP_1106	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base G (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTGCGTAGCGGGT  >  NC_002947/1264493‑1264780                                                                                                                                            |                                                                                                                                                     gcgcACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTagca                                                                                                                                             <  4:88163/149‑1 (MQ=255)    cACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAgcagca                                                                                                                                          >  7:144948/1‑149 (MQ=255)     aCCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCAt                                                                                                                                         >  7:171140/1‑149 (MQ=255)             cTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTagca                                                                                                                                             >  5:23756/1‑137 (MQ=255)             cTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTagca                                                                                                                                             <  6:23756/137‑1 (MQ=255)             cTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGc                                                                                                                                 <  7:24352/149‑1 (MQ=255)              tCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCt                                                                                                                                >  1:123518/1‑149 (MQ=255)                    gAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGaa                                                                                                                          <  6:361939/149‑1 (MQ=255)                     aGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAg                                                                                                                         >  3:169008/1‑149 (MQ=255)                     aGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAg                                                                                                                         <  2:320703/149‑1 (MQ=255)                         gtgagtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCAt                                                                                                                     <  5:309814/149‑1 (MQ=255)                         gtgagtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCAt                                                                                                                     <  6:18335/149‑1 (MQ=255)                          tgagtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATc                                                                                                                    <  3:371120/149‑1 (MQ=255)                              tgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTg                                                                                                                <  2:123518/149‑1 (MQ=255)                              tgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTg                                                                                                                >  3:348636/1‑149 (MQ=255)                                 ccATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTaa                                                                                                             >  5:282262/1‑149 (MQ=255)                                  cATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAAt                                                                                                            >  7:96000/1‑149 (MQ=255)                                                 aCCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCa                                                                                             >  1:187389/1‑149 (MQ=255)                                                         ggCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTa                                                                                     >  5:137573/1‑149 (MQ=255)                                                          gCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTAc                                                                                    >  5:145836/1‑149 (MQ=255)                                                            cAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGgc                                                                                                <  4:117051/135‑1 (MQ=255)                                                            cAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGgc                                                                                                >  3:117051/1‑135 (MQ=255)                                                            cAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTg                                                                                  >  2:74202/1‑149 (MQ=255)                                                                  cccAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCg                                                                            <  8:96000/149‑1 (MQ=255)                                                                     agagCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCt                                                                         >  1:214995/1‑149 (MQ=255)                                                                           ggccaggccaGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACggg                                                                   >  1:189527/1‑149 (MQ=255)                                                                            gccaggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCTGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACgggg                                                                  <  2:214995/149‑1 (MQ=255)                                                                            gccaggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGc                                                                                                         <  8:233870/110‑1 (MQ=255)                                                                            gccaggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGc                                                                                                         >  7:233870/1‑110 (MQ=255)                                                                            gccaggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGc                                                                                                         <  6:232133/110‑1 (MQ=255)                                                                            gccaggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGc                                                                                                         >  5:232133/1‑110 (MQ=255)                                                                              caggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGgc                                                                <  2:187389/149‑1 (MQ=255)                                                                                   cagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAg                                                           <  4:169008/149‑1 (MQ=255)                                                                                   cagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAg                                                           >  4:182563/1‑149 (MQ=255)                                                                                    agTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAgg                                                          <  6:282262/149‑1 (MQ=255)                                                                                    agTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAg                                                           >  6:152633/1‑148 (MQ=255)                                                                                      tCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGAc                                                        <  8:144948/149‑1 (MQ=255)                                                                                        ggCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCa                                                      <  5:152633/149‑1 (MQ=255)                                                                                           gTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCCCCGTCAGCCGGTGgcag                                                              <  3:74112/138‑3 (MQ=255)                                                                                           gTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCCCCGTCAGCCGGTGgcag                                                              >  4:74112/1‑136 (MQ=255)                                                                                           gTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTg                                                   <  6:137573/149‑1 (MQ=255)                                                                                             ggTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCACCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGcc                                                 >  7:335962/1‑149 (MQ=255)                                                                                               tCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTg                                               >  5:93174/1‑149 (MQ=255)                                                                                                 ggACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCGACTGCCACCGGCTGCAGGGGGCGCCAGGACCAGAGCCTGgc                                             >  4:280460/1‑149 (MQ=255)                                                                                                  gACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGgcg                                            <  2:308038/149‑1 (MQ=255)                                                                                                   aCAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGCACTGGTCCTGgcgc                                           <  3:182563/149‑1 (MQ=255)                                                                                                     aGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCgcag                                         >  4:91352/1‑149 (MQ=255)                                                                                                     aGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCgcag                                         <  3:91352/149‑1 (MQ=255)                                                                                                       cacaAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGc                                                  <  4:160224/138‑1 (MQ=255)                                                                                                       cacaAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGc                                                  >  3:160224/1‑138 (MQ=255)                                                                                                         caAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGt                                     <  6:242037/149‑1 (MQ=255)                                                                                                          aaaGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTg                                    >  6:321252/1‑149 (MQ=255)                                                                                                                ccGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCAt                                                                                                                     <  4:364616/62‑1 (MQ=255)                                                                                                                  ggTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTg                            <  1:74202/149‑1 (MQ=255)                                                                                                                    tCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTgcg                          <  7:95397/149‑1 (MQ=255)                                                                                                                       aCGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTCGCGCAGCAGTGGATCGTTGcgcg                        >  4:330435/1‑148 (MQ=255)                                                                                                                        ctacgcggcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGcagcag                                      >  1:8778/4‑133 (MQ=255)                                                                                                                        ctacgcggcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGcagcag                                      <  2:8778/130‑1 (MQ=255)                                                                                                                           cgcgccgcCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAata                   >  7:23678/1‑149 (MQ=255)                                                                                                                            gcggcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGCCGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAatat                  >  1:366299/1‑149 (MQ=255)                                                                                                                            gcggcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGt                                                    <  3:17395/115‑1 (MQ=255)                                                                                                                            gcggcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGt                                                    >  4:17395/1‑115 (MQ=255)                                                                                                                             cgccgcCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATc                 <  5:321252/149‑1 (MQ=255)                                                                                                                               ccgccgGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGc               >  5:195287/1‑149 (MQ=255)                                                                                                                                  ccgGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTg            <  6:145836/149‑1 (MQ=255)                                                                                                                                            ccgcGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACAAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTGCGTAGCggga  >  6:102697/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                                     GCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCGCCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTGCGTAGCGGGT  >  NC_002947/1264493‑1264780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑