Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	537,458	T→A	6.3%	noncoding (54/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,458	0	T	A	6.3%	185.6 / 2.6	63	noncoding (54/66 nt)	PP_mr07	P26 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base A (2/2);  total (30/33)
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.48e-01

CTCAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAG  >  NC_002947/537311‑537600                                                                                                                                                     |                                                                                                                                               ctcAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGAAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCt                                                                                                                                                >  2:353842/1‑148 (MQ=255)      gTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGCCTGATGGCTGGTGGCTTATGc                                                                                                                                           >  3:95484/1‑148 (MQ=255)         gTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCAc                                                                                                                                        <  5:360633/148‑1 (MQ=255)             cggAAGGGGTTTCGAACTAAGCCGCTGAAGACGCTAAGATCACGCTGGAAGAAGCAGGCGCTAAAGGCGAGCTTAAGTAATCTCTACCTTGCGATGCCAGCCT‑GAGCGTTGAGCTAAAGGCTGATGGCCGGTGGCATAAGCCACCAGc                                                                                                                                    <  8:140854/146‑1 (MQ=255)             cTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGcc                                                                                                                                   <  7:228586/149‑1 (MQ=255)             cTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCATAAGCCAc                                                                                                                                        <  4:309763/144‑1 (MQ=255)             cTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCATAAGCCAc                                                                                                                                        >  3:309763/1‑144 (MQ=255)             cTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGc                                                                                                                                    >  7:140854/1‑148 (MQ=255)                      tttCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGt                                                                                                                          >  5:144793/1‑149 (MQ=255)                        tCGAAAGAAGCCGCTTAAGACGCTATGAAGAAGCTGTAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGtt                                                                                                                         <  6:224790/148‑1 (MQ=255)                          gAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGTCACCGGCCTTTTTCCGTTAtt                                                                                                                      <  5:306226/149‑1 (MQ=255)                          gAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAtt                                                                                                                      >  6:185878/1‑149 (MQ=255)                          gAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAtt                                                                                                                      >  2:208544/1‑149 (MQ=255)                          gAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAtt                                                                                                                      <  2:22709/149‑1 (MQ=255)                                  ccGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCg                                                                                                              <  6:349799/149‑1 (MQ=255)                                    gCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGt                                                                                                            >  1:85052/1‑149 (MQ=255)                                        aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                        >  1:132798/1‑149 (MQ=255)                                         aGACGCTAAGAAGAAGCTGCAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                        <  1:353842/148‑1 (MQ=255)                                             gCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGt                                                                                                   >  4:330166/1‑149 (MQ=255)                                               tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTcc                                                                                                 >  6:68844/1‑149 (MQ=255)                                                aagaagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTccc                                                                                                >  5:324072/1‑149 (MQ=255)                                                 agaagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCg                                                                                               <  3:330166/149‑1 (MQ=255)                                                        gCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAAc                                                                                          >  8:244798/1‑147 (MQ=255)                                                        gCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGa                                                                                        <  6:324072/149‑1 (MQ=255)                                                             aagaagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGa                                                                                              <  7:308070/138‑1 (MQ=255)                                                             aagaagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGa                                                                                              >  8:308070/1‑138 (MQ=255)                                                               gaagCAGGCGCTAAAGGCGAGCTGAAGTAATCTCGACCTTGCGATCCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCaa                                                                                 <  7:114446/149‑1 (MQ=255)                                                                aagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGc                                                                                   >  2:130338/1‑146 (MQ=255)                                                                    aGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGGGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGacac                                                                            <  2:132798/149‑1 (MQ=255)                                                                       cgcTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCAGCCAGaaaatcg                                                                         >  6:149929/1‑143 (MQ=255)                                                                          tAAAGTCGATCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGATGATGGATGGTGGCTTATGCCACAGGCCTTATTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGATGCAAGACACCCAtgt                                                                      >  5:210829/1‑149 (MQ=255)                                                                          tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtg                                                                       >  1:69356/1‑149 (MQ=255)                                                                          tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                      >  1:293084/1‑149 (MQ=255)                                                                               tCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgggtg                                                                 <  2:85052/149‑1 (MQ=255)                                                                                 gAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCATGACACCCATGTGGGTGgc                                                               >  8:142807/1‑149 (MQ=255)                                                                                 gAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCACatgggtgtc                                                               <  4:196085/149‑3 (MQ=255)                                                                                   gCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGgcg                                                              >  2:133982/1‑148 (MQ=255)                                                                                          tAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCaa                                                                                 <  6:104920/122‑1 (MQ=255)                                                                                          tAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCaa                                                                                 >  5:104920/1‑122 (MQ=255)                                                                                           aaTCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCCGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                     >  6:167815/1‑149 (MQ=255)                                                                                               tcGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGt                                                 <  5:3225/149‑1 (MQ=255)                                                                                                 gACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGttt                                               <  1:133982/149‑1 (MQ=255)                                                                                                    cTTGCTATGCCATCCG‑GAGCGTTGAGCGAATGGCTGTTGTCTGGTGGCTAATGCCACCGGCCTTTTTCCGTTATTGGGGGCCGGTCAGGCCGGTGCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAATCCAAGGGGTTTg                                              <  8:326734/147‑1 (MQ=255)                                                                                                       gCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGa                                         >  2:174887/1‑149 (MQ=255)                                                                                                              cAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTg                                  <  5:185878/149‑1 (MQ=255)                                                                                                               aGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAATGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTgg                                 >  3:173642/1‑149 (MQ=255)                                                                                                               aGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTgg                                 <  7:244798/149‑1 (MQ=255)                                                                                                                 ccT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGgc                                                               <  2:200358/117‑1 (MQ=255)                                                                                                                 ccT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGgc                                                               >  1:200358/1‑117 (MQ=255)                                                                                                                  cT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgct                               >  5:43812/1‑148 (MQ=255)                                                                                                                      aGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTg                                              <  5:188605/130‑1 (MQ=255)                                                                                                                      aGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTg                                              >  6:188605/1‑130 (MQ=255)                                                                                                                       gCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTcc                          >  3:214463/1‑149 (MQ=255)                                                                                                                        cGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAATGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTccc                         <  4:173642/149‑1 (MQ=255)                                                                                                                         gTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCccg                        >  8:195902/1‑149 (MQ=255)                                                                                                                          ttGAGCGAAAGGCCGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCccgc                       <  5:167815/149‑1 (MQ=255)                                                                                                                             aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                    <  7:195902/149‑1 (MQ=255)                                                                                                                             aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                    <  2:293084/149‑1 (MQ=255)                                                                                                                               cGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGga                  <  1:208544/149‑1 (MQ=255)                                                                                                                                  aaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa               >  3:255903/1‑149 (MQ=255)                                                                                                                                  aaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa               <  4:255903/149‑1 (MQ=255)                                                                                                                                      cTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGTCa           <  1:174887/149‑1 (MQ=255)                                                                                                                                             cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc    <  7:142807/149‑1 (MQ=255)                                                                                                                                                gtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGTGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg  <  4:95484/148‑1 (MQ=255)                                                                                                                                                     ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa               >  8:145697/1‑130 (MQ=255)                                                                                                                                                     ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa               <  7:145697/130‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               CTCAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCT‑GAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAG  >  NC_002947/537311‑537600

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