Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,373,238	C→T	12.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	12.0%	149.8 / 9.1	58	noncoding (72/82 nt)	PP_mr55	P31 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base T (2/5);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

CGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGTT  >  NC_002947/5373100‑5373385                                                                                                                                           |                                                                                                                                                     cGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCActcc                                                                                                                                            <  3:226249/149‑1 (MQ=255)     ggCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtt                                                                                                                                        <  7:220505/149‑1 (MQ=39)      gCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCttt                                                                                                                                       >  2:24441/1‑149 (MQ=255)         ccGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGGAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCAt                                                                                                                                    <  6:71911/149‑1 (MQ=255)         ccGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttttt                                                                                                                                    <  1:296160/149‑4 (MQ=38)                 ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttttt                                                                                                                                    >  8:42628/1‑138 (MQ=38)                 ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttttt                                                                                                                                    <  7:42628/141‑4 (MQ=38)                 ttCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaa                                                                                                                            >  7:288909/1‑149 (MQ=255)                   cACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaat                                                                                                                           <  5:13109/148‑1 (MQ=255)                    aCCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaatt                                                                                                                          >  3:176683/1‑148 (MQ=255)                      cgcgGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattc‑aa                                                                                                                      >  1:183486/1‑149 (MQ=255)                      cgcgGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattc‑aa                                                                                                                      >  4:130581/1‑149 (MQ=255)                        cgGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTTaattc‑aatt                                                                                                                    <  1:24441/149‑1 (MQ=255)                          gAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCg                                                                                                                  >  3:152104/1‑149 (MQ=255)                           aaaaTACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGc                                                                                                                 >  1:53891/1‑149 (MQ=255)                                  tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttttt                                                                                                                                    <  4:98532/124‑4 (MQ=37)                                  tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttttt                                                                                                                                    >  3:98532/1‑121 (MQ=37)                                  tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTgcgcgc                                                                                                          >  7:200274/1‑149 (MQ=255)                                         ggCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAg                                                                                                    <  8:288909/148‑1 (MQ=255)                                              gggCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGaca                                                                                              <  2:53891/149‑1 (MQ=255)                                                     tATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTg                                                                                                  <  4:96113/138‑1 (MQ=255)                                                     tATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTg                                                                                                  >  3:96113/1‑138 (MQ=255)                                                          ctggctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTg                                                                                   <  1:235923/149‑1 (MQ=255)                                                          ctggctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGacac                                                                                             >  8:160880/1‑138 (MQ=255)                                                          ctggctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGacac                                                                                             <  7:160880/138‑1 (MQ=255)                                                             gctTTGTCAAGAAGACCTACCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGc                                                                               <  8:205583/149‑1 (MQ=255)                                                                     aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGt                                                                                                   >  4:129798/1‑121 (MQ=255)                                                                     aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGt                                                                                                   <  3:129798/121‑1 (MQ=255)                                                                     aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGaca                                                                                              >  6:64183/1‑126 (MQ=255)                                                                     aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGaca                                                                                              <  5:64183/126‑1 (MQ=255)                                                                     aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGa                                                                        >  5:235908/1‑148 (MQ=255)                                                                       gaagaCCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGaaaa                                                                     <  2:183486/149‑1 (MQ=255)                                                                             cTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGGTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgggtgg                                                               >  8:124960/1‑148 (MQ=255)                                                                             cTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgtgt                                                               >  8:145637/1‑149 (MQ=255)                                                                             cTCCAGTAAATAAAGCATGAAAGGGAGATCGCAGCTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCACGCTTTAATGTGCTCAATTC‑AATTCGCTGCGCGATTTTAGAGGACACACTAAATCTGAATCTTACAGAAAAgtgatt                                                               >  8:150186/1‑146 (MQ=255)                                                                              tCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGtt                                                              >  1:158527/1‑149 (MQ=255)                                                                              tCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGtt                                                              >  3:175902/1‑149 (MQ=255)                                                                                ccGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgtgt                                                               <  4:87311/146‑1 (MQ=255)                                                                                ccGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAgtgtgt                                                               >  3:87311/1‑146 (MQ=255)                                                                                ccGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGtttt                                                            <  5:95900/149‑1 (MQ=255)                                                                                    aaaTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGt                                                        >  5:345867/1‑149 (MQ=255)                                                                                      aTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGaaaa                                                                     <  4:288737/134‑1 (MQ=255)                                                                                      aTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGaaaa                                                                     >  3:288737/1‑134 (MQ=255)                                                                                       tCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTa                                                     >  6:149863/1‑149 (MQ=255)                                                                                         aaGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTAcc                                                   >  5:316885/1‑149 (MQ=255)                                                                                         aaGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTAc                                                    <  8:249640/148‑1 (MQ=255)                                                                                                aaaGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGAt                                             <  2:158527/148‑1 (MQ=255)                                                                                                    ggCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGtc                                        >  2:346717/1‑149 (MQ=255)                                                                                                     gCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGtct                                       <  6:235908/149‑1 (MQ=255)                                                                                                      cctgggCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGtctc                                      <  2:271881/144‑1 (MQ=255)                                                                                                        aTCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTctct                                     >  1:38394/1‑148 (MQ=255)                                                                                                         tCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGGTTTCGGTCTACCGCGGATGTGTGTCTa                                    >  6:266141/1‑148 (MQ=255)                                                                                                            cAGGTTGCAAAAAAGCGTAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTTTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAac                                <  4:176683/149‑1 (MQ=255)                                                                                                               gTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAAcaca                             <  6:345867/149‑1 (MQ=255)                                                                                                                ttGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACAt                            <  3:130581/149‑1 (MQ=255)                                                                                                                   cAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCaa                         >  2:290429/1‑149 (MQ=255)                                                                                                                     aaaaaGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATg                       >  4:230310/1‑149 (MQ=255)                                                                                                                                  aCGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACa          <  4:175902/149‑1 (MQ=255)                                                                                                                                          acacGTCACCCCTCTTTCATTTTCTTAATTC‑AATTCGCTTCTCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTgtgt   <  7:150186/148‑1 (MQ=255)                                                                                                                                          acacGTCACCCCTCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTgtgt   <  7:145637/148‑1 (MQ=255)                                                                                                                                          acacGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGa                    <  1:80600/131‑1 (MQ=255)                                                                                                                                          acacGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGa                    >  2:80600/1‑131 (MQ=255)                                                                                                                                          acacGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGtt  <  1:346717/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     CGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACATCTGTGTT  >  NC_002947/5373100‑5373385

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