Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,373,233	A→G	11.3%	noncoding (67/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,233	0	A	G	11.3%	159.4 / 10.3	62	noncoding (67/82 nt)	PP_mr55	P31 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base G (2/5);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

ACCAACGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACA  >  NC_002947/5373086‑5373377                                                                                                                                                    |                                                                                                                                                  aCCAACGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTAt                                                                                                                                                  <  7:63903/149‑1 (MQ=255)   cAACGGTCTGACCGTCGGCACCGATGGATTAACCGCGGAAAATACTGCCGTCGGCAAGTGCCAGTATGGCTGGCTATGTAAAGAACACCTAACGTAAATAAAGCATGAAAGGGCGATCGCAGCTTGCAAAAAAGCGTAGTGACGTATg                                                                                                                                                 >  6:13109/1‑148 (MQ=255)    aaCGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGa                                                                                                                                               >  8:108786/1‑149 (MQ=255)    aaCGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGa                                                                                                                                               <  7:108786/149‑1 (MQ=255)               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                                                                                                                            >  3:98532/1‑121 (MQ=37)                                                tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGGGTGACGTGTCCATACGTCACTCCGCtttttt                                                                                                                            <  4:98532/124‑4 (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                                                                                          >  3:96113/1‑138 (MQ=255)                                                                   tATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTg                                                                                          <  4:96113/138‑1 (MQ=255)                                                                        ctggctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTCAAATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTg                                                                           <  1:235923/149‑1 (MQ=255)                                                                        ctggctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGacac                                                                                     <  7:160880/138‑1 (MQ=255)                                                                        ctggctTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGacac                                                                                     >  8:160880/1‑138 (MQ=255)                                                                           gctTTGTCAAGAAGACCTACCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGc                                                                       <  8:205583/149‑1 (MQ=255)                                                                                   aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGt                                                                                           <  3:129798/121‑1 (MQ=255)                                                                                   aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGt                                                                                           >  4:129798/1‑121 (MQ=255)                                                                                   aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGaca                                                                                      <  5:64183/126‑1 (MQ=255)                                                                                   aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGaca                                                                                      >  6:64183/1‑126 (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                                                             >  3:288737/1‑134 (MQ=255)                                                                                                    aTCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGaaaa                                                             <  4:288737/134‑1 (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)                                                                                                                                                    |                                                                                                                                                  ACCAACGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATGCGACAGACTGACA  >  NC_002947/5373086‑5373377

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