Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,166,934	Δ1 bp	5.9%	coding (586/1212 nt)	creA ←	creatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,166,934	0	G	.	5.9%	255.0 / 9.6	84	coding (586/1212 nt)	creA	creatinase
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base . (2/3);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCG  >  NC_002947/4166786‑4167076                                                                                                                                                      |                                                                                                                                               aTCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGt                                                                                                                                               >  1:1076804/1‑150 (MQ=255) aTCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGt                                                                                                                                               >  2:1252735/1‑150 (MQ=255)   ctcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGc                                                                                                                                             <  1:1887208/150‑1 (MQ=255)   ctagCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGc                                                                                                                                             <  1:1175005/147‑1 (MQ=255)    tcGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGcc                                                                                                                                            >  1:1596963/1‑150 (MQ=255)        ccTTGCCCACACGGCGGCTGGTGACAGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGCCTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCt                                                                                                                                        >  2:829540/1‑150 (MQ=255)         cTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTc                                                                                                                                       >  1:1549821/1‑150 (MQ=255)          ttGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCACTGCA‑TGCCAcccag                                                                                                                                      <  1:1300716/150‑4 (MQ=255)          ttGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCg                                                                                                                                      >  1:320576/1‑150 (MQ=255)            gCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTa                                                                                                                                    >  2:2292316/1‑150 (MQ=255)                  acGGCGGCTGGTGAACGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTAGGGGACGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGGC‑TGCAGTGCCACCTCGTACTCggg                                                                                                                              >  2:1273142/1‑150 (MQ=255)                   cgccggcTGGTGACCGGGATGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGc                                                                                                                             <  2:1047472/147‑1 (MQ=255)                        gcTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATATCGCGCACCATCGTCGTTGTGGCA‑TGCCCTGCCCCCTCGTACGCGGGCACCTg                                                                                                                        >  1:1777214/1‑150 (MQ=255)                          tGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCTGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCt                                                                                                                      >  1:162862/1‑150 (MQ=255)                          tGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCt                                                                                                                      >  2:1315548/1‑150 (MQ=255)                          tGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCt                                                                                                                      >  1:1141477/1‑150 (MQ=255)                             tGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCAc                                                                                                                   >  1:1575938/1‑150 (MQ=255)                               aCCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGc                                                                                                                 >  1:1964541/1‑150 (MQ=255)                                ccGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCa                                                                                                                <  2:1596963/150‑1 (MQ=255)                                 cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCACTGCA‑TGCCAcccag                                                                                                                                      <  1:1190198/127‑4 (MQ=255)                                 cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCACTGCA‑TGCCAcccag                                                                                                                                      >  2:1190198/1‑124 (MQ=255)                                 cGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                               >  2:1473677/1‑150 (MQ=255)                                  gggTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAgg                                                                                                              <  1:1889470/150‑1 (MQ=255)                                      tgtgtgCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCAGCTCGTACTCGGGCACCTGCTCACGCAGGGCtt                                                                                                          <  1:168484/150‑1 (MQ=255)                                      tgtgtgCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCtt                                                                                                          <  1:784023/150‑1 (MQ=255)                                            cgcCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCa                                                                                                    >  1:892861/1‑150 (MQ=255)                                            cgcCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCa                                                                                                    <  2:1331283/150‑1 (MQ=255)                                               cATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTg                                                                                                 <  2:1575938/150‑1 (MQ=255)                                                aTCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgc                                                                                                <  1:1368453/150‑1 (MQ=255)                                                 tCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgcc                                                                                               >  1:1466912/1‑150 (MQ=255)                                                 tCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgcc                                                                                               >  2:1642781/1‑150 (MQ=255)                                                 tCGGTGTTGAGCCAGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgcc                                                                                               <  1:1317425/150‑1 (MQ=255)                                                  cGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgccg                                                                                              <  2:1141477/150‑1 (MQ=255)                                                      gtTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGccgcccc                                                                                          <  1:1642781/150‑1 (MQ=255)                                                        tGAGTCCGGACTGGAACCATGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGGTTGGCAATTTCGAGCACCATCGTCTGGCTGGTG‑TGCAGTGCCACCTCTTACTCAGGCTCCTGCTCACGCAGAGCTTCGACCACTGccgcaccgc                                                                                        <  2:230556/150‑1 (MQ=255)                                                           gCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGa                                                                                     <  2:1329870/150‑1 (MQ=255)                                                             ccGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgc                                                                                                >  2:166693/1‑137 (MQ=255)                                                             ccGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgc                                                                                                <  1:166693/137‑1 (MQ=255)                                                                gACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCg                                                                                <  2:320576/150‑1 (MQ=255)                                                                  cTGGAACCAGGTCCAGGTGTCCATTAGTTCTGCCTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCC‑CCGCGTAATCGAGGCTATGATGGCACCCGCTCACGCAGGGCTTCGACCACTGCCACCCCGCAGATGTCGGc                                                                              >  1:2026696/1‑150 (MQ=255)                                                                  cTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTAGGAGAAGCGTTTGGCAAATGCCCGACCCATCGCCTGGGTAGCC‑AGGTGAGATAGGACGTGCTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGc                                                                              >  1:1112879/1‑150 (MQ=255)                                                                   tGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGcc                                                                             >  2:1945627/1‑150 (MQ=255)                                                                    ggAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAggg                                                                                                             >  2:1375331/1‑117 (MQ=255)                                                                    ggAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAggg                                                                                                             <  1:1375331/117‑1 (MQ=255)                                                                     gAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCAc                                                                           >  1:2243504/1‑150 (MQ=255)                                                                     gAACCAGGTCCAGGTGGCCATCAGCTCGGCGTCGGGGACGCGTTGGGCAATTTCGCGGACCATCGTCTGGATGGCG‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGAGGTCCGCCAc                                                                           >  2:640051/1‑150 (MQ=255)                                                                      aaCCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                               <  2:178860/113‑1 (MQ=255)                                                                      aaCCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                                                               >  1:178860/1‑113 (MQ=255)                                                                        ccaggtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCtg                                                                        >  1:1943234/1‑150 (MQ=255)                                                                            gtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCa                                                                    <  1:829540/150‑1 (MQ=255)                                                                             tccaggtGTCCATCAGCTCGGCGTCGGGCAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCTATGCCGCCGATGTCGGCCACCTGTGCAc                                                                   <  1:652271/150‑1 (MQ=255)                                                                              ccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCACCCCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCAcc                                                                  <  1:1310619/150‑1 (MQ=255)                                                                                aggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGGCTGGGTGGCA‑TGCGGTGCCACCGCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGt                                                                >  2:78767/1‑150 (MQ=255)                                                                                aggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGt                                                                >  1:1251701/1‑150 (MQ=255)                                                                                   tgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGAc                                                             >  1:2265146/1‑150 (MQ=255)                                                                                       cATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt                                                         >  1:1116867/1‑150 (MQ=255)                                                                                       cATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt                                                         <  2:1616363/150‑1 (MQ=255)                                                                                             cTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCCTCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGc                                                   <  1:1291321/150‑1 (MQ=255)                                                                                             cTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGc                                                   <  1:1473677/150‑1 (MQ=255)                                                                                              tCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGcc                                                  >  1:2224383/1‑150 (MQ=255)                                                                                               ctgCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCt                                                 <  1:328574/148‑1 (MQ=255)                                                                                               cGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCt                                                 <  1:1268542/150‑1 (MQ=255)                                                                                                       gggAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCg                                         >  2:453546/1‑150 (MQ=255)                                                                                                        ggAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGctcct                                           >  2:1495951/1‑147 (MQ=255)                                                                                                        ggAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGctcct                                           <  1:1495951/147‑1 (MQ=255)                                                                                                              gTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACt                                  >  1:1843966/1‑150 (MQ=255)                                                                                                              gTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACt                                  <  2:2224383/150‑1 (MQ=255)                                                                                                                  ggCAATTTCGCGCATCATCGCCTGGGTGGCA‑TGCAGTGCCCCCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCCCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGAt                              <  1:1273142/150‑1 (MQ=255)                                                                                                                      aTTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATg                          <  2:2243504/150‑1 (MQ=255)                                                                                                                           gcgcACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCAtatgaaa                     >  2:1865434/1‑144 (MQ=255)                                                                                                                             gcACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATAc                   >  1:994201/1‑150 (MQ=255)                                                                                                                               aCCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGc                                       <  2:1658076/128‑1 (MQ=255)                                                                                                                               aCCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGc                                       >  1:1658076/1‑128 (MQ=255)                                                                                                                                 cATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACtt                                 >  2:196826/1‑132 (MQ=255)                                                                                                                                 cATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACtt                                 <  1:196826/132‑1 (MQ=255)                                                                                                                                 cATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCAt               >  2:633788/1‑150 (MQ=255)                                                                                                                                 cATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCAt               <  2:1251701/150‑1 (MQ=255)                                                                                                                                 cATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCAt               <  2:2240883/150‑1 (MQ=255)                                                                                                                                       cgaggTGGCACTGCA‑TGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgc                         >  1:907648/4‑134 (MQ=255)                                                                                                                                       cgaggTGGCACTGCA‑TGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATgc                         <  2:907648/131‑1 (MQ=255)                                                                                                                                       cTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACgg         <  2:2026696/150‑1 (MQ=255)                                                                                                                                       cTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACgg         <  2:1116867/150‑1 (MQ=255)                                                                                                                                        tGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGgc        <  2:1964541/150‑1 (MQ=255)                                                                                                                                              gCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCg  >  1:903522/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                               ATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCG  >  NC_002947/4166786‑4167076

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