Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,373,225	A→G	10.8%	noncoding (59/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,225	0	A	G	10.8%	167.9 / 10.4	65	noncoding (59/82 nt)	PP_mr55	P31 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base G (2/5);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

ACCTCACCAACGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATG  >  NC_002947/5373081‑5373364                                                                                                                                                 |                                                                                                                                             accTCACCAACGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTCGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGa                                                                                                                                          >  7:249640/1‑149 (MQ=255)  ccTCACCAACGGTCTGACCGCCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGAc                                                                                                                                         >  2:235923/1‑149 (MQ=255)      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                                                                                                 >  4:130581/1‑149 (MQ=255)                                         cgcgGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGaattc‑aa                                                                                                 >  1:183486/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)                                                     tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGa                                                                                                                                          <  6:138791/97‑1 (MQ=255)                                                     tGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGa                                                                                                                                          >  5:138791/1‑97 (MQ=255)                                                     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                                                                              <  3:129798/121‑1 (MQ=255)                                                                                        aagaagACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGt                                                                              >  4:129798/1‑121 (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                                                >  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)                                                                                                                                                 |                                                                                                                                             ACCTCACCAACGGTCTGACCGTCGGCACCGATGGCTTCACCGCGGAAAATACTGCCGTCGGCAAGGGCGAGTATGGCTGGCTTTGTCAAGAAGACCTCCCGTAAATCAAGCATGAAAGGGCGATCGCAGGTTGCAAAAAAGCGGAGTGACGTATGGACACGTCACCCCGCTTTCATGTGCTGAATTC‑AATTCGCTGCGCGCTTTTAGTGGACACACTAAAGCTGTAGCTTACAGAAAAGTGTGTTTTCGGTCTACCGCGGATGTGTCTCTAAAACACATCAATG  >  NC_002947/5373081‑5373364

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