Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,096,465	C→T	9.4%	V101I (GTC→ATC) ‡	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,465	0	C	T	9.4%	144.6 / 3.8	53	V101I (GTC→ATC) ‡	argT	lysine/arginine/ornithine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (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

AGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGATATC  >  NC_002947/5096317‑5096612                                                                                                                                                      |                                                                                                                                                    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                                                                                                                       <  5:253283/149‑1 (MQ=255)                               gAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACAt                                                                                                                       >  3:79606/1‑149 (MQ=255)                               gAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACAt                                                                                                                       <  2:153471/149‑1 (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                                                                                                                                                 <  3:317263/64‑1 (MQ=255)                                                                                           ggtgCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTcgaccgac                                                                                                                                                 >  4:317263/1‑64 (MQ=255)                                                                                            gtgCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTc                                                           <  7:94479/149‑1 (MQ=255)                                                                                              gCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCt                                                         >  5:283653/1‑149 (MQ=255)                                                                                                    cTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTgc                                                                                             >  2:118546/1‑107 (MQ=255)                                                                                                    cTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTgc                                                                                             <  1:118546/107‑1 (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                                                             <  7:306602/91‑1 (MQ=255)                                                                                                                                                    gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAAc                                                             >  8:306602/1‑91 (MQ=255)                                                                                                                                                    gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGatat   >  3:106731/1‑149 (MQ=255)                                                                                                                                                    gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGata    <  2:333477/148‑1 (MQ=255)                                                                                                                                                    gaccgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGACCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGatat   >  4:328585/1‑149 (MQ=255)                                                                                                                                                     accgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGAcc                                                 >  4:299347/1‑102 (MQ=255)                                                                                                                                                     accgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGAcc                                                 <  3:299347/102‑1 (MQ=255)                                                                                                                                                     accgacTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGATATc  <  4:106731/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACG‑GTGGTGCCGTCCTTCATGACCAGGCGCGCCGGAGTCAGGTAGTAGCGCTTGGTGAAGTCGACCGACTTCTTGCGGTCATCGGTGATCGACATGGACGACAGGATGGCGTCGATCTTGCGCACTTTCAGCGCCGGGATCAGGCCGTCGAACTCCTGCTCGACCCAGGTACACTTGGTTTTCATCTCTTCGCACAGCGCGTTGCCGATATC  >  NC_002947/5096317‑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: ‑