Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,543,124	A→C	6.8%	A45A (GCT→GCG)	rplI ←	50S ribosomal protein L9

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,543,124	0	A	C	6.8%	179.6 / 4.2	59	A45A (GCT→GCG)	rplI	50S ribosomal protein L9
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base C (2/2);  total (30/29)
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.15e-01

GAACCGAACAGCTTGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTCCATTTGGTAAAACCT  >  NC_002947/5542984‑5543270                                                                                                                                                |                                                                                                                                                   gaaccgaacAGCTTGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCCGCCaacct                                                                                                                                            >  6:350191/1‑145 (MQ=255)  aaccgaacAGCTTGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTg                                                                                                                                           >  1:63766/1‑149 (MQ=255)     cgaacAGCTTGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGCTGGGGGGCTTTAAGCCGAAGGTTTATTGTAAGCGGCTGCTTTTTACAGCTCGGCACGACGCAATTCGAACGCAGCCAGGTTg                                                                                                                                           >  1:108754/1‑149 (MQ=255)          aGCTTGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCCGCCAACCTGGCTGCGtt                                                                                                                                   <  3:140494/149‑1 (MQ=255)            cTTGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCCGCCAACCTGGCTGCGttcg                                                                                                                                 <  1:87299/149‑3 (MQ=255)             ttGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCAGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTgg                                                                                                                                 >  6:119594/1‑148 (MQ=255)                 cTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACg                                                                                                                            >  1:326973/1‑149 (MQ=255)                         ccAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCtt                                                                                                                    >  5:300988/1‑149 (MQ=255)                           aGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGc                                                                                                                  <  2:43564/149‑1 (MQ=255)                               gTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGaa                                                                                                              >  8:285034/1‑149 (MQ=255)                               gTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGaa                                                                                                              <  7:285034/149‑1 (MQ=255)                                tggctttGATGCTCACTTCCAGGTGCGCGAGTTGGGGTGC‑‑‑GCGGCTTTGAGCGGAAGCTTTCTTGTCAGCGGCTGATTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGGCTTGCCGAAt                                                                                                             <  4:324875/143‑1 (MQ=255)                                  gCAGTGATGGTCACTTCCAGCTCGGCCAGTTGTGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGAGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATgg                                                                                                           <  6:300988/149‑1 (MQ=255)                                  gCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATgg                                                                                                           <  2:326973/149‑1 (MQ=255)                                  gCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATgg                                                                                                           >  6:233587/1‑149 (MQ=255)                                   cAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGgc                                                                                                          <  5:233587/149‑1 (MQ=255)                                    aGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATgg                                                                                                           >  3:107360/1‑147 (MQ=255)                                           gTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAgg                                                                                                    >  6:356922/1‑147 (MQ=255)                                              aCTTCCAGCTCGGCCAGTTTTTCAGC‑‑‑TCGGCTTTCAGCCGAAGCTTTCTTCTCATCGGATGCTTTTTCCAGCTCTTCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGtt                                                                                               <  2:63766/149‑1 (MQ=255)                                                ttCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTAc                                                                                             >  3:1365/1‑149 (MQ=255)                                                            cAGTTTGGCTGA‑‑‑GCGTCTTTCAGCCGAAGCTTTCTTGTCAGGGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAg                                                                                 <  6:221768/149‑1 (MQ=255)                                                             aGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGc                                                                                <  5:356922/149‑1 (MQ=255)                                                             aGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCAAGACGCTCTTCGAACGCAGCAAGGTTGGCGGCTTTGGGAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTAAGGCCGTAAcctg                                                                                 >  5:180145/1‑146 (MQ=255)                                                              gTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGcc                                                                               <  4:102531/149‑1 (MQ=255)                                                              gTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGcc                                                                               <  2:108754/149‑1 (MQ=255)                                                                 gggCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCtt                                                                             >  7:77209/1‑148 (MQ=255)                                                                    cAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACt                                                                         <  8:77209/149‑1 (MQ=255)                                                                      gc‑‑‑gcGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACttt                                                                       <  4:107360/149‑1 (MQ=255)                                                                       c‑‑‑gcgGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGaa                                                                                                              >  2:255954/1‑109 (MQ=255)                                                                       c‑‑‑gcgGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGaa                                                                                                              <  1:255954/109‑1 (MQ=255)                                                                                   cAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACgg                                                                                           <  2:312753/119‑1 (MQ=255)                                                                                   cAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACgg                                                                                           >  1:312753/1‑119 (MQ=255)                                                                                                 ttGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAGCGCAGCCAGGTTTGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCgcccaggttgcccag                                               >  8:324978/1‑149 (MQ=255)                                                                                                   gTCAGCGGCTGCTTTTTCCAGCTCAGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCgcccaggttgcccag                                               <  5:119594/147‑1 (MQ=255)                                                                                                         ggCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGa                                       >  3:227679/1‑149 (MQ=255)                                                                                                           ctgctTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACt                                     >  3:105082/1‑149 (MQ=255)                                                                                                              ctTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACtttt                                  >  3:191758/1‑149 (MQ=255)                                                                                                              ctTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACtttt                                  >  1:192236/1‑149 (MQ=255)                                                                                                              ctTTTTCCAGCTCGGCACGACCCGCTTCGAACGCAGCCATGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACATTTACTTTGTCGCCCATGTTGCCCACGTTACCGACtttt                                  <  4:105082/149‑1 (MQ=255)                                                                                                               tttttCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACtttt                                                                      <  1:152362/112‑1 (MQ=255)                                                                                                               tttttCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACtttt                                                                      >  2:152362/1‑112 (MQ=255)                                                                                                               tttttCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACttt                                                                 >  2:236462/1‑117 (MQ=255)                                                                                                               tttttCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACttt                                                                 <  1:236462/117‑1 (MQ=255)                                                                                                                 tttCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTc                                <  4:227679/148‑1 (MQ=255)                                                                                                                   tCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCcag                             >  1:66078/1‑149 (MQ=255)                                                                                                                   tCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCca                              <  4:191758/148‑1 (MQ=255)                                                                                                                         tCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGc                                         >  2:231635/1‑131 (MQ=255)                                                                                                                         tCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGc                                         <  1:231635/131‑1 (MQ=255)                                                                                                                           ggCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGAt                       <  6:217741/147‑1 (MQ=255)                                                                                                                                   gCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTcc                               >  6:54916/1‑131 (MQ=255)                                                                                                                                   gCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTcc                               <  5:54916/131‑1 (MQ=255)                                                                                                                                   gCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTCCAtt             >  6:33732/1‑149 (MQ=255)                                                                                                                                     tctTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATc                      >  7:327061/1‑138 (MQ=255)                                                                                                                                     tctTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATc                      <  8:327061/138‑1 (MQ=255)                                                                                                                                     tctTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTcc                <  3:336526/144‑1 (MQ=255)                                                                                                                                     tctTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTcc                >  4:336526/1‑144 (MQ=255)                                                                                                                                       ttCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTCCATTTgg          >  3:313922/1‑148 (MQ=255)                                                                                                                                         ccaaCGCCGCCAACCTGGCTGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTCCATTTGGTa        >  4:209893/3‑148 (MQ=255)                                                                                                                                              gCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTCCATTTGGTAAAACct  <  4:1365/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GAACCGAACAGCTTGCCTTCGTCGCCAGCGGTGGCAGTGATGGTCACTTCCAGCTCGGCCAGTTGGGCAGC‑‑‑GCGGCTTTCAGCCGAAGCTTTCTTGTCAGCGGCTGCTTTTTCCAGCTCGGCACGACGCTCTTCGAACGCAGCCAGGTTGGCGGCGTTGGCAACGGTGGCCTTGCCGAATGGCAGCAGGAAGTTACGGCCGTAACCAGCCTTAACTTTTACTTTGTCGCCCAGGTTGCCCAGGTTAGCGACTTTTTCCAGCAGGATCAGTTCCATTTGGTAAAACCT  >  NC_002947/5542984‑5543270

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