Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,509,064	A→G	5.9%	C32R (TGC→CGC)	PP_2203 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,509,064	0	A	G	5.9%	197.6 / 2.5	68	C32R (TGC→CGC)	PP_2203	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/36);  new base G (2/2);  total (30/38)
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.79e-01

CGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCAGGCGGATACCGGCGACCAGGCGGGTCTGGC  >  NC_002947/2508920‑2509212                                                                                                                                                  |                                                                                                                                                     cgccgATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTgg                                                                                                                                                 <  2:1612269/150‑1 (MQ=255)  gccgATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACTGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTggg                                                                                                                                                <  2:1888105/150‑1 (MQ=255)   ccgATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTgggg                                                                                                                                               <  2:269447/150‑1 (MQ=255)       tCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGATAGGGGCGTTGACGGTAACACAACACATTATTGA‑GCGGCTGGTTGAGGCCTTGCAGCAGCCGCAGTTGGGGTcgc                                                                                                                                           >  1:434905/1‑150 (MQ=255)          ccccGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGt                                                                                                                                        <  2:1421490/150‑1 (MQ=255)           cccGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTTTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCGGTTGGGGTCGCGGTa                                                                                                                                       >  2:1919328/1‑150 (MQ=255)           cccGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTa                                                                                                                                       >  1:1404353/1‑150 (MQ=255)           cccGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTa                                                                                                                                       >  2:1710067/1‑150 (MQ=255)           cccGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTCTCACTCACCTGGTTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTa                                                                                                                                       <  1:37659/150‑1 (MQ=255)            ccGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTTTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCGGTTGGGGTCGCGGTAg                                                                                                                                      <  1:1919328/150‑1 (MQ=255)                    gCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGAt                                                                                                                              >  1:1227622/1‑150 (MQ=255)                             aaGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTcggc                                                                                                                     >  2:1657538/1‑150 (MQ=255)                              aGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCG‑TTGTCACGGAGGTGGCTGATCCAGGCCTTGGAGCAGCCGCAGTCAGGGTCGAACTCAACGTCGATGGTCTcggcg                                                                                                                    >  1:1069185/1‑150 (MQ=255)                                          aCGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAgg                                                                                                        >  2:882583/1‑150 (MQ=255)                                           cGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGc                                                                                                       <  2:1540436/150‑1 (MQ=255)                                               gAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGt                                                                                                   <  2:1758538/150‑1 (MQ=255)                                               gAGGCGGTGCTGGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCCGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGt                                                                                                   <  2:198489/150‑1 (MQ=255)                                                    tttGCTTGTCGGCGCGCGTGTTTGGCTCGTCGTTGTCGTTGACGGTAAAACCT‑TTGTCACGCAGGTGGCTGCCCCAGGGCGTGGAGCAGCCGCACTTGGGGTCGCGGTAGACGTCTATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACgg                                                                                              <  2:410728/149‑1 (MQ=255)                                                           gACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGt                                                                                                                                  >  1:1184642/1‑107 (MQ=255)                                                           gACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGt                                                                                                                                  <  2:1184642/107‑1 (MQ=255)                                                            aCGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAgg                                                                                      <  1:993507/150‑1 (MQ=255)                                                               gcgcTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTc                                                                                   >  2:740204/1‑150 (MQ=255)                                                                  cTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAgg                                                                                <  1:740204/150‑1 (MQ=255)                                                                    cATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGt                                                                              <  1:1277694/150‑1 (MQ=255)                                                                          tGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTcagcag                                                                        <  2:1227622/150‑1 (MQ=255)                                                                            gCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGgt                                                                      >  1:644888/1‑150 (MQ=255)                                                                              tCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGgtgt                                                                    >  1:457423/1‑150 (MQ=255)                                                                                gACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGttt                                                                  >  1:733349/1‑150 (MQ=255)                                                                                         gTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTggg                                                                                                                                                <  2:1822270/63‑1 (MQ=255)                                                                                         gTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTggg                                                                                                                                                >  1:1822270/1‑63 (MQ=255)                                                                                         gTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGt                                                                                                   >  1:1603528/1‑108 (MQ=255)                                                                                         gTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGt                                                                                                   <  2:1603528/108‑1 (MQ=255)                                                                                            gACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGt                                                      <  2:1263856/150‑1 (MQ=255)                                                                                            gACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGt                                                      >  2:753231/1‑150 (MQ=255)                                                                                               ggTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATg                                                        >  1:1776230/1‑145 (MQ=255)                                                                                               ggTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATg                                                        <  2:1776230/145‑1 (MQ=255)                                                                                               ggTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGaa                                                   <  2:188652/150‑1 (MQ=255)                                                                                               ggTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGaa                                                   <  2:1706460/150‑1 (MQ=255)                                                                                               ggTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGaa                                                   <  2:457423/150‑1 (MQ=255)                                                                                                gTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAg                                                  <  2:434905/150‑1 (MQ=255)                                                                                                 tAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAg                                                                                 <  1:400276/118‑1 (MQ=255)                                                                                                 tAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAg                                                                                 >  2:400276/1‑118 (MQ=255)                                                                                                   aaaCCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTcagcag                                                                        <  2:638397/125‑1 (MQ=255)                                                                                                   aaaCCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTcagcag                                                                        >  1:638397/1‑125 (MQ=255)                                                                                                    aaCCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTGTCGGTGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCGCAGGGTCAGCGGGTGTTTGCGCAGCATGGTGAAGCACt                                              >  1:619258/1‑150 (MQ=255)                                                                                                    aaCCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACt                                              >  2:1121176/1‑150 (MQ=255)                                                                                                     aCCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGTCCGGTGACGGGCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTc                                             <  2:1726029/150‑1 (MQ=255)                                                                                                     aCCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTc                                             <  1:1092963/150‑1 (MQ=255)                                                                                                      ccA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTcc                                            <  2:644888/150‑1 (MQ=255)                                                                                                       cA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGTCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCt                                           <  2:1404353/150‑1 (MQ=255)                                                                                                            gTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTaa                                       <  2:1069185/150‑1 (MQ=255)                                                                                                             tCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGGCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTaaa                                      >  2:1289782/1‑150 (MQ=255)                                                                                                              cACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTaaaa                                     <  2:482103/150‑1 (MQ=255)                                                                                                                 gCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGgt                                                                      >  1:921712/1‑114 (MQ=255)                                                                                                                 gCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGgt                                                                      <  2:921712/114‑1 (MQ=255)                                                                                                                 gCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATg                                                        >  2:1551461/1‑128 (MQ=255)                                                                                                                 gCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATg                                                        <  1:1551461/128‑1 (MQ=255)                                                                                                                 gCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGaa                                                   >  1:895574/1‑133 (MQ=255)                                                                                                                 gCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGaa                                                   <  2:895574/133‑1 (MQ=255)                                                                                                                 gCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACg                                  >  2:183978/1‑150 (MQ=255)                                                                                                                     gTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCAgg                              <  1:1121176/150‑1 (MQ=255)                                                                                                                                aGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCAGGCGGATACCGGc                   >  1:955009/1‑150 (MQ=255)                                                                                                                                 ggCCTTGCAGCAGCCGCGGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCa                                                                                                          >  2:445583/1‑62 (MQ=255)                                                                                                                                 ggCCTTGCAGCAGCCGCGGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCa                                                                                                          <  1:445583/62‑1 (MQ=255)                                                                                                                                             gcCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCa                                >  1:531135/1‑124 (MQ=255)                                                                                                                                             gcCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCa                                <  2:531135/124‑1 (MQ=255)                                                                                                                                                gCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCAGGCGGATACCGGCGACCAGGCGGGTCTgg   >  1:482101/1‑150 (MQ=255)                                                                                                                                                 cAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCAGGCGGATACCGGCGACCAGGCGGGTCTGGc  <  1:1348884/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACGGCGCTCATGTTTGGCTCGACGTGGTCGTTGACGGTAAAACCA‑TTGTCACGCAGGTGGCTGATCCAGGCCTTGCAGCAGCCGCAGTTGGGGTCGCGGTAGACGTCGATGGTCTCGGCGGCCTGGGCCAGGCCGGTGACGGCCAGCAGGCTCAGGGTCAGCAGGTGTTTGCGCAGCATGGTGAAGCACTCCTTTAAAAACGCAGGCGGATACCGGCGACCAGGCGGGTCTGGC  >  NC_002947/2508920‑2509212

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