Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,919	T→G	5.2%	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.2%	284.8 / 4.3	92	noncoding (116/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base T (40/47);  new base G (3/2);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCA  >  NC_002947/450775‑451064                                                                                                                                                   |                                                                                                                                                  aaGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑Tcggc                                                                                                                                               >  8:130485/1‑149 (MQ=255)   gATTGGATCCACTGTTAATCTTTACCTCATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑Tcggcgg                                                                                                                                             <  4:2415/149‑1 (MQ=255)   gATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑Tcggcgg                                                                                                                                             <  3:328198/149‑1 (MQ=255)     ttGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCggg                                                                                                                                            <  6:177510/148‑1 (MQ=255)        gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGttt                                                                                                                                        >  6:70719/1‑149 (MQ=255)         aTCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCC‑‑GAAGCGGGGtttt                                                                                                                                       <  7:211188/149‑1 (MQ=255)          tCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTc                                                                                                                                      <  3:334369/149‑1 (MQ=255)            cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCtt                                                                                                                                    <  6:283923/149‑1 (MQ=255)            cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCtt                                                                                                                                    >  4:78151/1‑149 (MQ=255)            cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGAAACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCtt                                                                                                                                    >  7:223279/1‑149 (MQ=255)             aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCttt                                                                                                                                   <  4:406718/149‑1 (MQ=255)             aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGAAAAACCCCGCC‑‑GAAGCGGGGTTTTcc                                                                                                                                     <  1:138357/147‑2 (MQ=255)             aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGAAAAACCCCGCC‑‑GAAGCGGGGTTTTcc                                                                                                                                     >  2:138357/1‑146 (MQ=255)               tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCttt                                                                                                                                   <  3:74787/147‑1 (MQ=255)               tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCttt                                                                                                                                   >  4:74787/1‑147 (MQ=255)               tGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgc                                                                                                                                 >  1:155127/1‑149 (MQ=255)                 ttAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgcgt                                                                                                                               >  5:217891/1‑149 (MQ=255)                 ttAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgcg                                                                                                                                <  5:188650/148‑1 (MQ=255)                    aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgcgtgc                                                                                                                             <  8:5767/148‑1 (MQ=255)                       tttACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgc                                                                                                                                 <  7:115256/141‑1 (MQ=255)                       tttACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgc                                                                                                                                 >  8:115256/1‑141 (MQ=255)                         tACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                           <  5:93092/145‑1 (MQ=255)                         tACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                           <  3:78151/145‑1 (MQ=255)                         tACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGttt                                                                                                                       >  3:273852/1‑149 (MQ=255)                           gcctgtGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                     <  5:54328/144‑1 (MQ=255)                                  gACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCa                                                                                                              <  8:223279/149‑1 (MQ=255)                                   aCGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgc                                                                                                                                 >  3:110660/1‑129 (MQ=255)                                   aCGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTtgc                                                                                                                                 <  4:110660/129‑1 (MQ=255)                                     gTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                  <  8:51363/142‑1 (MQ=255)                                     gTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAg                                                                                                                  >  7:51363/1‑142 (MQ=255)                                          gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTc                                                                                                                    >  2:267272/1‑135 (MQ=255)                                          gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTc                                                                                                                    <  1:267272/135‑1 (MQ=255)                                          gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACCTc                                                                                                       >  4:243905/1‑148 (MQ=255)                                          gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                       >  1:373033/1‑148 (MQ=255)                                          gcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATc                                                                                                       >  3:181525/1‑148 (MQ=255)                                            aGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAg                                                                                                    <  7:130485/149‑1 (MQ=255)                                                   gcAAAGGACAGATCAAGGAGGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAACCAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGAt                                                                                             >  2:404499/1‑149 (MQ=255)                                                   gcAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGGT‑‑GCGCCGGGGTTTTATTTGCGTCCGGGTTTTTCAGGCAAGCCCTACCAGCGTTGAt                                                                                             >  2:146917/1‑149 (MQ=255)                                                   gcAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGAt                                                                                             >  7:33907/1‑149 (MQ=255)                                                    cAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATg                                                                                            >  3:363433/1‑149 (MQ=255)                                                     aaaGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGc                                                                                           >  5:134687/1‑149 (MQ=255)                                                       aGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTc                                                                                         <  1:204673/149‑1 (MQ=255)                                                         gACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTG‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTcc                                                                                        <  1:323044/149‑1 (MQ=255)                                                             gATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACAACCAGCGGTGATGCTCCCAt                                                                                     >  2:79482/1‑147 (MQ=255)                                                               tCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGt                                                                                                <  4:299533/134‑1 (MQ=255)                                                               tCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGt                                                                                                >  3:299533/1‑134 (MQ=255)                                                                 aaGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGc                                                                                <  4:273852/148‑1 (MQ=255)                                                                        acgacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtcc                                                                        <  4:363433/149‑1 (MQ=255)                                                                            cgaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCAGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgtc                                                                     >  3:110915/1‑148 (MQ=255)                                                                             gaAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGtctc                                                                   >  7:384243/1‑149 (MQ=255)                                                                                ggACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtc                                                                >  6:317764/1‑149 (MQ=255)                                                                                    gaccTGCCTGGTTGGCGGGGCGATACGGATTTCAGAGGGAAACAGTCTGGAAAACCCCGAT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCATCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcctc                                                             <  4:110915/147‑1 (MQ=255)                                                                                    caccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCTAGGCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCAACCACATCCAGCGGTGATGATCCCATCAGGCGAGCGTCCGTCTCTtcct                                                              >  1:117436/1‑149 (MQ=255)                                                                                     accTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAACACATCCAGCGGt                                                                                                >  8:241982/1‑112 (MQ=255)                                                                                     accTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAACACATCCAGCGGt                                                                                                <  7:241982/112‑1 (MQ=255)                                                                                      ccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAg                                                          >  4:84942/1‑149 (MQ=255)                                                                                      ccTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAg                                                          <  8:33907/149‑1 (MQ=255)                                                                                          ccAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAAc                                                      >  3:305893/1‑149 (MQ=255)                                                                                           cAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCATCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACg                                                     <  5:70719/149‑1 (MQ=255)                                                                                            aGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGt                                                    <  2:155127/149‑1 (MQ=255)                                                                                                   cGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagc                                             <  2:117436/149‑1 (MQ=255)                                                                                                    ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTc                                                                                                                                      <  2:172568/59‑1 (MQ=255)                                                                                                    ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTgtgtgcgt                                                                                                                           <  4:46356/70‑1 (MQ=255)                                                                                                    ggggCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTgtgtgcgt                                                                                                                           >  3:46356/1‑70 (MQ=255)                                                                                                       gCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                  <  4:196779/124‑1 (MQ=255)                                                                                                       gCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                                  >  3:196779/1‑124 (MQ=255)                                                                                                          aTACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCca                                      <  4:367670/149‑1 (MQ=255)                                                                                                          aTACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCca                                      >  3:367670/1‑149 (MQ=255)                                                                                                            aCGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagc                                             >  1:190676/1‑140 (MQ=255)                                                                                                            aCGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCgcagc                                             <  2:190676/140‑1 (MQ=255)                                                                                                              ggATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCcaca                                    >  8:90657/1‑147 (MQ=255)                                                                                                               gATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                                 <  4:181525/149‑1 (MQ=255)                                                                                                                aTGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAgg                                >  7:222553/1‑149 (MQ=255)                                                                                                                   tCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCct                             >  3:288119/1‑149 (MQ=255)                                                                                                                     acaGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgc                           <  3:84942/149‑1 (MQ=255)                                                                                                                       aGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgctg                         <  5:317764/149‑1 (MQ=255)                                                                                                                        gggAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTgg                        >  5:160373/1‑149 (MQ=255)                                                                                                                        gggAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTgg                        <  4:288119/149‑1 (MQ=255)                                                                                                                              cAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg                  <  4:180639/149‑1 (MQ=255)                                                                                                                               aGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTggg                 <  8:384243/149‑1 (MQ=255)                                                                                                                                   tGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATc             >  1:361730/1‑149 (MQ=255)                                                                                                                                     gAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGtctc                                                                   >  8:87256/1‑93 (MQ=255)                                                                                                                                     gAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGtctc                                                                   <  7:87256/93‑1 (MQ=255)                                                                                                                                       aaaCCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCCCATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGa         <  1:404499/149‑1 (MQ=255)                                                                                                                                        aaCCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGa         >  3:115439/1‑148 (MQ=255)                                                                                                                                          ccccGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAAt      <  6:217891/149‑1 (MQ=255)                                                                                                                                          ccccGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAAt      <  7:110087/149‑1 (MQ=255)                                                                                                                                           cccGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATa     >  5:164068/1‑149 (MQ=255)                                                                                                                                            ccGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTggg                 >  2:287158/1‑136 (MQ=255)                                                                                                                                            ccGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTggg                 <  1:287158/136‑1 (MQ=255)                                                                                                                                             cGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATAcc   <  8:222553/149‑1 (MQ=255)                                                                                                                                              gCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCa  <  4:305893/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  AAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGATGTCACAGGGAAACAGTCTGGAAAACCCCGCT‑‑TCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGCACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCA  >  NC_002947/450775‑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: ‑