Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	537,461	T→A	8.0%	noncoding (57/66 nt)	PP_mr07 →	P26 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	537,461	0	T	A	8.0%	218.7 / 4.6	74	noncoding (57/66 nt)	PP_mr07	P26 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base A (3/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAA  >  NC_002947/537314‑537608                                                                                                                                                    |                                                                                                                                                       aGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAATCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTAt                                                                                                                                                        <  7:42751/148‑1 (MQ=255)        aGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCgg                                                                                                                                                >  2:299153/1‑149 (MQ=255)         gCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGc                                                                                                                                               >  8:14679/1‑149 (MQ=255)            gAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCAACGGCCtt                                                                                                                                            >  7:182915/1‑149 (MQ=255)                  gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑g                                                                                                                                     >  3:291607/1‑149 (MQ=255)                    ttCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                       >  6:301841/1‑146 (MQ=255)                    ttCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                       <  5:301841/146‑1 (MQ=255)                       taaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAATAAGCAGGCGATAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCAcca                                                                                                                                                 <  4:189734/132‑2 (MQ=255)                       gAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TAtt                                                                                                                               <  7:115775/149‑1 (MQ=255)                       gAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TAtt                                                                                                                               <  6:116126/149‑1 (MQ=255)                          agaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATtggt                                                                                                                            >  8:119935/1‑149 (MQ=255)                          agaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATtggt                                                                                                                            <  7:119935/149‑1 (MQ=255)                           gaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATtggtg                                                                                                                           >  7:173106/1‑149 (MQ=255)                             agCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATCCCACCGGCCTTTTTCC‑GT‑TATTGGTGGc                                                                                                                         <  1:174417/149‑1 (MQ=255)                                 gCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGt                                                                                                                     >  7:234809/1‑149 (MQ=255)                                  cTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGGTATATTTGTGGCCgg                                                                                                                      >  1:259682/1‑149 (MQ=255)                                     aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAg‑g                                                                                                                >  8:60745/1‑149 (MQ=255)                                     aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAg                                                                                                                  <  4:291607/148‑1 (MQ=255)                                       gACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑Gcc                                                                                                              >  8:179899/1‑149 (MQ=255)                                              agaagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTccc                                                                                                        >  2:63422/1‑148 (MQ=255)                                                aagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGa                                                                                                      >  2:194188/1‑148 (MQ=255)                                                  gaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATaa                                                                                                   >  8:113510/1‑149 (MQ=255)                                                  gaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCAccag                                                                                                                                                >  6:149328/1‑105 (MQ=255)                                                  gaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCAccag                                                                                                                                                <  5:149328/107‑3 (MQ=255)                                                      cTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGTTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAg                                                                                               >  3:176241/1‑149 (MQ=255)                                                      cTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAg                                                                                               <  1:132916/149‑1 (MQ=255)                                                         gaagaaGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGc                                                                                              <  5:168451/147‑1 (MQ=255)                                                           agaagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCa                                                                                          <  7:179899/149‑1 (MQ=255)                                                            gaagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCaa                                                                                         <  5:203096/149‑1 (MQ=255)                                                              agCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGa                                                                                       <  1:194188/149‑1 (MQ=255)                                                              agCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGGTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGa                                                                                       <  3:216648/149‑1 (MQ=255)                                                                cAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                       >  8:238145/1‑102 (MQ=255)                                                                cAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                       <  7:238145/102‑1 (MQ=255)                                                                    cgcTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCa                                                                                 >  5:100025/1‑149 (MQ=255)                                                                       tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                              >  6:6485/1‑149 (MQ=255)                                                                         aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                              <  1:299153/147‑1 (MQ=255)                                                                         aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgg                                                                            <  3:132934/149‑1 (MQ=255)                                                                              gAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGgc                                                                       <  8:106196/149‑1 (MQ=255)                                                                               aGCTGAAGTAATCTCGACCTTGCGATGCCAGCCCGACCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCACatgggtgtct                                                                      <  5:45879/149‑4 (MQ=255)                                                                                    aaGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TAtt                                                                                                                               >  5:135668/1‑88 (MQ=255)                                                                                    aaGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TAtt                                                                                                                               <  6:135668/88‑1 (MQ=255)                                                                                    aaGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAGTGCCGGTCCCGATAACGAGCTGCAAGACACCCACatgggtgtcttg                                                                    <  1:63422/147‑6 (MQ=255)                                                                                         aTCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAgg                                                            >  5:195860/1‑149 (MQ=255)                                                                                         aTCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                             <  8:173106/148‑1 (MQ=255)                                                                                             cGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCaa                                                                   >  2:76605/1‑138 (MQ=255)                                                                                             cGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCaa                                                                   <  1:76605/138‑1 (MQ=255)                                                                                              gACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGttt                                                       >  1:199782/1‑149 (MQ=255)                                                                                                 cTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTTATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCGGATAACGAGCTGCAAGACACGCATGTGGGTGGCGCAAACCAAGGGGTTTGCa                                                    <  7:60745/149‑1 (MQ=255)                                                                                                 cTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAACCCAAGGGGTTTGCa                                                    >  6:58364/1‑149 (MQ=255)                                                                                                 cTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCa                                                    >  4:288711/1‑149 (MQ=255)                                                                                                 cTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCa                                                    <  2:199782/149‑1 (MQ=255)                                                                                                    gCGATGCCAGCCTGAGCGTTGAGCGAAAGGGTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCAc                                                   <  4:176241/147‑1 (MQ=255)                                                                                                     cGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGTCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGCCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGa                                                 <  5:81803/148‑1 (MQ=255)                                                                                                         gCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTc                                            >  5:35827/1‑149 (MQ=255)                                                                                                            aGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCTCCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTgg                                         >  3:275350/1‑149 (MQ=255)                                                                                                            aGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTgg                                         <  7:119319/149‑1 (MQ=255)                                                                                                             gCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgc                                        <  7:14679/149‑1 (MQ=255)                                                                                                                     gTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCccg                                >  1:166281/1‑149 (MQ=255)                                                                                                                         aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCcggcgg                            >  6:108502/1‑149 (MQ=255)                                                                                                                         aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                            <  5:6485/149‑1 (MQ=255)                                                                                                                          gCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                            >  8:48620/1‑148 (MQ=255)                                                                                                                           cGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                            >  6:107056/1‑147 (MQ=255)                                                                                                                            gAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTAC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                            >  4:285499/1‑146 (MQ=255)                                                                                                                              aaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                       <  2:249731/149‑1 (MQ=255)                                                                                                                              aaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                       <  8:234809/149‑1 (MQ=255)                                                                                                                              aaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                       <  6:195860/149‑1 (MQ=255)                                                                                                                                     aTGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGCCag                  <  5:108502/147‑2 (MQ=255)                                                                                                                                       ggCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAAc                >  4:254877/1‑147 (MQ=255)                                                                                                                                       ggCCGGTGGCATAAGCCACCAGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCaaa                 <  1:271579/146‑1 (MQ=255)                                                                                                                                       ggCCGGTGGCATAAGCCACCAGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCaaa                 >  2:271579/1‑146 (MQ=255)                                                                                                                                        gccggtggCATAAGCCACCAGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATCTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGCAGAAGCCAAACaac             >  6:101969/4‑148 (MQ=255)                                                                                                                                            gtggCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg          <  5:107056/148‑1 (MQ=255)                                                                                                                                                cTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAc     <  6:35827/149‑1 (MQ=255)                                                                                                                                                 ttATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCC‑GGAGAAGCCAAACAAGCAGGTGAcc    >  5:112169/1‑148 (MQ=255)                                                                                                                                                   aTGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAcc    <  3:288711/147‑1 (MQ=255)                                                                                                                                                    tGCCACCGGCCTTTTTGC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAATGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCaa  <  8:182915/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       AGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCC‑GT‑TATTGGTGGCCGGTCAG‑GCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAA  >  NC_002947/537314‑537608

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