Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,325,511	T→A	10.8%	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	10.8%	168.5 / 11.3	65	noncoding (133/149 nt)	PP_mr52	PrrF ncRNA
Reads supporting (aligned to +/- strand):  ref base T (33/25);  new base A (3/4);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

GGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCT  >  NC_002947/5325369‑5325659                                                                                                                                               |                                                                                                                                                       ggggCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGcc                                                                                                                                                   <  5:216911/148‑1 (MQ=255) ggggCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                                  <  7:282568/149‑1 (MQ=255)    gCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCGGGt                                                                                                                                               <  4:287960/149‑1 (MQ=255)     cAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTc                                                                                                                                              <  2:265698/149‑1 (MQ=255)        aGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCGGGTcaat                                                                                                                                           <  5:346711/149‑4 (MQ=255)                   gCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCtgtg                                                                                                                                >  3:108299/1‑149 (MQ=255)                      aaaTGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGccg                                                                                                                             >  6:351120/1‑148 (MQ=255)                           agaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCtttt                                                                                                                          >  8:226681/1‑147 (MQ=255)                           agaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTCCCGGGTCTTTTTTTGCCTGTGCCTTTTc                                                                                                                         <  1:220221/148‑1 (MQ=255)                             aaTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAg                                                                                                                       >  3:27203/1‑148 (MQ=255)                              aTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGa                                                                                                                     <  3:206097/149‑1 (MQ=255)                                  ttattattaTCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCa                                                                                                                  <  6:214362/148‑1 (MQ=255)                                     ttattatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAgcg                                                                                                               <  4:27203/148‑1 (MQ=255)                                        ttatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGAATCTCAGATTATCTCCTCATCAGGCTAATAACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgccg                                                                                                            >  7:247383/1‑148 (MQ=255)                                          atCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGc                                                                                                                                     <  4:332497/121‑1 (MQ=255)                                          atCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGc                                                                                                                                     >  3:332497/1‑121 (MQ=255)                                           ttgCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGc                                                                                                        >  3:262070/3‑149 (MQ=255)                                             gCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTg                                                                                                      >  1:381342/1‑149 (MQ=255)                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCa                                                                                                                  <  7:210267/125‑1 (MQ=255)                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCa                                                                                                                  >  8:210267/1‑125 (MQ=255)                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑Tc                                                                                        <  8:14358/149‑1 (MQ=255)                                                          agaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑Tcc                                                                                       <  4:177705/149‑1 (MQ=255)                                                           gaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCt                                                                                      <  4:108299/149‑1 (MQ=255)                                                                gTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCtac                                                                                 >  5:173727/1‑147 (MQ=255)                                                                gTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTcc                                                                                 >  2:314087/1‑149 (MQ=255)                                                                 ttGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCg                                                                                <  7:226681/149‑1 (MQ=255)                                                                     aTAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGACTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCGCA‑ATAA‑TCCTGCTCCGGCAt                                                                            >  4:281682/1‑149 (MQ=255)                                                                      tAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATc                                                                           >  5:336131/1‑149 (MQ=255)                                                                      tAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATc                                                                           >  6:206790/1‑149 (MQ=255)                                                                      tAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATc                                                                           >  4:302106/1‑149 (MQ=255)                                                                         gCTAAGAGCCCCTGTTTCGTTCTCTCAGATTATCTCCTCATCAGGCTAATGTCGGTTATTGACCCGGCAATTTGCCTGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCGGCAGCTGCGCACA‑ATAA‑TCGTGCTCCGGCATcgcc                                                                        <  7:31680/149‑1 (MQ=38)                                                                         gCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATcgcc                                                                        >  7:131469/1‑149 (MQ=255)                                                                           tAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCGCGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATcgccgc                                                                      >  8:235855/1‑149 (MQ=255)                                                                              gagCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCgt                                                                   >  3:190486/1‑149 (MQ=255)                                                                                           ggACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATa                                                       <  5:206790/148‑1 (MQ=255)                                                                                           ggACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATa                                                       <  3:302106/148‑1 (MQ=255)                                                                                                 tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGcc                                                >  7:17915/1‑149 (MQ=255)                                                                                                   aGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCATTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTc                                              <  1:124590/149‑1 (MQ=255)                                                                                                    gATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCg                                             >  1:115692/1‑149 (MQ=255)                                                                                                            cctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTc                                     >  2:247298/1‑149 (MQ=255)                                                                                                             ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCg                                    <  8:194689/149‑1 (MQ=255)                                                                                                             ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCg                                    >  1:342500/1‑149 (MQ=255)                                                                                                              tcatcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCgg                                   >  6:334470/1‑149 (MQ=255)                                                                                                                atcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAAt                                                     >  5:100906/1‑129 (MQ=255)                                                                                                                atcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAAt                                                     <  6:100906/129‑1 (MQ=255)                                                                                                                atcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGtt                                 >  6:149716/1‑149 (MQ=255)                                                                                                                     gCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCGTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctc                            >  5:279255/1‑149 (MQ=255)                                                                                                                      cTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctcg                           >  1:172532/1‑149 (MQ=255)                                                                                                                         aTCACGGTTATTGACCCGGCAATTTTCCGGGTCTTTTTTTTCCGGTGCCCTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAATTCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGcc                         <  1:314087/149‑1 (MQ=255)                                                                                                                            gcGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCAc                      <  7:235855/147‑1 (MQ=255)                                                                                                                            aCGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATc                                                    >  7:11536/1‑118 (MQ=255)                                                                                                                            aCGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATc                                                    <  8:11536/118‑1 (MQ=255)                                                                                                                              ggTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCt                    >  1:105159/1‑148 (MQ=255)                                                                                                                               gTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTca                  <  8:247383/149‑1 (MQ=255)                                                                                                                                ttATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTcag                 <  8:131469/149‑1 (MQ=255)                                                                                                                                   ttGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGACTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTcagc                >  2:352328/1‑147 (MQ=255)                                                                                                                                   aagACCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCt                                                                                      >  6:332208/3‑77 (MQ=255)                                                                                                                                   aagACCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCt                                                                                      >  6:205774/3‑77 (MQ=255)                                                                                                                                   aagACCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCt                                                                                      <  5:205774/75‑1 (MQ=255)                                                                                                                                   aagACCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCt                                                                                      <  5:332208/75‑1 (MQ=255)                                                                                                                                     gACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACcaca                                                           >  8:331736/1‑102 (MQ=255)                                                                                                                                     gACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACcaca                                                           <  7:331736/102‑1 (MQ=255)                                                                                                                                        ccGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCcagcag          >  7:38435/1‑148 (MQ=255)                                                                                                                                         cGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAGTTAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCcagcag          >  1:314643/1‑148 (MQ=255)                                                                                                                                          ggCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCaa       >  1:64876/1‑149 (MQ=255)                                                                                                                                               tttGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCt  <  2:105159/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                       GGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACA‑ATAA‑TCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCT  >  NC_002947/5325369‑5325659

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