Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	537,468	G→A	7.2%	noncoding (64/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,468	0	G	A	7.2%	200.4 / 4.2	69	noncoding (64/66 nt)	PP_mr07	P26 ncRNA
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base A (3/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

GCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGG  >  NC_002947/537322‑537615                                                                                                                                                   |                                                                                                                                                          gCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGc                                                                                                                                                         >  3:255530/1‑149 (MQ=255)     aaGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCttt                                                                                                                                                     >  4:303633/1‑149 (MQ=255)     aaGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGccat                                                                                                                                                      >  2:178694/1‑146 (MQ=255)      aGGCGTTTCGAAAGAAGCCGCTGAAGACGCAAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGccatca                                                                                                                                                    <  1:178694/149‑5 (MQ=255)        gCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCCAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGGTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTc                                                                                                                                                  <  5:314581/149‑1 (MQ=255)         cGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                                 <  7:78202/149‑1 (MQ=255)          gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCATAAGCCACCAGc                                                                                                                                                         <  2:94242/140‑1 (MQ=255)          gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCATAAGCCACCAGc                                                                                                                                                         >  1:94242/1‑140 (MQ=255)                aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAtt                                                                                                                                           <  2:277236/148‑1 (MQ=255)                aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtg                                                                                                                                          <  7:118143/149‑1 (MQ=255)                 aagaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                                         <  4:148973/149‑1 (MQ=255)                 aagaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                                         >  7:259654/1‑149 (MQ=255)                     agCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTGGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTCCCGTTATTGGTGGc                                                                                                                                     <  6:277470/149‑1 (MQ=255)                        cGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCgg                                                                                                                                  >  8:220715/1‑149 (MQ=255)                             aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                                             >  2:82245/1‑149 (MQ=255)                             aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                                             >  6:88554/1‑149 (MQ=255)                             aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATCCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAgg                                                                                                                             >  7:12784/1‑149 (MQ=255)                              aGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGc                                                                                                                                     >  4:308772/1‑140 (MQ=255)                              aGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGc                                                                                                                                     <  3:308772/140‑1 (MQ=255)                              aGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGc                                                                                                                            <  6:50173/149‑1 (MQ=255)                                  gCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGt                                                                                                                        <  3:91438/149‑1 (MQ=255)                                      agaagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTccc                                                                                                                     >  8:202056/1‑148 (MQ=255)                                       gaagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGa                                                                                                                   <  1:216915/149‑1 (MQ=255)                                        aagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGAt                                                                                                                  <  5:88554/149‑1 (MQ=255)                                          gaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                                >  2:304148/1‑149 (MQ=255)                                            agCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCTAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAAc                                                                                                               <  6:57052/148‑1 (MQ=255)                                                     aagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGc                                                                                                        >  8:86234/1‑146 (MQ=255)                                                        cAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGaca                                                                                                  >  1:321689/1‑149 (MQ=255)                                                               tAACGTCGTGCTGACTTATTCTCGACGTTGCGTTGCCTGCCTGAGCGTTGCGCGAACGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTTTTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCaa                                                                                                      <  7:211650/138‑1 (MQ=255)                                                               tAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                                           >  3:127503/1‑149 (MQ=255)                                                               tAAAGTCGAGCTGAAGTAATCTCGACCATGCGAAGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTAAGCCACCGGCCATTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCCAAAACGACCGGCAAGacaacggcag                                                                                           >  8:211650/1‑142 (MQ=255)                                                                 aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAt                                                                                                                                            <  1:60898/98‑1 (MQ=255)                                                                 aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                                           <  1:90966/147‑1 (MQ=255)                                                                 aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgg                                                                                         <  1:304148/149‑1 (MQ=255)                                                                 aaGTCGAGCTGAAGTAATCTCGACCTTGCGATGCAAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAt                                                                                                                                            >  2:60898/1‑98 (MQ=255)                                                                    tCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgggtg                                                                                      >  1:91888/1‑149 (MQ=255)                                                                          tGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCa                                                                               >  3:97326/1‑148 (MQ=255)                                                                            aaGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑c                                                                          <  7:220715/149‑1 (MQ=255)                                                                                 aTCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAg                                                                      <  4:127503/148‑1 (MQ=255)                                                                                  tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAggg                                                                    >  7:11968/1‑149 (MQ=255)                                                                                      gACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑ttt                                                              <  4:255530/149‑1 (MQ=255)                                                                                        ccTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGc                                                            <  3:303633/149‑1 (MQ=255)                                                                                            gCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGa                                                        >  6:13830/1‑149 (MQ=255)                                                                                              gATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGAtt                                                      >  6:313157/1‑149 (MQ=255)                                                                                              gATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGAtt                                                      >  6:326136/1‑149 (MQ=255)                                                                                                  ccAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCt                                                  >  7:66377/1‑149 (MQ=255)                                                                                                  ccAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCAACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGTTCCCGATAACGAGCTGCaa                                                                                                      >  3:152024/1‑103 (MQ=255)                                                                                                   cAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATCCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTg                                                 <  8:12784/149‑1 (MQ=255)                                                                                                      ccTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGgct                                              >  1:244032/1‑149 (MQ=255)                                                                                                        tGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGgctgc                                            >  2:177957/1‑149 (MQ=255)                                                                                                        tGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGgctg                                             >  5:105909/1‑148 (MQ=255)                                                                                                         gAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGgct                                              <  6:133916/146‑1 (MQ=255)                                                                                                         gAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGgct                                              >  5:133916/1‑146 (MQ=255)                                                                                                              ttGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTccctc                                      >  5:296372/1‑147 (MQ=255)                                                                                                                 aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCgg                                   <  8:259654/149‑1 (MQ=255)                                                                                                                 aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCgg                                   <  5:313157/149‑1 (MQ=255)                                                                                                                   cGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTTCTCCcggcgg                                   >  6:20084/1‑147 (MQ=255)                                                                                                                         gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGcc                           <  1:82245/149‑1 (MQ=255)                                                                                                                                gccggtggCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg                    >  6:155587/4‑149 (MQ=255)                                                                                                                                gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg                    >  3:144915/1‑149 (MQ=255)                                                                                                                                 ccggtggCATAAGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg                    >  8:155671/3‑148 (MQ=255)                                                                                                                                 cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGc                   <  7:202056/149‑1 (MQ=255)                                                                                                                                         ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCaaa                                                                             <  5:161768/87‑1 (MQ=255)                                                                                                                                         ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCaaa                                                                             >  6:161768/1‑87 (MQ=255)                                                                                                                                         ttATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAc            >  7:237460/1‑148 (MQ=255)                                                                                                                                            tGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAg        >  3:241135/1‑149 (MQ=255)                                                                                                                                             gCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGc       >  6:10931/1‑149 (MQ=255)                                                                                                                                               cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTg     <  2:321689/149‑1 (MQ=255)                                                                                                                                                  cGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCCTGTAGAGGCCAAGGGTGAGTTGCACGAATTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGAAGGTGACCaa         >  1:106682/1‑148 (MQ=255)                                                                                                                                                   ggCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgggg  >  2:265462/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                          GCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGC‑‑GCAAA‑‑CCAAGGG‑G‑TTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGG  >  NC_002947/537322‑537615

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