Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,918	T→C	9.7%	noncoding (115/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,918	0	T	C	9.7%	195.6 / 9.5	73	noncoding (115/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base C (3/4);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

AGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGA  >  NC_002947/450776‑451057                                                                                                                                               |                                                                                                                                            aGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcg                                                                                                                                       >  8:160798/1‑149 (MQ=255)  gATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggcgg                                                                                                                                      <  5:225239/149‑1 (MQ=255)   aTTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                     <  3:265844/149‑1 (MQ=255)    ttGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCgggg                                                                                                                                    >  7:301083/1‑149 (MQ=255)       gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGttt                                                                                                                                 <  5:141985/149‑1 (MQ=255)       gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGttt                                                                                                                                 >  1:92426/1‑149 (MQ=255)       gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGttt                                                                                                                                 <  1:372283/149‑1 (MQ=255)       gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                 >  2:341893/1‑149 (MQ=255)         tCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTc                                                                                                                               <  7:295642/149‑1 (MQ=255)           cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTcca                                                                                                                             <  7:349932/149‑3 (MQ=255)            aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAATGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGTATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                            <  5:389835/149‑1 (MQ=255)            aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTc                                                                                                                               >  8:75580/1‑146 (MQ=255)             cTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTccaga                                                                                                                           <  7:75580/149‑5 (MQ=255)              tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                     >  3:250670/1‑138 (MQ=255)              tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                     <  4:250670/138‑1 (MQ=255)                   aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTAGAAAACCCCGCCGAAGCGGGGTTTTcca                                                                                                                             >  2:290255/1‑139 (MQ=255)                   aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGAAACGGATGTCACAGGGAAACAGTCTAGAAAACCCCGCCGAAGCGGGGTTTTcca                                                                                                                             <  1:290255/141‑3 (MQ=255)                    tCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                    >  8:135829/1‑149 (MQ=255)                      tttACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGAAGTCACAGGCAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCgtgt                                                                                                                  >  4:235222/1‑149 (MQ=255)                         aCCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttt                                                                                                                <  1:23175/148‑1 (MQ=255)                           cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTc                                                                                                             <  7:160798/149‑1 (MQ=255)                            acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                              >  8:213825/1‑147 (MQ=255)                            acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                            <  4:259660/149‑1 (MQ=255)                             caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgt                                                                                                                        >  2:309643/1‑136 (MQ=255)                             caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgt                                                                                                                        <  1:309643/136‑1 (MQ=255)                             caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                           <  2:92426/149‑1 (MQ=255)                             caTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                           <  1:378343/149‑1 (MQ=255)                               tGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGc                                                                                                         >  3:179143/1‑149 (MQ=255)                                 gACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCa                                                                                                       <  4:179143/149‑1 (MQ=255)                                  aCGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAg                                                                                                      <  3:323197/149‑1 (MQ=255)                                     ttGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgcg                                                                                                                     <  5:300980/131‑1 (MQ=255)                                     ttGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgcg                                                                                                                     >  6:300980/1‑131 (MQ=255)                                             gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCg                                                                                           >  4:125933/1‑149 (MQ=255)                                             gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCg                                                                                           <  4:169857/149‑1 (MQ=255)                                              aaGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCgg                                                                                          <  1:173837/149‑1 (MQ=255)                                                      aGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTc                                                                                  <  3:125933/149‑1 (MQ=255)                                                        gACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgc                                                                                                                          >  1:227597/1‑107 (MQ=255)                                                        gACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTtgc                                                                                                                          <  2:227597/107‑1 (MQ=255)                                                        gACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATg                                                                                     <  4:136684/144‑1 (MQ=255)                                                        gACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATg                                                                                     >  3:136684/1‑144 (MQ=255)                                                        gACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTcc                                                                                 >  1:42867/1‑148 (MQ=255)                                                                         gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgt                                                               <  4:324601/149‑1 (MQ=255)                                                                         gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgt                                                               >  7:81096/1‑149 (MQ=255)                                                                           cgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCa                                                                            >  4:345301/1‑134 (MQ=255)                                                                           cgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCa                                                                            <  3:345301/134‑1 (MQ=255)                                                                            gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGtctc                                                            >  3:365198/1‑149 (MQ=255)                                                                               ggACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtc                                                         <  2:50401/149‑1 (MQ=255)                                                                                  ccaccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcctc                                                      >  6:322207/1‑149 (MQ=255)                                                                                        gCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                            >  1:238920/1‑89 (MQ=255)                                                                                        gCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                            <  2:238920/89‑1 (MQ=255)                                                                                         ccAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGtctc                                                            >  2:134653/1‑136 (MQ=255)                                                                                         ccAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAAc                                               <  8:301083/149‑1 (MQ=255)                                                                                         ccAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAAc                                               >  4:376396/1‑149 (MQ=255)                                                                                              aTGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTgcg                                          >  5:100894/1‑149 (MQ=255)                                                                                                  cGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagc                                      <  4:365198/149‑1 (MQ=255)                                                                                                     ggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCagcagc                                   <  7:135829/149‑1 (MQ=255)                                                                                                      gCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcag                                       <  3:174691/144‑1 (MQ=255)                                                                                                      gCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcag                                       >  4:174691/1‑144 (MQ=255)                                                                                                        gATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTc                                 >  7:23280/1‑148 (MQ=255)                                                                                                         aTACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCca                               >  7:182089/1‑149 (MQ=255)                                                                                                              gATGTCACAGGGAAACAGTCTGGAAAACCCCGCATCGGCCGGGTATTCGTTGCGTGCGTGTTTTTCAGGCCAGCACATCCACCGGAGATGCTCGCATCAGGCGAGCGTCCGTCTCTTCCTCGAGCAACGTGCGCAGCAGCTCCAAagag                          >  6:159204/1‑149 (MQ=255)                                                                                                                    acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTGTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCAATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGAAGCTCCACAGAGGCctgc                    >  3:81152/1‑149 (MQ=255)                                                                                                                    acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgg                    >  5:333979/1‑148 (MQ=255)                                                                                                                    acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGAGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgc                    >  7:82484/1‑149 (MQ=255)                                                                                                                          aaaCAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgcgc              >  2:296072/1‑149 (MQ=255)                                                                                                                             cAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg           <  6:100894/149‑1 (MQ=255)                                                                                                                             cAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg           <  8:182089/149‑1 (MQ=255)                                                                                                                             cAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg           <  8:81096/149‑1 (MQ=255)                                                                                                                                    gAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCac                                                                           <  4:28373/78‑2 (MQ=255)                                                                                                                                    gAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCac                                                                           >  3:28373/1‑77 (MQ=255)                                                                                                                                      aaaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGa  <  1:296072/149‑1 (MQ=255)                                                                                                                                           ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                <  6:298883/50‑1 (MQ=255)                                                                                                                                           ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                >  5:298883/1‑50 (MQ=255)                                                                                                                                               |                                                                                                                                            AGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGA  >  NC_002947/450776‑451057

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