Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,354,671	G→A	6.2%	T7I (ACC→ATC)  noncoding (4/97 nt)	rpsO ← PP_mr54 ←	30S ribosomal protein S15 S15 regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,354,671	0	G	A	6.2%	194.3 / 2.4	64	T7I (ACC→ATC)  noncoding (4/97 nt)	rpsO PP_mr54	30S ribosomal protein S15 S15 regulatory region
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base A (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CTTGTTGATGTTGGCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGCCGACGCGGCGCAATGCG  >  NC_002947/5354523‑5354819                                                                                                                                                     |                                                                                                                                                     cTTGTTGATGTTGGCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCgggg                                                                                                                                                      <  1:85204/149‑1 (MQ=255)    caggATGTTGGCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTc                                                                                                                                                    <  3:97340/145‑1 (MQ=255)     ttgATGTTGGCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGACCCCGCATGCGGGATccg                                                                                                                                                  <  6:323789/149‑2 (MQ=255)     ttgATGTTGGCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGACCCCGCATGCGGGATcc                                                                                                                                                   <  7:240246/148‑1 (MQ=255)         tGTTGGCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAg                                                                                                                                              >  6:24072/1‑149 (MQ=255)              gCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCa                                                                                                                                          <  2:143675/148‑1 (MQ=255)              gCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAg                                                                                                                                         <  8:35501/149‑1 (MQ=255)              gCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAg                                                                                                                                         <  2:231597/149‑1 (MQ=255)                ggTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATa                                                                                                                                               <  6:128205/141‑1 (MQ=255)                ggTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATa                                                                                                                                               >  5:128205/1‑141 (MQ=255)                  tCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGAc                                                                                                                                     <  6:177369/149‑1 (MQ=255)                  tCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGAc                                                                                                                                     <  7:273383/149‑1 (MQ=255)                  tCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGAc                                                                                                                                     <  7:161956/149‑1 (MQ=255)                      cagAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAAt                                                                                                                                 >  3:200523/1‑149 (MQ=255)                        gagCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATcc                                                                                                                               <  2:44402/149‑1 (MQ=255)                         agCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGa                                                                                                                                      >  1:69320/1‑141 (MQ=255)                         agCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGa                                                                                                                                      <  2:69320/141‑1 (MQ=255)                         agCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATccc                                                                                                                              <  5:222059/149‑1 (MQ=255)                          gCAACCTGCACTTGCGGGCTACCGGCATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGGTTTTTCTTCAACGCTGAGGCCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCt                                                                                                                             <  6:204218/149‑1 (MQ=255)                          gCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCt                                                                                                                             <  5:218676/149‑1 (MQ=255)                           cAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTg                                                                                                                            <  4:330660/149‑1 (MQ=255)                           cAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTg                                                                                                                            >  4:169290/1‑149 (MQ=255)                              ccTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACaa                                                                                                                                   <  3:315072/139‑1 (MQ=255)                              ccTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACaa                                                                                                                                   >  4:315072/1‑139 (MQ=255)                                tGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTAtt                                                                                                                        >  2:53658/1‑148 (MQ=255)                                 gCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATtaa                                                                                                                      >  4:307032/1‑149 (MQ=255)                                     ttCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCt                                                                                                                             <  2:217006/138‑1 (MQ=255)                                     ttCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCt                                                                                                                             >  1:217006/1‑138 (MQ=255)                                       ccGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATaaaa                                                                                                                 <  1:53658/148‑1 (MQ=255)                                           gCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGgt                                                                                                            <  7:128416/149‑1 (MQ=255)                                           gCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGgt                                                                                                            <  3:307032/149‑1 (MQ=255)                                               ccGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGgtgt                                                                                                          <  1:240230/147‑1 (MQ=255)                                                  gTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCg                                                                                                     >  3:71035/1‑149 (MQ=255)                                                   tATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGt                                                                                                    <  3:70784/149‑1 (MQ=255)                                                       gCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGcc                                                                                                 <  3:146527/148‑1 (MQ=255)                                                        ccGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGCGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTa                                                                                               <  2:87040/149‑1 (MQ=255)                                                                ttGCTGATATTCGGCAACGATCTGAGCTTCTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGc                                                                                       >  8:147709/1‑149 (MQ=255)                                                                      atatTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCt                                                                                 >  3:321386/1‑149 (MQ=255)                                                                      atatTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCt                                                                                 >  5:271802/1‑149 (MQ=255)                                                                               aaCGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATg                                                                            <  6:48728/145‑1 (MQ=255)                                                                               aaCGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATg                                                                            >  5:48728/1‑145 (MQ=255)                                                                                      tGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTc                                                                 >  2:300079/1‑149 (MQ=255)                                                                                        aGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCg                                                                <  4:200523/148‑1 (MQ=255)                                                                                          cTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGc                                                              >  3:948/1‑148 (MQ=255)                                                                                              ttcttcAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACt                                                          <  4:321386/148‑1 (MQ=255)                                                                                                 ttcAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCt                                                      >  5:310308/1‑149 (MQ=255)                                                                                                     aCGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGAcc                                                  <  3:169290/149‑1 (MQ=255)                                                                                                            gggCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGt                                           <  5:24072/149‑1 (MQ=255)                                                                                                                cATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCg                                       <  4:948/149‑1 (MQ=255)                                                                                                                cATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCg                                       >  2:275068/1‑149 (MQ=255)                                                                                                                cATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCg                                       <  7:147709/149‑1 (MQ=255)                                                                                                                 aTGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCgg                                      <  4:71035/149‑1 (MQ=255)                                                                                                                   gtgGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGt                                     >  7:242353/1‑148 (MQ=255)                                                                                                                        ttCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCATACTGACAGCCACCCTTGATGCAGTGGGGATCGGCAACTGACTGACCACGCGGTTTCGGTCATTc                                >  3:35675/1‑148 (MQ=255)                                                                                                                           tcCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTccgac                            >  6:126884/1‑149 (MQ=255)                                                                                                                            ccTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGAAAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTccgacc                           >  1:259887/1‑149 (MQ=255)                                                                                                                              tGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTccgaccg                          >  1:9103/1‑148 (MQ=255)                                                                                                                                 tattggACCCCGCATGGGGGATCCAATAGGCCAGGGACAAATGCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCCCCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAAt                       >  2:154132/5‑148 (MQ=255)                                                                                                                                 tattggACCCCGCATGCGGGATCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAAt                       >  4:116298/5‑148 (MQ=255)                                                                                                                                 tAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAAt                       >  4:48416/1‑148 (MQ=255)                                                                                                                                 tAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGATCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAAt                       >  1:168782/1‑148 (MQ=255)                                                                                                                                      gATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGcc                  >  4:178898/1‑148 (MQ=255)                                                                                                                                      gATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGcc                  >  1:277428/1‑148 (MQ=255)                                                                                                                                       aTCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGCCGa                >  1:221779/1‑149 (MQ=255)                                                                                                                                                gCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGCCGACGCGgcgc        >  1:326920/1‑148 (MQ=255)                                                                                                                                                     gTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGCCGACGCGGCGCAATgcg  >  5:309938/1‑149 (MQ=255)                                                                                                                                                     gTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGCCGACGCGGCGCAATgcg  <  1:147068/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CTTGTTGATGTTGGCGGTCAGCAGAGCAACCTGCACTTCCGGGCTACCGGTATCGCCGGCGGCTTGCTGATATTCGGCAACGATCTGAGCTTTTTCTTCAACGCTGAGGGCCATGTGGCTTCTCCTGATAACGGATCCCGCATGCGGGGTCCAATAGGCCAGGGACAAATCCCTGTATTAATAAAAAAGGTGTGACCGTGCCTACTGACAGCCACCCTTGATGCCGTGGGGCTCGGCAACTGCCTGACCCCGCGGTTTCGGTCATTCCGACCGAATCAGCCGACGCGGCGCAATGCG  >  NC_002947/5354523‑5354819

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