Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,592,255	C→T	5.7%	D67D (GAC→GAT)	PP_4920 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,592,255	0	C	T	5.7%	210.2 / 2.3	70	D67D (GAC→GAT)	PP_4920	lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base T (2/2);  total (36/34)
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.42e-01

CCTGTTCACGCCCAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCAAGCCAAGGTACGCGA  >  NC_002947/5592117‑5592398                                                                                                                                           |                                                                                                                                                ccTGTTCACGCCCAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCAtcgtc                                                                                                                                        >  5:169013/1‑148 (MQ=255)     ttCACGCCCAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGaaa                                                                                                                                    <  5:178213/148‑1 (MQ=255)        aCGCCCAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGATGGCATCGTCGAAAAACg                                                                                                                                >  7:95965/1‑149 (MQ=255)        aCGCCCAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAAc                                                                                                                                 <  4:172830/148‑1 (MQ=255)         cGCCCAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGc                                                                                                                               <  7:199175/149‑1 (MQ=255)            ccAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCctg                                                                                                                            >  5:168905/1‑149 (MQ=255)                   ccGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAAc                                                                                                                     >  8:128951/1‑149 (MQ=255)                            gCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCg                                                                                                               >  8:219154/1‑146 (MQ=255)                                gAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTTAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACaa                                                                                                        <  6:97640/149‑1 (MQ=255)                                  ggTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGCCAATc                                                                                                      >  2:76805/1‑149 (MQ=255)                                  ggTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATc                                                                                                      <  6:159207/149‑1 (MQ=255)                                    tCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAg                                                                                                    <  4:254465/149‑1 (MQ=255)                                     cAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAg                                                                                                    >  4:112038/1‑148 (MQ=255)                                      aaGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCa                                                                                                  >  5:267056/1‑149 (MQ=255)                                       aGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCTACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCAt                                                                                                 <  3:340221/149‑1 (MQ=255)                                       aGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCATGCCCAAGCTCTACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCAt                                                                                                 <  5:240306/149‑1 (MQ=255)                                        gcgcGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATg                                                                                                >  7:68256/1‑149 (MQ=255)                                        gcgcGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATg                                                                                                <  6:168905/149‑1 (MQ=255)                                        gcgcGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATg                                                                                                <  7:219154/149‑1 (MQ=255)                                               gCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGgcaccg                                                                                          >  6:296285/1‑148 (MQ=255)                                               gCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGgcaccg                                                                                          >  2:32797/1‑148 (MQ=255)                                                  tGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGgcaccgcacc                                                                                      >  4:72940/1‑149 (MQ=255)                                                    ggAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAg                                                                                    <  7:264775/149‑1 (MQ=255)                                                    ggAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTCACCGAAGACAATCAGCATGGCACCGCACCa                                                                                     <  1:32797/148‑1 (MQ=255)                                                        cATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCAc                                                                                <  6:267056/149‑1 (MQ=255)                                                         aTGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCAcc                                                                               <  5:43272/149‑1 (MQ=255)                                                         aTGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCAca                                                                               >  6:43272/1‑148 (MQ=255)                                                            tCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCt                                                                             <  7:128951/148‑1 (MQ=255)                                                            tCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTg                                                                            >  6:36409/1‑149 (MQ=255)                                                            tCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTg                                                                            >  6:48332/1‑149 (MQ=255)                                                             cAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTg                                                                            >  6:164199/1‑148 (MQ=255)                                                              aaaCCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCa                                                                          <  5:296285/149‑1 (MQ=255)                                                                 cccGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCCGCCACCCTGCAGGc                                                                       >  4:352946/1‑149 (MQ=255)                                                                 cccGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCCCCCTGCAGGc                                                                       >  3:306149/1‑149 (MQ=255)                                                                  ccGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAg                                                                                                            <  5:141660/111‑1 (MQ=255)                                                                  ccGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAg                                                                                                            >  6:141660/1‑111 (MQ=255)                                                                   cGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCt                                                                     <  5:36409/149‑1 (MQ=255)                                                                         cAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGATGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAAc                                                               <  8:95965/149‑1 (MQ=255)                                                                                 gCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAatat                                                       <  5:164199/149‑1 (MQ=255)                                                                                        cccGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGc                                                  >  4:249959/1‑147 (MQ=255)                                                                                         ccGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAA‑GCCTGCTGCAACTGACCGAAGACAATAAGCATGGCAACGCACAAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGt                                              >  6:385966/1‑149 (MQ=255)                                                                                           gCTGGTGAACATCGACATGATCCACTTCCCGGTCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTg                                             <  2:57972/149‑1 (MQ=255)                                                                                               gTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGa                                         >  4:197298/1‑149 (MQ=255)                                                                                                         aCATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTaa                                >  2:188408/1‑148 (MQ=255)                                                                                                         aCATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTaa                                >  7:46136/1‑148 (MQ=255)                                                                                                         aCATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGACATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTg                                                                            >  1:248480/1‑104 (MQ=255)                                                                                                         aCATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGACATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTg                                                                            <  2:248480/104‑1 (MQ=255)                                                                                                            tGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAAcagc                            <  3:112038/149‑1 (MQ=255)                                                                                                            tGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAAcag                             >  7:248884/1‑148 (MQ=255)                                                                                                             gATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAAcag                             >  7:42438/1‑147 (MQ=255)                                                                                                              aTCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAAcagcag                          <  5:48332/149‑1 (MQ=255)                                                                                                                ccACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGcc                                           >  5:52432/1‑130 (MQ=255)                                                                                                                ccACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGcc                                           <  6:52432/130‑1 (MQ=255)                                                                                                                 cACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTa                       >  4:110614/1‑149 (MQ=255)                                                                                                                  aCTTCCCGGCCCTGCCCAAGCTCGATGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTAc                      >  6:123596/1‑149 (MQ=255)                                                                                                                   cTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTACCGACAAACGTAACAGCAGCTACt                     >  2:282476/1‑149 (MQ=255)                                                                                                                        cGGCCCTGCCCAAGCTCGATGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCaa                <  5:123596/149‑1 (MQ=255)                                                                                                                        cGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCaa                <  3:197298/149‑1 (MQ=255)                                                                                                                        cGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTACCGACAAACGTAACAGCAGCTACTTGCaa                <  1:282476/149‑1 (MQ=255)                                                                                                                         ggCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTTGCCAGTGc                                            >  4:274209/1‑120 (MQ=255)                                                                                                                         ggCCCTGCCCAAGATCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTTGCCAGTGc                                            <  3:274209/120‑1 (MQ=255)                                                                                                                            ccTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCAAGcc             >  8:140515/1‑148 (MQ=255)                                                                                                                               gCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACGCAAGACAATCAGCATGGCGCCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCAAGCCAAgg         >  1:263313/1‑149 (MQ=255)                                                                                                                                cccAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCAAGCCAAg          <  3:249959/147‑1 (MQ=255)                                                                                                                                cccAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCAAGCCAAg          <  3:72940/147‑1 (MQ=255)                                                                                                                                  cAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCt                                                                     <  7:172565/86‑1 (MQ=255)                                                                                                                                  cAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCt                                                                     >  8:172565/1‑86 (MQ=255)                                                                                                                                   aaGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAAcagc                            <  7:76372/126‑1 (MQ=255)                                                                                                                                   aaGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAAcagc                            >  8:76372/1‑126 (MQ=255)                                                                                                                                      cTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCGCCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCAAGCCAAGGTACGCGa  <  2:263313/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                CCTGTTCACGCCCAACTACCGGGCCCCGCTCGAGGTCAAGCGCGACGCCTGGGAGCATGTCAAACCCGGCTGCAGCGAGAGCGACTGCCCGCTGGTGAACATCGACATGATCCACTTCCCGGCCCTGCCCAAGCTCGACGGCATCGTCGAAAAACGCCTGCTGCAACTGACCGAAGACAATCAGCATGGCACCGCACCAGCCACCCTGCAGGCCTACGAACAGCAATATCTGGCCAGTGCCGACAAACGTAACAGCAGCTACTTGCAAGCCAAGGTACGCGA  >  NC_002947/5592117‑5592398

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