Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,919	T→G	6.7%	noncoding (116/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,919	0	T	G	6.7%	183.3 / 4.6	60	noncoding (116/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base T (22/33);  new base G (2/2);  total (25/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 = 9.40e-01

TTCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCA  >  NC_002947/450771‑451064                                                                                                                                                     |                                                                                                                                                  ttCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCtt                                                                                                                                                   >  4:269982/1‑149 (MQ=255)    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                                                                                                                                       <  6:131697/141‑1 (MQ=255)                     ttAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGtttt                                                                                                                                       >  5:131697/1‑141 (MQ=255)                        aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgcg                                                                                                                            >  1:347373/1‑149 (MQ=255)                                 acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                     >  2:340059/1‑147 (MQ=255)                                 acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                   <  6:107232/149‑1 (MQ=255)                                 acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                   <  6:120789/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                                                                  >  7:280404/1‑149 (MQ=255)                                                                                  aaGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                  <  3:253354/149‑1 (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                                               <  3:207900/147‑1 (MQ=255)                                                                                                       cGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgca                                               >  4:207900/1‑147 (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)                                                                                                                                                     |                                                                                                                                                  TTCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCA  >  NC_002947/450771‑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: ‑