Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	450,921	G→A	9.2%	noncoding (118/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,921	0	G	A	9.2%	178.2 / 7.1	65	noncoding (118/127 nt)	PP_mr05	RsmY ncRNA
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base A (4/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

TAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGC  >  NC_002947/450774‑451067                                                                                                                                                     |                                                                                                                                                    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                                                          >  1:88135/1‑122 (MQ=255)                                                                                                                      cacaGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCaa                                                          <  2:88135/122‑1 (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)                                                                                                                                                    cggcggGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGc  <  4:168848/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TAAGATTGGATCCACTGTTAATCTTTACCACATGGACGTTGCGCAGGAAGCGCAAAGGACAGATCAAGGACGACGACGAAGGACCACCTGCCAGGATGGCGGGGCGATACGG‑ATGTCACAGGGAAACAGTCTGGAAAACCCCGCTTCGGCGGGGTTTTCTTTGCGTGCGTGTTTTTCAGGCCAG‑CACATCCAGCGGTGATGCTCCCATCAGGCGAGCGTCCGTCTCTTCCTCCAGCAACGTGCGCAGCAGCTCCACAGAGGCCTGCTGGCGCTGGGCATCGCGAAATACCAGGC  >  NC_002947/450774‑451067

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