Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,919	T→G	8.8%	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	8.8%	196.3 / 10.8	68	noncoding (116/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base G (4/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

TTCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAA  >  NC_002947/450771‑451059                                                                                                                                                      |                                                                                                                                              ttCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGACTGGAAAACCCCGCtt                                                                                                                                               >  4:74441/1‑149 (MQ=255)    tAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcgg                                                                                                                                            >  8:219565/1‑149 (MQ=255)        aTTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                        >  5:27642/1‑149 (MQ=255)         ttGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                        <  6:107871/148‑1 (MQ=255)         ttGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                        <  7:17133/148‑1 (MQ=255)         ttGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCggg                                                                                                                                        <  2:126105/148‑1 (MQ=255)            gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGttt                                                                                                                                    >  8:307535/1‑149 (MQ=255)            gATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGttt                                                                                                                                    >  2:81280/1‑149 (MQ=255)             aTCCACTGTTAATCTTTACCACATGGACGTTGCGCATGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGtttt                                                                                                                                   <  7:307535/149‑1 (MQ=255)                cACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCtt                                                                                                                                <  7:98321/149‑1 (MQ=255)                 aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCGAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggc                                                                                                                                           >  2:68828/1‑137 (MQ=255)                 aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGGCATGTCACAGGGAAACAGTCTGGAAAACCCCGCCGAAGCGGGGTTTTcca                                                                                                                                <  1:81280/149‑3 (MQ=255)                 aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTcggc                                                                                                                                           <  1:68828/137‑1 (MQ=255)                 aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                               >  2:35112/1‑149 (MQ=255)                 aCTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTCGAAAACCCCGCTTCGGCGGGGTTTTCttt                                                                                                                               >  6:59476/1‑149 (MQ=255)                        aTCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGTACGACGACGAAGGACCACCTGCCAGGATGGCTGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGTGGTTTTCTTtgcgtgcg                                                                                                                        <  2:222039/149‑1 (MQ=255)                         tCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTgcgtgcgt                                                                                                                       >  2:21248/1‑149 (MQ=255)                               ccACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                 <  1:35112/149‑1 (MQ=255)                               ccACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                 <  3:336614/149‑1 (MQ=255)                                cacaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGttttt                                                                                                                 >  5:275912/1‑148 (MQ=255)                                 acaTGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCa                                                                                                               <  7:167172/149‑1 (MQ=255)                                    tGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGc                                                                                                            >  3:168848/1‑149 (MQ=255)                                      gACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCCGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCAGCTTCGGCGGGGTTTTCTTTGCGTCCGTGTTATTCAGGc                                                                                                            >  2:147873/1‑147 (MQ=255)                                        cGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑c                                                                                                       <  3:74441/149‑1 (MQ=255)                                        cGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑c                                                                                                       <  8:100546/149‑1 (MQ=255)                                         gTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑ca                                                                                                      >  1:286938/1‑149 (MQ=255)                                             cgcAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATc                                                                                                  <  2:21210/149‑1 (MQ=255)                                                        cAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATg                                                                                       <  4:317815/149‑1 (MQ=255)                                                        cAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATg                                                                                       >  4:237024/1‑149 (MQ=255)                                                         aaaGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGc                                                                                      >  7:226934/1‑149 (MQ=255)                                                           aGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTc                                                                                    <  4:171444/149‑1 (MQ=255)                                                                 gATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCAt                                                                                >  2:58414/1‑147 (MQ=255)                                                                  aTCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAg                                                                             <  1:58414/149‑1 (MQ=255)                                                                  aTCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAg                                                                             >  8:310239/1‑149 (MQ=255)                                                                     aaGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAGTCACATCCAGCGGTGATGCTCCCATCAgg                                                                            <  1:21248/148‑1 (MQ=255)                                                                           gacgacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATACTCCCATCAGGCGAGcgtc                                                                    >  4:170793/1‑149 (MQ=255)                                                                              gacgaAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGcgtccgt                                                                 <  7:219565/149‑1 (MQ=255)                                                                                 gaAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGa                                                                         >  8:37306/1‑138 (MQ=255)                                                                                 gaAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGa                                                                         <  7:37306/138‑1 (MQ=255)                                                                                 gaAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGG‑CAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTctct                                                             <  8:296425/149‑1 (MQ=255)                                                                                     gACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTtcc                                                          <  3:237024/149‑1 (MQ=255)                                                                                                   aTGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTgcg                                            <  7:310239/149‑1 (MQ=255)                                                                                                              aTACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCca                                 >  3:52070/1‑149 (MQ=255)                                                                                                              aTACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCca                                 <  8:226934/149‑1 (MQ=255)                                                                                                                   g‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACagag                            <  1:147873/149‑1 (MQ=255)                                                                                                                         cacaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCaa                                                  <  2:88135/122‑1 (MQ=255)                                                                                                                         cacaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCaa                                                  >  1:88135/1‑122 (MQ=255)                                                                                                                         cacaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctg                       >  5:2405/1‑149 (MQ=255)                                                                                                                          acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgc                      >  1:27122/1‑149 (MQ=255)                                                                                                                          acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctg                       <  6:275912/148‑1 (MQ=255)                                                                                                                          acaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGACctgc                      >  5:139992/1‑149 (MQ=255)                                                                                                                            aGGGAAACAGTCTGGAAAACCCCGCTTCGGCGTGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAT‑CACATCCAGCGGTGATGCTCCCATCAGGCGa                                                                         <  5:194916/96‑1 (MQ=255)                                                                                                                            aGGGAAACAGTCTGGAAAACCCCGCGTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCctgctg                    >  3:43977/1‑149 (MQ=255)                                                                                                                             gggAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTgg                   <  2:286938/149‑1 (MQ=255)                                                                                                                                   cAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg             <  6:27642/149‑1 (MQ=255)                                                                                                                                     gcaaagAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTggg            >  6:25488/6‑148 (MQ=255)                                                                                                                                     gTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGcg                                                                      >  1:290219/1‑90 (MQ=255)                                                                                                                                     gTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGcg                                                                      <  2:290219/90‑1 (MQ=255)                                                                                                                                         ggAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcg       >  3:166738/1‑149 (MQ=255)                                                                                                                                         ggAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTACGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcg       >  1:27895/1‑149 (MQ=255)                                                                                                                                         agAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg             >  3:140745/2‑143 (MQ=255)                                                                                                                                         agAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTgg             <  4:140745/142‑1 (MQ=255)                                                                                                                                          gAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTc                                   <  6:347853/120‑1 (MQ=255)                                                                                                                                          gAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTc                                   >  5:347853/1‑120 (MQ=255)                                                                                                                                          gAAAACCCCGCCGAAGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATcgc      >  6:308058/1‑149 (MQ=255)                                                                                                                                             aaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCAt                                                                                <  6:232416/72‑1 (MQ=255)                                                                                                                                             aaCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGATCCCAt                                                                                >  5:232416/1‑72 (MQ=255)                                                                                                                                              aCCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGACCTGCTGGCGCTGGGCATCGCGaaa  <  6:139992/149‑1 (MQ=255)                                                                                                                                               ccccGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGaaa  >  7:17121/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                              TTCTAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAA  >  NC_002947/450771‑451059

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