Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,373,236	G→C	9.6%	noncoding (70/82 nt)	PP_mr55 →	P31 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,373,236	0	G	C	9.6%	209.7 / 12.5	73	noncoding (70/82 nt)	PP_mr55	P31 ncRNA
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base C (3/4);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-02

CCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGT  >  NC_002947/5373099‑5373384                                                                                                                                          |                                                                                                                                                     ccGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCAccc                                                                                                                                           >  5:256923/1‑149 (MQ=255)  cGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCTAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCActcc                                                                                                                                          <  4:76593/149‑1 (MQ=255)       gCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCttt                                                                                                                                     <  4:36868/149‑1 (MQ=39)           cGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATAGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCAtg                                                                                                                                 >  4:46981/1‑149 (MQ=255)                  ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttt                                                                                                                                    <  7:67119/139‑2 (MQ=38)                  ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttt                                                                                                                                    >  8:67119/1‑138 (MQ=38)                  ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaa                                                                                                                          >  6:188620/1‑149 (MQ=255)                  ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCGGaa                                                                                                                          >  4:46015/1‑149 (MQ=255)                    cACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaa                                                                                                                          <  5:205607/147‑1 (MQ=255)                      ccGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTACCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGAGGTATGGACACTTCACCCCGCTTTCATGAACTGaattc                                                                                                                       >  8:184291/1‑148 (MQ=255)                       cgcgGAAAAAACTGCCGTCGGCAAGGGCCAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGCAGTGACGTATGGACACGTCACCCCGCTTTCATGTGATGaattcaa                                                                                                                     >  6:268379/1‑149 (MQ=255)                        gcgGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGacac                                                                                                                                                  <  7:280418/119‑1 (MQ=255)                        gcgGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGAGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGacac                                                                                                                                                  >  8:280418/1‑119 (MQ=255)                                tACTGCAGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTgcg                                                                                                            >  2:49449/1‑149 (MQ=255)                                  cTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTgcgcg                                                                                                          <  7:173884/149‑1 (MQ=255)                                   tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCtt                                                                                                                                      >  7:220922/1‑120 (MQ=255)                                   tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCtt                                                                                                                                      <  8:220922/120‑1 (MQ=255)                                      cGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCttt                                                                                                      <  8:70689/149‑1 (MQ=255)                                          ggCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGt                                                                                                  >  5:279192/1‑149 (MQ=255)                                           gCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTg                                                                                                 >  5:242746/1‑149 (MQ=255)                                            cAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCCCGCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGGTTTTAGTg                                                                                                 >  5:25313/1‑148 (MQ=255)                                            cAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGACACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTgg                                                                                                <  2:179339/149‑1 (MQ=255)                                               gggCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTg                                                                                                 <  1:138160/145‑1 (MQ=255)                                               gggCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTg                                                                                                 >  2:138160/1‑145 (MQ=255)                                                ggCGAGTTTGGCTGGCTTTGTCAAGAAGACGTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTAGGGACACGTCAGCCGGCTGTCAGGTGCTGAATTCAATTAGCTGCGCGCTTTTAGTGGacac                                                                                            >  8:129239/1‑149 (MQ=255)                                                   gAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGacacac                                                                                          >  2:308766/1‑148 (MQ=255)                                                     gTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTa                                                                                        >  4:247360/1‑148 (MQ=255)                                                      tATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTaaa                                                                                      <  4:262445/149‑1 (MQ=255)                                                          gctggctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCt                                                                                   <  6:242746/148‑1 (MQ=255)                                                              gctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTa                                                                                <  1:49449/147‑1 (MQ=255)                                                                      aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGa                                                                       >  4:97136/1‑148 (MQ=255)                                                                         aagaCCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAg                                                                                     >  6:10531/1‑131 (MQ=255)                                                                         aagaCCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAg                                                                                     <  5:10531/131‑1 (MQ=255)                                                                           gaCCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGc                                                                              >  3:135611/1‑136 (MQ=255)                                                                           gaCCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGc                                                                              <  4:135611/136‑1 (MQ=255)                                                                             ccTCCCGTAAATCAAGCATGAAAGGGCGATCTCAGGTTTCAAAAAAGCGGAGTGACGTATTGACACTTGACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgtg                                                               <  5:268379/149‑1 (MQ=255)                                                                                cccGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattcaat                                                                                                                    <  8:160203/93‑1 (MQ=255)                                                                                cccGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattcaat                                                                                                                    >  7:160203/1‑93 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgt                                                                <  2:22149/142‑1 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgt                                                                >  1:22149/1‑142 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                  >  3:234037/1‑140 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                  <  4:234037/140‑1 (MQ=255)                                                                                     aaaTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTAc                                                                          >  2:59516/1‑130 (MQ=255)                                                                                     aaaTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTAc                                                                          <  1:59516/130‑1 (MQ=255)                                                                                      aaTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGt                                                       <  3:247360/148‑1 (MQ=255)                                                                                        tCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTa                                                    >  1:278397/1‑149 (MQ=255)                                                                                         cAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACAACGTAAAGATTTACATTAATGAAAAGTGTTTTTGCGTTATAc                                                   >  1:185973/1‑149 (MQ=255)                                                                                            gCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACACCGCTTTCATGTGCTGAATACAATTCGCTGCGCGC‑TTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACcgcg                                               >  1:46945/1‑149 (MQ=255)                                                                                             cATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACcgcg                                               >  5:603/1‑149 (MQ=255)                                                                                                gAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattc                                                                                                                       <  8:269015/74‑1 (MQ=255)                                                                                                gAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattc                                                                                                                       >  7:269015/1‑74 (MQ=255)                                                                                                gAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCACGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGAt                                            >  2:107757/1‑149 (MQ=255)                                                                                                 aaaGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGAt                                            <  7:285996/148‑1 (MQ=255)                                                                                                 aaaGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGAt                                            <  6:279192/148‑1 (MQ=255)                                                                                                  aaGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGAGACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGAtgt                                          >  6:27284/1‑149 (MQ=255)                                                                                                  aaGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGATTACAGAAAAGTGTGTTTTCGGTCTACCGCGGAtgt                                          >  2:78229/1‑149 (MQ=255)                                                                                                    gggCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATgtgt                                        <  7:163047/149‑1 (MQ=255)                                                                                                           cGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTaa                                  >  4:266264/1‑148 (MQ=255)                                                                                                             cAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTaaa                                                                                      <  6:43196/94‑1 (MQ=255)                                                                                                             cAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTaaa                                                                                      >  5:43196/1‑94 (MQ=255)                                                                                                             cAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAac                               <  5:188620/149‑1 (MQ=255)                                                                                                               ggTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAacac                             <  2:278397/149‑1 (MQ=255)                                                                                                                         aaGCGGACTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTCTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCTATGCGa                   <  6:25313/149‑1 (MQ=255)                                                                                                                             ggggTGACGTGTCCATACGTCACTCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGa                                             <  2:97097/116‑1 (MQ=38)                                                                                                                             ggggTGACGTGTCCATACGTCACTCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGa                                             >  1:97097/4‑119 (MQ=38)                                                                                                                             ggggTGACGTGTCCATACGTCACTCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGa               >  8:60690/4‑149 (MQ=255)                                                                                                                                 tGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGa           >  2:216540/1‑149 (MQ=255)                                                                                                                                     gTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCa                         >  6:269831/1‑131 (MQ=255)                                                                                                                                     gTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCa                         <  5:269831/131‑1 (MQ=255)                                                                                                                                       aTGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCt      <  4:27575/148‑1 (MQ=255)                                                                                                                                       aTGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCt      >  3:27575/1‑148 (MQ=255)                                                                                                                                        tGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                  <  4:31591/87‑1 (MQ=255)                                                                                                                                        tGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                  >  3:31591/1‑87 (MQ=255)                                                                                                                                          gACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTgtgt  >  5:158373/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                     CCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGT  >  NC_002947/5373099‑5373384

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