Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,545,948	G→A	7.0%	R350H (CGC→CAC)	PP_2238 →	metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,545,948	0	G	A	7.0%	253.0 / 5.7	85	R350H (CGC→CAC)	PP_2238	metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base A (3/3);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

CCTGCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACTGAGC  >  NC_002947/2545800‑2546096                                                                                                                                                     |                                                                                                                                                     ccTGCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCg                                                                                                                                                      >  8:7062/1‑149 (MQ=255)        gggCGCGGCTGCTGGGGGACCCGACCCGGGTCACCTGCGGCCGGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGCGAATATGCACCAGGTc                                                                                                                                               <  8:242219/149‑1 (MQ=255)            gcgGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCAGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGTATATTCGCCAGGTCATg                                                                                                                                            >  1:258232/1‑148 (MQ=255)              ggCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAg                                                                                                                                         >  8:39771/1‑149 (MQ=255)               gctgctGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGc                                                                                                                                        <  6:178563/149‑1 (MQ=255)                   ctgGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAg                                                                                                                                      <  8:101516/147‑1 (MQ=255)                    tggtggACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGAc                                                                                                                                   >  6:327801/1‑149 (MQ=255)                     ggtggACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACg                                                                                                                                  >  6:70196/1‑149 (MQ=255)                           cccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAg                                                                                                                            >  8:157739/1‑149 (MQ=255)                           cccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACAACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAg                                                                                                                            >  8:30391/1‑149 (MQ=255)                            ccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGc                                                                                                                           >  3:51516/1‑149 (MQ=255)                            ccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGc                                                                                                                           >  6:26642/1‑149 (MQ=255)                             cgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGcc                                                                                                                          >  1:78662/1‑149 (MQ=255)                              gacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCt                                                                                                                         >  8:40568/1‑149 (MQ=255)                               acccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCtt                                                                                                                        <  7:30391/149‑1 (MQ=255)                               acccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCtt                                                                                                                        >  6:314759/1‑149 (MQ=255)                                   gggTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGc                                                                                                                    <  5:314759/149‑1 (MQ=255)                                     gTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGc                                                                                                                    >  6:343602/1‑147 (MQ=255)                                          cTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAACCCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGtttt                                                                                                             >  5:180640/1‑149 (MQ=255)                                             cGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTc                                                                                                           <  5:343602/148‑1 (MQ=255)                                                  tgctgGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGt                                                                                                                <  2:93650/138‑1 (MQ=255)                                                  tgctgGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGt                                                                                                                >  1:93650/1‑138 (MQ=255)                                                    ctgGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTgg                                                                                                     <  7:157739/147‑1 (MQ=255)                                                      ggACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTc                                                                                                  >  7:240207/1‑148 (MQ=255)                                                           aCCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGc                                                                                             >  6:182830/1‑148 (MQ=255)                                                               ggcggcAGAGGTGGTTCTTGTCGAGGTGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGaa                                                                                          >  8:165016/1‑147 (MQ=255)                                                               ggcggcAGAGGTGGTTCTTGTCGAGGTGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGaa                                                                                          <  7:165016/147‑1 (MQ=255)                                                                 cggcAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGCCTCGAAccag                                                                                      <  4:312369/149‑3 (MQ=255)                                                                 cggcAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGCCTCGAAccag                                                                                      <  5:26642/149‑3 (MQ=255)                                                                       agGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCgg                                                                                <  5:70196/149‑1 (MQ=255)                                                                        ggtggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCggg                                                                               <  5:327801/149‑1 (MQ=255)                                                                        ggtggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCgg                                                                                <  5:182830/148‑1 (MQ=255)                                                                          tggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGa                                                                              <  1:232469/148‑1 (MQ=255)                                                                          tggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAg                                                                             <  7:40568/149‑1 (MQ=255)                                                                           ggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAgg                                                                            >  7:237831/1‑149 (MQ=255)                                                                            gtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGTGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGt                                                                           >  6:186746/1‑149 (MQ=255)                                                                               cTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtcat                                                                        >  2:187145/1‑149 (MQ=255)                                                                                 tGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtcatca                                                                      >  5:251747/1‑149 (MQ=255)                                                                                  gTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGTGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtcatc                                                                       <  5:186746/147‑1 (MQ=255)                                                                                  gTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAc                                                                     >  2:77655/1‑149 (MQ=255)                                                                                   tCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAcc                                                                    <  1:187145/149‑1 (MQ=255)                                                                                   tCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAcc                                                                    >  1:260944/1‑149 (MQ=255)                                                                                   tCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCATGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAcc                                                                    >  8:262902/1‑149 (MQ=255)                                                                                      aGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAg                                                                 >  2:180347/1‑149 (MQ=255)                                                                                       ggggCTAAACACCACCACACTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGCAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCTAGGCGAACCTCGGCCGGGAGGTCATCACCGAGc                                                                >  3:193189/1‑149 (MQ=255)                                                                                            tAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACg                                                           >  2:26476/1‑149 (MQ=255)                                                                                               aCCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGtt                                                        >  1:303901/1‑149 (MQ=255)                                                                                               aCCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGtt                                                        <  2:303901/149‑1 (MQ=255)                                                                                                  ccaccaCGCTTAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGAc                                                     <  4:51516/149‑1 (MQ=255)                                                                                                  ccaccaCGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGAc                                                     <  6:251747/149‑1 (MQ=255)                                                                                                  ccaccaCGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGAc                                                     <  1:77655/149‑1 (MQ=255)                                                                                                  ccaccaCGCTGAGCAAGTCCTACAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAACGAGCTGACGGTTGAc                                                     >  8:220694/1‑149 (MQ=255)                                                                                                    accacGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGa                                                   <  7:262902/149‑1 (MQ=255)                                                                                                     ccacGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAg                                                  >  2:305434/1‑149 (MQ=255)                                                                                                        cGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTCGAACCAGGCGaaa                                                                                         >  5:326990/1‑106 (MQ=255)                                                                                                        cGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTCGAACCAGGCGaaa                                                                                         <  6:326990/107‑2 (MQ=255)                                                                                                         gCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAAt                                              >  2:116059/1‑149 (MQ=255)                                                                                                         gCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAAt                                              <  7:7062/149‑1 (MQ=255)                                                                                                          cTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAAtt                                             <  7:39771/149‑1 (MQ=255)                                                                                                            gAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtc                                                                          >  5:82982/1‑118 (MQ=255)                                                                                                            gAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtc                                                                          <  6:82982/118‑1 (MQ=255)                                                                                                             aGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGa                                           <  8:240207/148‑1 (MQ=255)                                                                                                             aGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGa                                           <  1:180347/148‑1 (MQ=255)                                                                                                             aGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGa                                           <  2:78662/148‑1 (MQ=255)                                                                                                               cAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTg                                        <  1:305434/149‑1 (MQ=255)                                                                                                                         aaGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                              <  1:26476/149‑1 (MQ=255)                                                                                                                                 cgacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                      <  2:260944/149‑1 (MQ=255)                                                                                                                                 ccatgacCTGGCGAATATGCACCAGGTCTTGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                      >  1:192324/5‑149 (MQ=255)                                                                                                                                  gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                     <  7:220694/149‑1 (MQ=255)                                                                                                                                  gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                     >  6:58505/1‑149 (MQ=255)                                                                                                                                  gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                     >  2:300469/1‑149 (MQ=255)                                                                                                                                  gacgacCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGACCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGcc                     >  3:336193/1‑149 (MQ=255)                                                                                                                                        cTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAAttt               >  7:303373/1‑149 (MQ=255)                                                                                                                                         tGGAGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGGTGACGGTTGACGAGCAATTGAGTGCCGCTCGTGCTCGACGGCCTAAttac              >  6:989/1‑147 (MQ=255)                                                                                                                                          ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                              <  7:91336/127‑1 (MQ=255)                                                                                                                                          ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                              <  7:86608/127‑1 (MQ=255)                                                                                                                                          ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                              >  8:86608/6‑132 (MQ=255)                                                                                                                                          ggcgaatatGCACCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGc                              >  8:91336/6‑132 (MQ=255)                                                                                                                                          ggTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCg             >  5:44466/1‑149 (MQ=255)                                                                                                                                           gTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTc              <  8:303373/147‑1 (MQ=255)                                                                                                                                           gTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCgt            <  6:44466/149‑1 (MQ=255)                                                                                                                                              cATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                      <  2:24306/136‑1 (MQ=255)                                                                                                                                              cATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                      >  1:24306/1‑136 (MQ=255)                                                                                                                                                tatTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTc         >  5:55883/1‑147 (MQ=255)                                                                                                                                                tatTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTc         <  6:55883/147‑1 (MQ=255)                                                                                                                                                tatTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCAc       >  6:228805/1‑149 (MQ=255)                                                                                                                                                     gCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACTGAGc  <  1:300469/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCTGCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGCAGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACTGAGC  >  NC_002947/2545800‑2546096

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