Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,921	G→A	6.9%	noncoding (118/127 nt)	PP_mr05 →	RsmY ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	450,921	0	G	A	6.9%	169.9 / 2.4	58	noncoding (118/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base G (21/33);  new base A (2/2);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

TAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCA  >  NC_002947/450774‑451064                                                                                                                                                    |                                                                                                                                                tAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcgg                                                                                                                                                >  4:334405/1‑149 (MQ=255)   aGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcg                                                                                                                                              >  6:369865/1‑149 (MQ=255)     aTTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                            <  5:152589/149‑1 (MQ=255)      ttGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                            <  6:76988/148‑1 (MQ=255)          aTCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGtttt                                                                                                                                       <  5:25183/149‑1 (MQ=255)           tCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTc                                                                                                                                      <  1:189590/149‑1 (MQ=255)             cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCtt                                                                                                                                    <  6:165811/149‑1 (MQ=255)             cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCtt                                                                                                                                    >  5:178522/1‑149 (MQ=255)              aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                   <  4:24140/149‑1 (MQ=255)                  ttAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGtttt                                                                                                                                       >  5:131697/1‑141 (MQ=255)                  ttAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGtttt                                                                                                                                       <  6:131697/141‑1 (MQ=255)                     aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgcg                                                                                                                            >  1:347373/1‑149 (MQ=255)                              acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                     >  2:340059/1‑147 (MQ=255)                              acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                   <  6:120789/149‑1 (MQ=255)                              acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                   <  6:107232/149‑1 (MQ=255)                               caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                  >  5:252995/1‑149 (MQ=255)                               caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                  <  3:269982/149‑1 (MQ=255)                               caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                  <  3:98286/149‑1 (MQ=255)                                aTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAgg                                                                                                                 >  4:253354/1‑149 (MQ=255)                                     cGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGc                                                                                                            <  6:178522/149‑1 (MQ=255)                                      gTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGca                                                                                                           >  5:343103/1‑149 (MQ=255)                                           gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                       <  2:53593/148‑1 (MQ=255)                                           gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                       >  1:53593/1‑148 (MQ=255)                                                              gATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCa                                                                                   <  3:334405/149‑1 (MQ=255)                                                              gATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCa                                                                                   <  6:6418/149‑1 (MQ=255)                                                                     gacgacgacGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGc                                                                            <  5:369865/149‑1 (MQ=255)                                                                               aaGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                  <  5:305947/149‑1 (MQ=255)                                                                               aaGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                  <  3:253354/149‑1 (MQ=255)                                                                               aaGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                  >  7:280404/1‑149 (MQ=255)                                                                                aGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                  >  4:183426/1‑148 (MQ=255)                                                                                 ggACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtc                                                                >  2:5760/1‑149 (MQ=255)                                                                                   accaccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcct                                                              <  6:343103/149‑1 (MQ=255)                                                                                    ccaccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcctc                                                             <  6:380332/149‑1 (MQ=255)                                                                                     caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcctc                                                             >  8:15732/1‑148 (MQ=255)                                                                                      accTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCa                                                           <  2:347373/149‑1 (MQ=255)                                                                                           ccAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAAc                                                      >  8:480/1‑149 (MQ=255)                                                                                             aGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGt                                                    <  3:183426/149‑1 (MQ=255)                                                                                              ggATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTg                                                   <  4:256293/149‑1 (MQ=255)                                                                                                   gCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcag                                              <  7:15732/149‑1 (MQ=255)                                                                                                    cGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgca                                               >  4:207900/1‑147 (MQ=255)                                                                                                    cGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgca                                               <  3:207900/147‑1 (MQ=255)                                                                                                         cGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTc                                        <  7:480/149‑1 (MQ=255)                                                                                                               ggATGTCACAGGGAAACAGTCTGGAAAACCCCGCTACGGCGGGGTTTTCTTTGCGTGAGTGTTTTTCAGGCCAGCACATCCAGCGGAGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCcaca                                    >  6:224467/1‑147 (MQ=255)                                                                                                                gATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                                 <  6:252995/149‑1 (MQ=255)                                                                                                                 aTGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAgg                                <  8:280404/149‑1 (MQ=255)                                                                                                                      acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCca                                      >  2:7169/1‑138 (MQ=255)                                                                                                                      acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCca                                      <  1:7169/138‑1 (MQ=255)                                                                                                                          ggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTcc                                                                                        <  6:380509/84‑1 (MQ=255)                                                                                                                          ggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgc                       >  4:296535/1‑149 (MQ=255)                                                                                                                               cAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg                  <  1:5760/149‑1 (MQ=255)                                                                                                                                      gAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcgc           >  5:279738/1‑149 (MQ=255)                                                                                                                                      gAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATc             <  1:340059/147‑1 (MQ=255)                                                                                                                                      gAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcgc                     <  4:244692/139‑1 (MQ=255)                                                                                                                                      gAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcgc                     >  3:244692/1‑139 (MQ=255)                                                                                                                                        aaaCCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcgcg          >  2:157086/1‑148 (MQ=255)                                                                                                                                          aCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaaa       <  6:8019/149‑1 (MQ=255)                                                                                                                                            cccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATa     >  3:299631/1‑149 (MQ=255)                                                                                                                                             ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAc    >  1:13320/1‑149 (MQ=255)                                                                                                                                              cGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAc    >  4:188374/1‑148 (MQ=255)                                                                                                                                               gCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCa  <  6:279738/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCA  >  NC_002947/450774‑451064

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