Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,931	C→G	5.9%	A156A (GCC→GCG)	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,931	0	C	G	5.9%	271.0 / 7.0	85	A156A (GCC→GCG)	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.54e-01

ATCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGA  >  NC_002947/5378782‑5379080                                                                                                                                                       |                                                                                                                                                         aTCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGcc                                                                                                                                                          >  1:1484920/1‑150 (MQ=255) aTCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACAGCCAAAGA‑ACGTCGTGTGGTTGACGAGTATCCCGGCGCCGTCGACCTGGcc                                                                                                                                                          >  2:1468654/1‑150 (MQ=255)  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                                                                                                                                                 <  1:2408417/95‑2 (MQ=255)                                                              tcgacaGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtc                                                                                                                                                 >  2:2408417/4‑97 (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                                                                                  <  1:538468/150‑1 (MQ=255)                                                                      cagcACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCg                                                                                  <  2:28083/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:333354/1‑150 (MQ=255)                                                                                 ccAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGagca                                                                       >  1:842954/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                                                                                 <  2:1423246/133‑1 (MQ=255)                                                                                        cTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCgg                                                                                 >  1:1423246/1‑133 (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        <  1:765712/150‑1 (MQ=255)                                                                                                                                                 ccTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGa        <  2:1694959/150‑1 (MQ=255)                                                                                                                                                 ccTGCCCCGCCAGGTACAGCTCGCACCCAAGCGCTC‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACACGGGCGACGAGCAGCGCGCCCGCCACCTGCTGCTCAGTTACCAGCTGTAAGAACTGGACAACCTCGGCCTGGGCGa        <  2:522895/150‑1 (MQ=255)                                                                                                                                                  cTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAg       <  1:2327423/150‑1 (MQ=255)                                                                                                                                                  cTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAg       >  2:439462/1‑150 (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  >  2:645821/1‑150 (MQ=255)                                                                                                                                                       ccGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGa  >  1:707066/1‑150 (MQ=255)                                                                                                                                                       ccGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGa  >  1:1212115/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                         ATCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCA‑CCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTG‑GA‑AT‑CAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCACGA  >  NC_002947/5378782‑5379080

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