Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,096,464	A→G	9.4%	V101A (GTC→GCC) ‡	argT ←	lysine/arginine/ornithine ABC transporter substrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,096,464	0	A	G	9.4%	142.0 / 5.4	53	V101A (GTC→GCC) ‡	argT	lysine/arginine/ornithine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base G (2/3);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGATATC  >  NC_002947/5096316‑5096612                                                                                                                                                      |                                                                                                                                                     cAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcga                                                                                                                                                      <  8:130830/149‑1 (MQ=255) cAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcga                                                                                                                                                      <  7:201744/149‑1 (MQ=255)  aGCACTTCCTTGGCGAAGCGGTCGTGAATGTAGCCGCGCTGCAAGCAGATCTTCTTGCCTTTCATTGCTTCCAGGCTGTCGCTGACT‑GTGGTGCCGTCCTTCATGACCAGTCTCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcgac                                                                                                                                                     <  4:191812/149‑1 (MQ=255)  aGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACTCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcgac                                                                                                                                                     <  5:57266/149‑1 (MQ=255)   gCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcgacc                                                                                                                                                    <  3:282646/149‑1 (MQ=255)       ttCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGAct                                                                                                                                                >  3:135394/1‑149 (MQ=255)         ccTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGGTCGActt                                                                                                                                               >  2:202603/1‑148 (MQ=255)           ttGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGActtctt                                                                                                                                            >  1:153471/1‑149 (MQ=255)             ggCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTg                                                                                                                                           <  4:128838/148‑1 (MQ=255)               ccaaGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCgg                                                                                                                                        <  5:323317/147‑1 (MQ=255)                gAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGGGCTTGGTGAAGTCGGTCGActtcaccaagc                                                                                                                                       <  1:31123/149‑8 (MQ=255)                gAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGGTCGActtcaccaagc                                                                                                                                       <  1:202603/149‑8 (MQ=255)                     ggTCGTGATTGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCg                                                                                                                                  <  6:100699/149‑1 (MQ=255)                       tCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGt                                                                                                                                >  6:20604/1‑149 (MQ=255)                            aaTGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATc                                                                                                                            <  1:322045/148‑1 (MQ=255)                              tGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGTGTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGTTGAATACGACCGACTTCTTGCGGTCATCGGTGATCg                                                                                                                           >  1:290476/1‑148 (MQ=255)                                gAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACAt                                                                                                                       <  2:153471/149‑1 (MQ=255)                                gAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACAt                                                                                                                       <  5:253283/149‑1 (MQ=255)                                gAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACAt                                                                                                                       >  3:79606/1‑149 (MQ=255)                                          gCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAg                                                                                                              >  3:296682/1‑148 (MQ=255)                                                       cttGCCTTTCAGATCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCt                                                                                                >  6:25261/1‑149 (MQ=255)                                                        ttGCCAATCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCtt                                                                                               >  8:94479/1‑149 (MQ=255)                                                          gCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTgc                                                                                             >  5:157698/1‑149 (MQ=255)                                                               tCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCAc                                                                                         <  6:157698/148‑1 (MQ=255)                                                                       tCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATg                                                                                                          >  3:28763/1‑123 (MQ=255)                                                                       tCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATg                                                                                                          <  4:28763/123‑1 (MQ=255)                                                                       tCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAgcgc                                                                                >  6:284804/1‑149 (MQ=255)                                                                       tCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAgcgc                                                                                >  7:335746/1‑149 (MQ=255)                                                                                  gCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAgg                                                                      <  4:296682/148‑1 (MQ=255)                                                                                    tGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCg                                                                   >  7:167351/1‑149 (MQ=255)                                                                                          gtggtgCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGaa                                                              <  4:135394/149‑1 (MQ=255)                                                                                            ggtgCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcgaccgac                                                                                                                                                 >  4:317263/1‑64 (MQ=255)                                                                                            ggtgCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcgaccgac                                                                                                                                                 <  3:317263/64‑1 (MQ=255)                                                                                             gtgCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTc                                                           <  7:94479/149‑1 (MQ=255)                                                                                               gCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCt                                                         >  5:283653/1‑149 (MQ=255)                                                                                                     cTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTgc                                                                                             <  1:118546/107‑1 (MQ=255)                                                                                                     cTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTgc                                                                                             >  2:118546/1‑107 (MQ=255)                                                                                                     cTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGa                                                   >  2:325850/1‑149 (MQ=255)                                                                                                          tGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAg                                              <  8:335746/149‑1 (MQ=255)                                                                                                              cAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATg                                                                                                          >  1:134437/1‑85 (MQ=255)                                                                                                              cAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATg                                                                                                          <  2:134437/85‑1 (MQ=255)                                                                                                                   gcgcCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTg                                     >  7:349020/1‑149 (MQ=255)                                                                                                                     gcCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGt                                   <  2:290476/149‑1 (MQ=255)                                                                                                                      ccGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGtt                                  >  2:228949/1‑149 (MQ=255)                                                                                                                         gAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTc                               <  1:228949/149‑1 (MQ=255)                                                                                                                              aGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGc                                                       >  5:162326/1‑120 (MQ=255)                                                                                                                              aGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGc                                                       <  6:162326/120‑1 (MQ=255)                                                                                                                                gtagtaGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTc                               <  7:215448/142‑1 (MQ=255)                                                                                                                                gtagtaGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGACGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTc                               >  8:215448/1‑142 (MQ=255)                                                                                                                                 tagtagCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTc                       <  5:284804/149‑1 (MQ=255)                                                                                                                                   gtagCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGACCTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCg                      >  4:189852/1‑148 (MQ=255)                                                                                                                                      gcgcTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGcaca                  <  5:20604/149‑1 (MQ=255)                                                                                                                                         cTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAgcg               >  6:347348/1‑149 (MQ=255)                                                                                                                                          ttGGTGAAGTCGACCGACTTCTAACGGTCATCGGTGATCGACATGGTCGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTAATCTCTTCGCACAgcgc              >  7:232102/1‑149 (MQ=255)                                                                                                                                            ggTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAgcgc              <  4:79606/147‑1 (MQ=255)                                                                                                                                            caagAAGTCGGTCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTc                               <  3:131164/127‑1 (MQ=255)                                                                                                                                            caagAAGTCGGTCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTc                               >  4:131164/4‑130 (MQ=255)                                                                                                                                             gTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGtt           <  5:25261/149‑1 (MQ=255)                                                                                                                                                  gTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTCCACGTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGa      >  4:98479/1‑149 (MQ=255)                                                                                                                                                  gTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGa      >  8:10070/1‑149 (MQ=255)                                                                                                                                                     gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAAc                                                             >  8:306602/1‑91 (MQ=255)                                                                                                                                                     gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAAc                                                             <  7:306602/91‑1 (MQ=255)                                                                                                                                                     gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGatat   >  3:106731/1‑149 (MQ=255)                                                                                                                                                     gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGata    <  2:333477/148‑1 (MQ=255)                                                                                                                                                     gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGACCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGatat   >  4:328585/1‑149 (MQ=255)                                                                                                                                                      accgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGAcc                                                 <  3:299347/102‑1 (MQ=255)                                                                                                                                                      accgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGAcc                                                 >  4:299347/1‑102 (MQ=255)                                                                                                                                                      accgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGATATc  <  4:106731/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGATATC  >  NC_002947/5096316‑5096612

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