Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,919	T→G	5.5%	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	5.5%	291.3 / 5.7	91	noncoding (116/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base T (38/48);  new base G (3/2);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAC  >  NC_002947/450774‑451062                                                                                                                                                   |                                                                                                                                                tAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTcgg                                                                                                                                              >  7:73291/1‑149 (MQ=255)     aTTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCggg                                                                                                                                          <  6:121818/149‑1 (MQ=255)             cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCtt                                                                                                                                  <  5:26956/149‑1 (MQ=255)              aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCAGAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTGCttt                                                                                                                                 >  2:279485/1‑149 (MQ=255)                 gTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCCGAAGCGGGGTTTTcca                                                                                                                                  >  8:11868/1‑143 (MQ=255)                 gTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCCGAAGCGGGGTTTTcca                                                                                                                                  <  7:11868/145‑3 (MQ=255)                  ttAATCTTTACCACATGGACGTTGCGCAGGAAGCACAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTtgc                                                                                                                               <  6:151951/147‑1 (MQ=255)                     aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGGCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgcg                                                                                                                          >  8:177109/1‑149 (MQ=255)                     aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGtt                                                                                                                                       >  8:91954/1‑136 (MQ=255)                     aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGtt                                                                                                                                       <  7:91954/136‑1 (MQ=255)                     aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgcg                                                                                                                          >  1:41755/1‑149 (MQ=255)                        tttACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCgtgt                                                                                                                       >  8:248215/1‑149 (MQ=255)                          tACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                         <  6:242086/145‑1 (MQ=255)                            ccACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTtgcg                                                                                                                              <  1:173702/138‑1 (MQ=255)                            ccACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTtgcg                                                                                                                              >  2:173702/1‑138 (MQ=255)                            ccACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                   >  6:99904/1‑149 (MQ=255)                             cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                   >  6:111057/1‑148 (MQ=255)                             cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTc                                                                                                                  <  4:290893/149‑1 (MQ=255)                             cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTc                                                                                                                  <  1:208421/149‑1 (MQ=255)                             cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCCGAAGCGGGGtt                                                                                                                                       >  1:144302/1‑128 (MQ=255)                             cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCCGAAGCGGGGtt                                                                                                                                       <  2:144302/128‑1 (MQ=255)                                aTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAgg                                                                                                               <  8:149389/149‑1 (MQ=255)                                aTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTATGCGTGCGTGTTTTTCAgg                                                                                                               >  7:144658/1‑149 (MQ=255)                                 tGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGc                                                                                                              >  2:817/1‑149 (MQ=255)                                  ggACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGc                                                                                                              >  5:326121/1‑148 (MQ=255)                                          cgcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACAt                                                                                                      <  3:230828/148‑1 (MQ=255)                                          cgcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                     <  8:73291/149‑1 (MQ=255)                                               gtaGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgc                                                                                                                           <  1:232016/120‑1 (MQ=255)                                               gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTtgcgtgc                                                                                                                           >  2:232016/1‑122 (MQ=255)                                               gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                         >  5:428/1‑124 (MQ=255)                                               gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                         <  6:428/124‑1 (MQ=255)                                               gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCg                                                                                                >  8:119225/1‑149 (MQ=255)                                               gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCg                                                                                                <  6:326121/149‑1 (MQ=255)                                               gAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCg                                                                                                <  4:159843/149‑1 (MQ=255)                                                aaGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCgg                                                                                               <  6:26262/149‑1 (MQ=255)                                                     cAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCggg                                                                                                                                          >  6:153115/1‑101 (MQ=255)                                                     cAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCggg                                                                                                                                          <  5:153115/101‑1 (MQ=255)                                                      aaaGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGc                                                                                         >  2:215683/1‑149 (MQ=255)                                                        aGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTc                                                                                       <  3:257809/149‑1 (MQ=255)                                                        aGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTGTTCAGGCCAGCACATCCAGCGGTGATGCTc                                                                                       >  5:276948/1‑149 (MQ=255)                                                         ggACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTcc                                                                                      <  1:254549/149‑1 (MQ=255)                                                              gATCAAGGACGACGACGAAGGACCACCTGCCAGGATGTCGGGGCGATACGGATGTCACAGGGAAACAGTCTT‑GAAAACCCCGCTTCGGCTGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCa                                                                                 <  5:224759/149‑1 (MQ=255)                                                              gATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCa                                                                                 <  5:99904/149‑1 (MQ=255)                                                              gATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCa                                                                                 <  1:215683/149‑1 (MQ=255)                                                              gATAAAGGACGACGACGAAGGACCACCTGCCAGGATGGAGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACACCGCTTCGGCGGGGTTTTCTTAGCGTGCCTGTATTTCAGGCCAGCACATCCAGAGGTGATGCTACCAt                                                                                   >  1:249976/1‑147 (MQ=255)                                                                 cAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                 <  1:273054/114‑1 (MQ=255)                                                                 cAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                                 >  2:273054/1‑114 (MQ=255)                                                                    ggACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAg                                                                                <  4:175236/144‑1 (MQ=255)                                                                    ggACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAg                                                                                >  3:175236/1‑144 (MQ=255)                                                                      acgacgacgAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGGTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAg                                                                                >  4:26449/1‑142 (MQ=255)                                                                      acgacgacgAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGGTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAg                                                                                <  3:26449/142‑1 (MQ=255)                                                                        gacgacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                >  1:307423/1‑108 (MQ=255)                                                                        gacgacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                <  2:307423/108‑1 (MQ=255)                                                                              gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTccc                                                                                     >  4:155888/1‑129 (MQ=255)                                                                              gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTccc                                                                                     <  3:155888/129‑1 (MQ=255)                                                                                  gACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGCGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                 <  1:24111/82‑1 (MQ=255)                                                                                  gACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                                 >  2:24111/1‑81 (MQ=255)                                                                                   accaccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTc                                                                                       <  1:192370/122‑1 (MQ=255)                                                                                   accaccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTc                                                                                       >  2:192370/1‑122 (MQ=255)                                                                                    ccaccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcctc                                                           <  4:28542/149‑1 (MQ=255)                                                                                    ccaccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcctc                                                           <  8:44214/149‑1 (MQ=255)                                                                                     caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcctcc                                                          <  1:148960/149‑1 (MQ=255)                                                                                      accTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCa                                                         <  5:303152/149‑1 (MQ=255)                                                                                           ccAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAAc                                                    >  5:18269/1‑149 (MQ=255)                                                                                                      gggCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGcagcag                                         >  3:71676/1‑149 (MQ=255)                                                                                                            tACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATcc                                                                                                    >  3:290306/1‑84 (MQ=255)                                                                                                            tACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATcc                                                                                                    <  4:290306/84‑1 (MQ=255)                                                                                                               ggATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCcaca                                  >  7:186890/1‑147 (MQ=255)                                                                                                                gATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                               <  3:201310/149‑1 (MQ=255)                                                                                                                gATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                               <  5:111057/149‑1 (MQ=255)                                                                                                                    tCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagc                                           >  5:292925/1‑133 (MQ=255)                                                                                                                    tCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagc                                           <  6:292925/133‑1 (MQ=255)                                                                                                                     cacaGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctg                          >  7:239590/1‑149 (MQ=255)                                                                                                                     cacaGGGAAACAATCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctg                          >  6:259156/1‑149 (MQ=255)                                                                                                                      acaGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGc                                                                                         >  6:202866/1‑85 (MQ=255)                                                                                                                      acaGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGc                                                                                         <  5:202866/85‑1 (MQ=255)                                                                                                                      acaGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgc                         >  1:209566/1‑149 (MQ=255)                                                                                                                        aGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgctg                       >  5:210974/1‑149 (MQ=255)                                                                                                                          ggAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGgc                     >  6:111047/1‑149 (MQ=255)                                                                                                                             aaCAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCt                  >  6:292554/1‑149 (MQ=255)                                                                                                                             aaCAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCt                  <  1:817/149‑1 (MQ=255)                                                                                                                              aCAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTg                 <  3:42175/149‑1 (MQ=255)                                                                                                                              aCAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTg                 <  8:239590/149‑1 (MQ=255)                                                                                                                               cAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg                <  2:209566/149‑1 (MQ=255)                                                                                                                                      agAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcg          >  1:314280/2‑149 (MQ=255)                                                                                                                                       gAATACCCAGGTTCGGCGTGGTTTTGTTTGCGTGCGTGTTTTTCAGGCCAGCAGATCCAGCGGTGATGCTCCCATCAGGCGATCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcgc         <  8:203103/149‑1 (MQ=255)                                                                                                                                       gAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcgc         <  4:71676/149‑1 (MQ=255)                                                                                                                                          aaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaa      <  7:119225/149‑1 (MQ=255)                                                                                                                                          aaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaa      <  2:41755/149‑1 (MQ=255)                                                                                                                                           aCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaaa     <  7:247905/149‑1 (MQ=255)                                                                                                                                              ccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAc  <  6:18269/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                TAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTG‑GAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAC  >  NC_002947/450774‑451062

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