Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,373,238	C→T	10.0%	noncoding (72/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,238	0	C	T	10.0%	199.2 / 9.3	72	noncoding (72/82 nt)	PP_mr55	P31 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base T (3/4);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGTTG  >  NC_002947/5373099‑5373386                                                                                                                                            |                                                                                                                                                     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                                                                                                                                      >  8:67119/1‑138 (MQ=38)                  ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttt                                                                                                                                      <  7:67119/139‑2 (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                                                                                <  4:135611/136‑1 (MQ=255)                                                                           gaCCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGc                                                                                >  3:135611/1‑136 (MQ=255)                                                                             ccTCCCGTAAATCAAGCATGAAAGGGCGATCTCAGGTTTCAAAAAAGCGGAGTGACGTATTGACACTTGACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgtg                                                                 <  5:268379/149‑1 (MQ=255)                                                                                cccGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattcaat                                                                                                                      <  8:160203/93‑1 (MQ=255)                                                                                cccGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattcaat                                                                                                                      >  7:160203/1‑93 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgt                                                                  >  1:22149/1‑142 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgt                                                                  <  2:22149/142‑1 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                    >  3:234037/1‑140 (MQ=255)                                                                                   gTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                    <  4:234037/140‑1 (MQ=255)                                                                                     aaaTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTAc                                                                            <  1:59516/130‑1 (MQ=255)                                                                                     aaaTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTAc                                                                            >  2:59516/1‑130 (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                           <  5:269831/131‑1 (MQ=255)                                                                                                                                     gTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCa                           >  6:269831/1‑131 (MQ=255)                                                                                                                                       aTGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCt        <  4:27575/148‑1 (MQ=255)                                                                                                                                       aTGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCt        >  3:27575/1‑148 (MQ=255)                                                                                                                                        tGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                    >  3:31591/1‑87 (MQ=255)                                                                                                                                        tGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgt                                                                    <  4:31591/87‑1 (MQ=255)                                                                                                                                          gACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTgtgt    >  5:158373/1‑149 (MQ=255)                                                                                                                                           acacGTCACCCGGATTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGtt   <  7:282401/149‑1 (MQ=255)                                                                                                                                           acacGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGAGCGCTTTTAGTCGACACACTAAAGCTGTAGCTTACAGAAAACTGTGTTTTCGGTCTACCGCCGATGTGTCTCTAAAACACATCAATGCGACAGAc                >  3:48705/1‑136 (MQ=255)                                                                                                                                            cacGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGTTg  <  1:107757/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     CCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGTTG  >  NC_002947/5373099‑5373386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑