Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,429,250	G→A	6.6%	R194C (CGC→TGC)	PP_4768 ←	exonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,429,250	0	G	A	6.6%	226.9 / 4.6	76	R194C (CGC→TGC)	PP_4768	exonuclease
Reads supporting (aligned to +/- strand):  ref base G (42/29);  new base A (2/3);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

TCAGTTTAACGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCGGTGCGCGGTCCCC  >  NC_002947/5429110‑5429398                                                                                                                                             |                                                                                                                                                       tCAGTTTAACGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCg                                                                                                                                               >  8:155210/1‑149 (MQ=255)    gTTTAACGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAg                                                                                                                                            >  6:77981/1‑149 (MQ=255)    gTTTAACGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAg                                                                                                                                            >  3:77145/1‑149 (MQ=255)       tAACGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCg                                                                                                                                         >  5:99610/1‑149 (MQ=255)        aaCGGACTGTCTGGGTGAAGGGACGTTAAAGCACGGGTACAATTGAGCGGGACAATAGTCGTGCGTTAGAAGTAGGCGCCCCTCCACCCGACGTACACTGGATAATCCCTGCAGAACTAGGCCCGCGCTTCGCGACG‑CGTCGCAGCAgg                                                                                                                                        >  8:344524/1‑149 (MQ=255)         aCGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTcgcg                                                                                                                                                        >  8:82241/1‑133 (MQ=255)         aCGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTcgcg                                                                                                                                                        <  7:82241/133‑1 (MQ=255)           ggCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGt                                                                                                                                     <  7:197173/149‑1 (MQ=255)              cTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGtttt                                                                                                                                  >  8:244568/1‑149 (MQ=255)              cTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGtttt                                                                                                                                  >  1:318409/1‑149 (MQ=255)               tGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGATAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTc                                                                                                                                 >  1:352388/1‑149 (MQ=255)                gTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCt                                                                                                                                >  8:158938/1‑149 (MQ=255)                 tCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCt                                                                                                                                <  8:88776/148‑1 (MQ=255)                 acTGGGTCAAGGCACGTTACAGAACGGGTTCAATTGAGAGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCACGACGTTCACTGGCTAATCCCTGCAGTACATGGCCAGCCCTACGCGCAG‑CTTCGCATCCCGGGTTTTCTg                                                                                                                               >  2:156391/2‑149 (MQ=255)                   tGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTg                                                                                                                               >  6:164118/1‑147 (MQ=255)                   tGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTg                                                                                                                               >  6:322213/1‑147 (MQ=255)                         aaGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCgg                                                                                                                        <  3:191178/148‑1 (MQ=255)                           gggACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCaa                                                                                                                     <  5:77981/149‑1 (MQ=255)                            ggACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCa                                                                                                                      >  3:209018/1‑147 (MQ=255)                             gACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAag                                                                                                                   <  7:155210/149‑1 (MQ=255)                                gTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAagcag                                                                                                                >  5:229864/1‑149 (MQ=255)                                     aGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGt                                                                                                           <  5:322213/149‑1 (MQ=255)                                      gAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTg                                                                                                          >  7:280167/1‑149 (MQ=255)                                       aaCGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGTGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGc                                                                                                         >  7:15990/1‑149 (MQ=255)                                         cGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGa                                                                                                       >  2:173966/1‑149 (MQ=255)                                         cGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGa                                                                                                       >  2:168401/1‑149 (MQ=255)                                          gggTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGAc                                                                                                      <  5:164118/149‑1 (MQ=255)                                           ggTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACa                                                                                                     <  4:77145/149‑1 (MQ=255)                                                aaTTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTg                                                                                                 <  2:150837/148‑1 (MQ=255)                                                 aTTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCt                                                                                               >  7:184895/1‑149 (MQ=255)                                                 aTTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCt                                                                                               <  7:158938/149‑1 (MQ=255)                                                 aTTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCt                                                                                               >  5:32908/1‑149 (MQ=255)                                                  ttGAGCGGCAGAATTCTGGTGCGTTCGACGTTGGCGACCCTCCCCCCGACGTTCACTGGCTAATGCCTGCAGTACTTGGCCAGCGCTTCGGGGAG‑CTTCGCTGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTTCt                                                                                               <  1:156391/148‑1 (MQ=255)                                                       cGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAAccc                                                                                         >  2:91715/1‑149 (MQ=255)                                                          cacaATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGt                                                                                      <  2:318409/149‑1 (MQ=255)                                                           acaATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTg                                                                                     <  6:32908/149‑1 (MQ=255)                                                               ttCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTcc                                                                                 >  5:375874/1‑149 (MQ=255)                                                               ttCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTcc                                                                                 >  5:9921/1‑149 (MQ=255)                                                               ttCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTcc                                                                                 >  1:163218/1‑149 (MQ=255)                                                               ttCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTcc                                                                                 <  4:209018/149‑1 (MQ=255)                                                                  tGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTgg                                                                              <  7:244568/149‑1 (MQ=255)                                                                         ttCGACGTAGGCGCCCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCCGCTCGTGCGACAGTTGGTGAACCCCGTGGTCCTGGCGCAATa                                                                       >  6:60610/1‑149 (MQ=255)                                                                               gTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCg                                                                >  4:260537/1‑149 (MQ=255)                                                                                 aGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCACAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGc                                                               >  3:279382/1‑148 (MQ=255)                                                                                  ggCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCa                                                             <  1:168401/149‑1 (MQ=255)                                                                                  ggCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCa                                                             <  1:173966/149‑1 (MQ=255)                                                                                   gCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATaca‑c                                                                   >  6:296875/1‑142 (MQ=255)                                                                                   gCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATaca‑c                                                                   <  5:296875/142‑1 (MQ=255)                                                                                        ccTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACt                                                       >  1:2371/1‑149 (MQ=255)                                                                                        ccTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCACAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTg                                                                                     <  4:215239/120‑1 (MQ=255)                                                                                        ccTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCACAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTg                                                                                     >  3:215239/1‑120 (MQ=255)                                                                                        ccTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCACAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTg                                                                                     <  4:210462/120‑1 (MQ=255)                                                                                           ccccccGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCg                                                                                                                                         <  3:360491/65‑1 (MQ=255)                                                                                           ccccccGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCg                                                                                                                                         >  4:360491/1‑65 (MQ=255)                                                                                            cccccGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCACAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTgc                                                     <  4:279382/147‑1 (MQ=255)                                                                                                cGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAgg                                               >  4:309466/1‑149 (MQ=255)                                                                                                          tGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGc                                                               >  5:332800/1‑123 (MQ=255)                                                                                                          tGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGc                                                               <  6:332800/123‑1 (MQ=255)                                                                                                          tGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCa                                     >  7:155953/1‑149 (MQ=255)                                                                                                            gCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGgcgc                                                                           <  4:97109/110‑1 (MQ=255)                                                                                                            gCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGgcgc                                                                           >  3:97109/1‑110 (MQ=255)                                                                                                            gCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCAtt                                   >  4:173403/1‑149 (MQ=255)                                                                                                            gCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCAtt                                   >  3:51424/1‑149 (MQ=255)                                                                                                             cTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCAttt                                  <  6:99610/149‑1 (MQ=255)                                                                                                              tAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTc                                 >  5:283843/1‑149 (MQ=255)                                                                                                                 tCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCgg                               >  7:134014/1‑148 (MQ=255)                                                                                                                    cTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATc                           >  1:230489/1‑149 (MQ=255)                                                                                                                      gCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCgg                               >  3:23867/1‑143 (MQ=255)                                                                                                                      gCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCgg                               <  4:23867/143‑1 (MQ=255)                                                                                                                       cAGTACTTGGCCAGCGCTTCGCGCAGACTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATa                                                         >  1:191945/1‑117 (MQ=255)                                                                                                                       cAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATa                                                         <  2:191945/116‑1 (MQ=255)                                                                                                                         gTACTTGGCCAGTGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCa                      <  8:15990/149‑1 (MQ=255)                                                                                                                             ttGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCCGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGc                  >  5:106187/1‑149 (MQ=255)                                                                                                                               ggCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGc                  <  3:309466/147‑1 (MQ=255)                                                                                                                                gCCAGCGCTGCGCGCAG‑CTTCGCAGCCGGGGTTGTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGTCGGCCATTTCGGCAGCGGCCATGGCCCg               <  7:344524/149‑1 (MQ=255)                                                                                                                                 ccAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCgg              >  2:219292/1‑149 (MQ=255)                                                                                                                                  cAGCGCTTCGCGCAG‑CTTCGCTGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCGGGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCGGt             <  2:352388/149‑1 (MQ=255)                                                                                                                                  cAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACACCCCTCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCgg              <  1:91715/149‑1 (MQ=255)                                                                                                                                           gcgcAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCCGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCGGTGCGCGGTcc    <  6:106187/149‑1 (MQ=255)                                                                                                                                             gcAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCGGTGCGCGGTcccc  <  2:163218/149‑1 (MQ=255)                                                                                                                                             gcAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCGGTGCGCGGTcccc  >  7:182834/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                       TCAGTTTAACGGCCTGTCTGGGTGAAGGGACGTTACAGAACGGGTTCAATTGAGCGGCACAATTCTGGTGCGTTCGACGTAGGCGACCCTCCCCCCGACGTTCACTGGCTAATCCCTGCAGTACTTGGCCAGCGCTTCGCGCAG‑CTTCGCAGCCGGGGTTTTCTGCAGCCGGCAAAGCAGCTCGTGCGACAGTTGCTGAACCCCGTGGTCCTGGCGCAATACA‑CCCGCCAGATACTGCAGCAGGTTGGCGGCCATTTCGGCATCGGCCATGGCCCGGTGCGCGGTCCCC  >  NC_002947/5429110‑5429398

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