Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,325,511	T→A	5.9%	noncoding (133/149 nt)	PP_mr52 →	PrrF ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,325,511	0	T	A	5.9%	209.4 / 2.5	68	noncoding (133/149 nt)	PP_mr52	PrrF ncRNA
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base A (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCT  >  NC_002947/5325370‑5325659                                                                                                                                              |                                                                                                                                                      gggCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCgg                                                                                                                                                <  2:362550/149‑1 (MQ=255)  ggCAAAGGGATTTACGTCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCggg                                                                                                                                               <  5:338547/149‑1 (MQ=255)  ggCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCggg                                                                                                                                               >  4:332553/1‑149 (MQ=255)      aaGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCGGG‑t                                                                                                                                             <  8:336115/146‑1 (MQ=255)      aaGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCGGG‑t                                                                                                                                             >  7:336115/1‑146 (MQ=255)       aGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                                 <  5:185004/142‑1 (MQ=255)       aGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                                 >  6:185004/1‑142 (MQ=255)           aTTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCttttttt                                                                                                                                     <  6:410469/149‑1 (MQ=255)           aTTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCttttttt                                                                                                                                     >  5:76329/1‑149 (MQ=255)                 ggatAAATGAGAATTGTTATTATTATAGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCtgt                                                                                                                               <  8:123968/146‑1 (MQ=255)                 ggCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCtgt                                                                                                                               >  7:43132/1‑149 (MQ=255)                   cTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTAGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCATGTGc                                                                                                                             <  5:284322/149‑1 (MQ=255)                          agaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                      <  1:296911/149‑1 (MQ=255)                           gaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                      <  7:331479/148‑1 (MQ=255)                           gaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGAAAAGCTAAGAGCCCCTGGGTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                      <  3:290011/148‑1 (MQ=255)                             aTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGCTCTTTTTTTGCCTGTGCCTTTTCAgg                                                                                                                    <  1:5187/149‑1 (MQ=255)                              ttGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGAt                                                                                                                  >  3:89929/1‑149 (MQ=255)                                 ttattattaTCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                      >  6:265002/1‑142 (MQ=255)                                 ttattattaTCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                      <  5:265002/142‑1 (MQ=255)                                 ttattattaTCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAg                                                                                                               >  5:335946/1‑149 (MQ=255)                                    ttattatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAgcgc                                                                                                            >  2:42278/1‑149 (MQ=255)                                      attatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgccg                                                                                                          >  6:350625/1‑149 (MQ=255)                                       ttatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgccg                                                                                                          >  5:20461/1‑148 (MQ=255)                                         atCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTc                                                                                                                       <  3:58766/133‑1 (MQ=255)                                         atCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTc                                                                                                                       >  4:58766/1‑133 (MQ=255)                                         atCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCa                                                                                                        <  5:350625/148‑1 (MQ=255)                                           cGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCt                                                                                                     >  8:218683/1‑149 (MQ=255)                                            gCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTg                                                                                                    <  7:47453/149‑1 (MQ=255)                                             cAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTgc                                                                                                   >  6:412408/1‑149 (MQ=255)                                                cTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTgcgc                                                                                                 >  8:291491/1‑148 (MQ=255)                                                     cgcgAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACaata                                                                                           <  3:332553/149‑1 (MQ=255)                                                                ttGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCg                                                                                <  5:343270/149‑1 (MQ=255)                                                                    aTAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCAt                                                                            <  8:43132/149‑1 (MQ=255)                                                                     tAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATc                                                                           >  8:138715/1‑149 (MQ=255)                                                                       aGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTTAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATcg                                                                          <  1:287338/148‑1 (MQ=255)                                                                           aaGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCgccgcc                                                                     >  4:317950/1‑149 (MQ=255)                                                                             gagCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgt                                                                   <  6:335946/149‑1 (MQ=255)                                                                             gagCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgt                                                                   <  1:42278/149‑1 (MQ=255)                                                                              agCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgtg                                                                  >  8:115033/1‑149 (MQ=255)                                                                               gCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgtgt                                                                 <  6:76329/149‑1 (MQ=255)                                                                                ccccTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTa                                                                <  4:89929/149‑1 (MQ=255)                                                                                 cccTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTAc                                                               <  7:291491/149‑1 (MQ=255)                                                                                   cTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACca                                                             <  7:218683/149‑1 (MQ=255)                                                                                       ttCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACcacaca                                                         <  5:28223/149‑1 (MQ=255)                                                                                          ggACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATaa                                                      <  6:20461/149‑1 (MQ=255)                                                                                            aCTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATc                                                    <  3:11402/149‑1 (MQ=255)                                                                                              tctcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGCG‑TCTTTTTTTGCCTGAGCCTTTTCAGGATCACCGCCGCAGCTGCGAACAATAATCCTTCTCCGGCATCGCCGCCGTGTACCACACATAATCAg                                                  >  5:71456/1‑149 (MQ=38)                                                                                                tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGcc                                                >  1:209392/1‑149 (MQ=255)                                                                                                tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGc                                                 >  5:144272/1‑148 (MQ=255)                                                                                                    aTTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGc                                            >  7:321274/1‑149 (MQ=255)                                                                                                     ttATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGcc                                           <  8:321274/149‑1 (MQ=255)                                                                                                           cctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTc                                     >  8:277448/1‑149 (MQ=255)                                                                                                            ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCg                                    >  8:361422/1‑149 (MQ=255)                                                                                                            ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCg                                    <  2:209392/149‑1 (MQ=255)                                                                                                             tcatcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCgg                                   >  6:11323/1‑149 (MQ=255)                                                                                                             tcatcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCgg                                   >  2:220781/1‑149 (MQ=255)                                                                                                                  aGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgc                              <  8:8089/149‑1 (MQ=255)                                                                                                                  aGGCTAATCACGGTTATTGACCCGGCAATTTGACGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgc                              >  1:254913/1‑149 (MQ=255)                                                                                                                    gCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctc                            <  5:412408/149‑1 (MQ=255)                                                                                                                    gCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctc                            >  5:90005/1‑149 (MQ=255)                                                                                                                     cTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctcg                           <  3:317950/149‑1 (MQ=255)                                                                                                                              gTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTca                  <  7:277448/149‑1 (MQ=255)                                                                                                                              aaaaaagACCCGGCAAATTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCt                    >  1:263474/7‑147 (MQ=255)                                                                                                                                     aCCCGGCAATTTGCCGGG‑TCTTTTTTTGCCCGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCCGGCGCGGTTCGCGCGCCCACCTCAGCcagca           >  8:146914/1‑149 (MQ=255)                                                                                                                                       ccGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCagcagc         <  6:90005/149‑1 (MQ=255)                                                                                                                                        cGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCagcagc         >  2:1875/1‑148 (MQ=255)                                                                                                                                          gCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAAt      <  1:1875/149‑1 (MQ=255)                                                                                                                                             aTTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATAcc   <  7:138715/149‑1 (MQ=255)                                                                                                                                              tttGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCt  >  8:371701/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                      GGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGG‑TCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCT  >  NC_002947/5325370‑5325659

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