Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,863,240	T→A	8.1%	L309Q (CTG→CAG)	PP_1666 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,863,240	0	T	A	8.1%	177.8 / 5.8	62	L309Q (CTG→CAG)	PP_1666	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base A (3/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

TGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGA  >  NC_002947/1863093‑1863385                                                                                                                                                    |                                                                                                                                                  tGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctg                                                                                                                                                  <  6:70945/149‑1 (MQ=255)      gTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctggcc                                                                                                                                               >  1:79575/1‑147 (MQ=255)         cAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCt                                                                                                                                          >  2:39704/1‑149 (MQ=255)             gCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAg                                                                                                                                      >  3:286163/1‑149 (MQ=255)               aTGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAgg                                                                                                                                     >  6:317404/1‑148 (MQ=255)                 gAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAg                                                                                                                                  >  6:263440/1‑149 (MQ=255)                        cAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGt                                                                                                                           <  1:39704/149‑1 (MQ=255)                        cAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGt                                                                                                                           >  8:40798/1‑149 (MQ=255)                         aGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTc                                                                                                                          <  5:317404/149‑1 (MQ=255)                            gcTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAg                                                                                                                       <  8:220868/149‑1 (MQ=255)                                tGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCa                                                                                                                     <  7:144475/147‑1 (MQ=255)                                          gCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGTACCGGCGGt                                                                                                         >  7:52242/1‑149 (MQ=255)                                               cgTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCAc                                                                                                                <  1:135313/137‑1 (MQ=255)                                                 tGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGcc                                                                                                  <  3:130151/149‑1 (MQ=255)                                                   cTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCGGGCACCGGCGGTGCCGGCCAt                                                                                                <  8:127588/149‑1 (MQ=255)                                                   cTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCa                                                                                                 <  5:263440/148‑1 (MQ=255)                                                   cTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCAt                                                                                                <  2:3605/149‑1 (MQ=255)                                                    ttGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATg                                                                                               >  2:55205/1‑149 (MQ=255)                                                    ttGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATg                                                                                               >  4:214135/1‑149 (MQ=255)                                                      gAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGc                                                                                              >  6:23539/1‑148 (MQ=255)                                                      gAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGc                                                                                              <  8:148769/148‑1 (MQ=255)                                                           cTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAc                                                                                        <  5:23539/149‑1 (MQ=255)                                                                  ccACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGt                                                                                 <  1:55205/149‑1 (MQ=255)                                                                   cACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCgg                                                                                     >  2:273784/1‑144 (MQ=255)                                                                   cACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCgg                                                                                     <  1:273784/144‑1 (MQ=255)                                                                          gCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGa                                                                              >  6:226203/1‑144 (MQ=255)                                                                          gCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGa                                                                              <  5:226203/144‑1 (MQ=255)                                                                           cAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGTACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGAcgccg                                                                         <  8:52242/148‑1 (MQ=255)                                                                            aGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGAcgccg                                                                         >  4:325751/1‑147 (MQ=255)                                                                            aGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCa                                                                       <  3:325751/149‑1 (MQ=255)                                                                             ggcggcTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAg                                                                      >  8:111345/1‑149 (MQ=255)                                                                              gcggcTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGc                                                                     >  4:58498/1‑149 (MQ=255)                                                                                 gcTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCa                                                                                                 >  2:22854/1‑118 (MQ=255)                                                                                 gcTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCcgc                                                                  >  7:296850/1‑149 (MQ=255)                                                                                      aGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGc                                                             <  8:296850/149‑1 (MQ=255)                                                                                        gCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGccc                                                           >  3:61569/1‑149 (MQ=255)                                                                                          aCTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCcgc                                                         >  7:87418/1‑149 (MQ=255)                                                                                               ccTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTg                                                    >  4:30465/1‑149 (MQ=255)                                                                                                 tGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTgcg                                                  >  3:57098/1‑149 (MQ=255)                                                                                                  gACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTgcg                                                  >  6:208688/1‑148 (MQ=255)                                                                                                     cTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGAcgccgc                                                                        >  4:274662/1‑123 (MQ=255)                                                                                                     cTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGAcgccgc                                                                        >  4:239827/1‑123 (MQ=255)                                                                                                     cTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagc                                              <  3:214135/149‑1 (MQ=255)                                                                                                       aCGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagc                                              >  7:295492/1‑147 (MQ=255)                                                                                                                cGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCaa                                    <  4:286163/148‑1 (MQ=255)                                                                                                                cGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAg                                   <  8:87418/149‑1 (MQ=255)                                                                                                                    cAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGt                               >  7:279271/1‑149 (MQ=255)                                                                                                                     aaTCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGtt                              <  7:111345/149‑1 (MQ=255)                                                                                                                              cAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCtt                     <  1:2194/149‑1 (MQ=255)                                                                                                                              cAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCtt                     >  8:139766/1‑149 (MQ=255)                                                                                                                               aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTCGAACGGCTTg                    <  3:58498/149‑1 (MQ=255)                                                                                                                                gCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTCGAACGGCTTgg                   <  7:40798/149‑2 (MQ=255)                                                                                                                                 cTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCg                  >  3:279148/1‑149 (MQ=255)                                                                                                                                    cGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGca                                                <  5:279026/116‑1 (MQ=255)                                                                                                                                    cGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGca                                                >  6:279026/1‑116 (MQ=255)                                                                                                                                    cGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCATCACCCAAGCCGTCGAACGGCTTGGGt                 <  6:93099/147‑1 (MQ=255)                                                                                                                                    cGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGttt               <  4:61569/149‑1 (MQ=255)                                                                                                                                       ctgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCcgc                                                         >  6:260392/4‑104 (MQ=255)                                                                                                                                       ctgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCcgc                                                         <  5:260392/101‑1 (MQ=255)                                                                                                                                        tgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAg                                   <  4:149929/123‑1 (MQ=255)                                                                                                                                        tgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAg                                   >  3:149929/3‑125 (MQ=255)                                                                                                                                        tgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGc           >  7:232355/3‑149 (MQ=255)                                                                                                                                         ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt          <  5:208688/149‑1 (MQ=255)                                                                                                                                           aGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTgg        >  4:29055/1‑149 (MQ=255)                                                                                                                                                 gcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGa  <  2:163895/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGA  >  NC_002947/1863093‑1863385

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