Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,932	C→G	5.9%	R157G (CGC→GGC)	recN →	DNA repair/recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,932	0	C	G	5.9%	271.0 / 7.0	85	R157G (CGC→GGC)	recN	DNA repair/recombination protein
Reads supporting (aligned to +/- strand):  ref base C (41/39);  new base G (3/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

TCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGAA  >  NC_002947/5378783‑5379081                                                                                                                                                       |                                                                                                                                                         tCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGccc                                                                                                                                                          >  1:2157352/1‑150 (MQ=255)     cGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCa                                                                                                                                                      <  2:1202638/150‑1 (MQ=255)       ncACGCCCTGCCCGCTCGGCGACCTCAANGCGCTGGGCGAACTGCTGNTCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAgg                                                                                                                                                    >  1:391030/2‑150 (MQ=255)         acacccTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTa                                                                                                                                                  >  1:2408757/4‑150 (MQ=255)         aCGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACAGCCG‑GCGTAGCCTTCCTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTa                                                                                                                                                  <  1:1624805/150‑1 (MQ=255)           gCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACa                                                                                                                                                >  1:1991935/1‑150 (MQ=255)            cccTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAg                                                                                                                                               <  2:330152/150‑1 (MQ=255)            cccTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCACTGCTCGA‑CCAGGCGGACATCGGACTCCCGCTTGACGGCTAGGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAg                                                                                                                                               <  1:1854294/150‑1 (MQ=255)             ccTGCCCGCTCGTCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGc                                                                                                                                              >  2:538468/1‑150 (MQ=255)               tGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCtt                                                                                                                                            >  2:857808/1‑150 (MQ=255)                  cgtcTCGGCCACCTAAAGGTGCTGGGCGAATTGCTCATCCATCTCCACAGCCAGGACGAGCACCAGTCCCTGCTCACGACACTCACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGCCCTGGCCCGCCAGGTAGAGCTTGCa                                                                                                                                         <  2:2309498/147‑1 (MQ=255)                       cGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCaa                                                                                                                                    <  1:207117/150‑1 (MQ=255)                       cGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCaa                                                                                                                                    <  2:1162119/150‑1 (MQ=255)                            aCCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCt                                                                                                                               <  2:2157352/150‑1 (MQ=255)                            aCCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCt                                                                                                                               <  1:1000783/150‑1 (MQ=255)                              cTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTg‑g                                                                                                                            >  1:2376761/1‑150 (MQ=255)                               tCAAGGCGCTGGGCGAACTGCTGATCGATACCCACACCCAGCACGAGCACCAGTCCCTCCTCAAGACCGACCCCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑Ga                                                                                                                           <  1:1186761/150‑1 (MQ=255)                                 aaGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑At                                                                                                                        >  2:1910525/1‑150 (MQ=255)                                    gcgcTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAg                                                                                                                    <  1:241066/150‑1 (MQ=255)                                      gcTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCTCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGAc                                                                                                                  <  1:1468654/150‑1 (MQ=255)                                        ggggCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGAccc                                                                                                                <  1:2230238/149‑1 (MQ=255)                                          ggCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑c                                                                                                                      >  1:917135/1‑142 (MQ=255)                                          ggCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑c                                                                                                                      <  2:917135/142‑1 (MQ=255)                                                 tgctgATAGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCctggagc                                                                                                       <  1:1910525/150‑1 (MQ=255)                                                  gctgATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCctggagct                                                                                                      <  1:339001/150‑1 (MQ=255)                                                    tgATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCctggagctgg                                                                                                    >  1:1391340/1‑150 (MQ=255)                                                    tgATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCctggagctgg                                                                                                    >  1:424014/1‑150 (MQ=255)                                                      aTCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCctggagctggag                                                                                                  <  2:1604842/150‑1 (MQ=255)                                                      aTCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCctggagctggag                                                                                                  >  2:300118/1‑150 (MQ=255)                                                             tcgacaGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtc                                                                                                                                                  >  2:2408417/4‑97 (MQ=255)                                                             tcgacaGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtc                                                                                                                                                  <  1:2408417/95‑2 (MQ=255)                                                              ccataggcagcACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCtctc                                                                                          >  2:2327423/8‑150 (MQ=255)                                                              ccACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCtctc                                                                                          <  1:1809513/150‑1 (MQ=255)                                                              ccACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCtctc                                                                                          >  2:1809513/1‑150 (MQ=255)                                                               cacaGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTcc                                                                                         >  1:28083/1‑150 (MQ=255)                                                                    ccagcACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTc                                                                                    >  1:910925/1‑150 (MQ=255)                                                                     cagcACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCg                                                                                   <  2:28083/150‑1 (MQ=255)                                                                     cagcACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCg                                                                                   <  1:538468/150‑1 (MQ=255)                                                                       gcACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCggg                                                                                 <  2:1383920/150‑1 (MQ=255)                                                                       gcACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCggg                                                                                 >  1:1809631/1‑150 (MQ=255)                                                                             gCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCaa                                                                                       >  2:611339/1‑138 (MQ=255)                                                                             gCACCAGTCCCTGCCCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCaa                                                                                       <  1:611339/138‑1 (MQ=255)                                                                                ccAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGagca                                                                        >  1:842954/1‑150 (MQ=255)                                                                                ccAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGagca                                                                        >  1:333354/1‑150 (MQ=255)                                                                                      ccTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCGAGCGCTG‑GA‑AT‑CAGa                                                                                                                   >  2:1626994/1‑101 (MQ=255)                                                                                      ccTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCGAGCGCTG‑GA‑AT‑CAGa                                                                                                                   <  1:1626994/101‑1 (MQ=255)                                                                                       cTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCgg                                                                                  >  1:1423246/1‑133 (MQ=255)                                                                                       cTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCgg                                                                                  <  2:1423246/133‑1 (MQ=255)                                                                                       cTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGcgcc                                                                 >  1:748962/1‑150 (MQ=255)                                                                                       cTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGcgcc                                                                 >  1:1401669/1‑150 (MQ=255)                                                                                           tCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGcgcccgcc                                                             >  1:1947563/1‑150 (MQ=255)                                                                                                  cGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGc                                                                                <  1:1599294/124‑1 (MQ=255)                                                                                                  cGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGc                                                                                >  2:1599294/1‑124 (MQ=255)                                                                                                    acacCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAgctgc                                                    <  2:1484920/150‑1 (MQ=255)                                                                                                    acacCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAgctgc                                                    <  2:1809631/150‑1 (MQ=255)                                                                                                      accca‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTc                                                  >  2:551169/1‑150 (MQ=255)                                                                                                       ccca‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCa                                                 <  2:1368068/150‑1 (MQ=255)                                                                                                        cca‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAg                                                >  1:995152/1‑150 (MQ=255)                                                                                                          a‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCctct                                                                                           >  2:746000/1‑105 (MQ=255)                                                                                                          a‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCctct                                                                                           <  1:746000/105‑1 (MQ=255)                                                                                                          a‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGcgacga                                                                           <  1:1895696/121‑1 (MQ=255)                                                                                                          a‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGcgacga                                                                           >  2:1895696/1‑121 (MQ=255)                                                                                                          a‑ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGtt                                              <  2:1969511/150‑1 (MQ=255)                                                                                                               tCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCa                                          >  2:1921568/1‑150 (MQ=255)                                                                                                                   cTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACCAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTg                                      <  2:1391340/150‑1 (MQ=255)                                                                                                                    tGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTgg                                     <  2:2408757/150‑1 (MQ=255)                                                                                                                       ttGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGgaag                                  >  2:1313465/1‑150 (MQ=255)                                                                                                                        tGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGgaaga                                 >  2:765712/1‑150 (MQ=255)                                                                                                                          aCGAGTATGCCGGCGCCGTCGAACTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAAc                               <  2:910925/150‑1 (MQ=255)                                                                                                                           cgcgtctGCCGGCCCCGGCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGGTGCGATATCCAGGTACG‑‑TGGTG‑GGGCGCGCCTCTCGACCTCGGGCGACGCACGGTCCGGCTGCGGCCAGCTGCTCAGTTACCAGCTGGAAGAACt                              <  1:551169/144‑1 (MQ=255)                                                                                                                            gAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTg                             >  2:785535/1‑150 (MQ=255)                                                                                                                              gTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGa                           >  1:2258475/1‑150 (MQ=255)                                                                                                                                aTGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACa                         >  1:1874873/1‑150 (MQ=255)                                                                                                                                     ggCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTc                    >  1:1651813/1‑150 (MQ=255)                                                                                                                                           gTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTg              >  1:1403898/1‑150 (MQ=255)                                                                                                                                            tgtaCCTGGCGGGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGagc                                                                         <  1:588629/87‑1 (MQ=255)                                                                                                                                            tgtaCCTGGCGGGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGagc                                                                         >  2:588629/4‑90 (MQ=255)                                                                                                                                              taCCTGGCGGGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGc           >  1:1437812/2‑150 (MQ=255)                                                                                                                                               aCCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCg          >  1:853219/1‑150 (MQ=255)                                                                                                                                                ccTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGa         <  2:1694959/150‑1 (MQ=255)                                                                                                                                                ccTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGa         <  1:765712/150‑1 (MQ=255)                                                                                                                                                ccTGCCCCGCCAGGTACAGCTCGCACCCAAGCGCTC‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACACGGGCGACGAGCAGCGCGCCCGCCACCTGCTGCTCAGTTACCAGCTGTAAGAACTGGACAACCTCGGCCTGGGCGa         <  2:522895/150‑1 (MQ=255)                                                                                                                                                 cTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAg        >  2:439462/1‑150 (MQ=255)                                                                                                                                                 cTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAg        <  1:2327423/150‑1 (MQ=255)                                                                                                                                                  tGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGc       >  2:173952/1‑150 (MQ=255)                                                                                                                                                   ggCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTAGACAACCTCGGCCTGGGCGAGCa      >  1:408733/1‑150 (MQ=255)                                                                                                                                                      ccGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGa   >  1:1212115/1‑150 (MQ=255)                                                                                                                                                      ccGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGa   >  2:645821/1‑150 (MQ=255)                                                                                                                                                      ccGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGa   >  1:707066/1‑150 (MQ=255)                                                                                                                                                       cGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGaa  <  1:173952/150‑1 (MQ=255)                                                                                                                                                       cGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGaa  >  2:1580621/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                         TCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGAA  >  NC_002947/5378783‑5379081

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