Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,102,340	(C)6→5	coding (513/720 nt)	trmB ←	tRNA m(7)G46 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,102,340	0	C	.	97.4%	170.5 / ‑2.5	39	coding (513/720 nt)	trmB	tRNA m(7)G46 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (17/21);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTTACGCAGACGACGGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGC  >  NC_000913/3102101‑3102557                                                                                                                                                                                                                                                |                                                                                                                                                                                                                          tttACGAAGACGACGGCTACGGTCCTATGCCATTATTTCACCCTCTCGAACATTAAGTCCCATCCTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCCGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTg                                                                                                                                                                                                               >  2:53781/1‑251 (MQ=255)        aGACGACGGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCttttt                                                                                                                                                                                                         >  1:235739/1‑250 (MQ=255)        aGACGACGGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATGTCGTCTCCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACGTCACTTCAAGCATATGTTCCGCATAGGGTTCCCAGTCGGTCGCCGTATGGAATGCGCCCGC‑AGCTGCATTTTGCTTTTTa                                                                                                                                                                                                        <  2:37223/251‑1 (MQ=255)            gacgGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTAcc                                                                                                                                                                                                      >  1:151613/1‑249 (MQ=255)                                       cACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATAc                                                                                                                                                                          <  2:151613/250‑1 (MQ=255)                                        aCCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACgg                                                                                                                                                                        >  1:200162/1‑251 (MQ=255)                                             ctcGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGttt                                                                                                                                                                   >  1:76875/1‑251 (MQ=255)                                                      tAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgcg                                                                                                                                                          >  1:170337/1‑251 (MQ=255)                                                               tACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTtgtg                                                                                                                                                 >  2:101916/1‑251 (MQ=255)                                                                               cGATGACCACGTTGTTCAAATTTCGTCACCGTACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACTGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGTGaaa                                                                                                                                 <  1:100736/251‑1 (MQ=255)                                                                               cGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaa                                                                                                                                  >  1:224009/1‑250 (MQ=255)                                                                                gATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATTCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATCCGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCTTTTATTATGGCGCGCGTTGTGCCACGGGTCAGGGaaaa                                                                                                                                <  2:45781/251‑1 (MQ=255)                                                                                       aCGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTg                                                                                                                         <  2:224009/251‑1 (MQ=255)                                                                                           tgttCAGATTTCGTCACCGGACGTGATGCCGGTCGCGGTACGTAATCATTGCGCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGAATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGTGTCAGGGAAAAAGATCTGCAc                                                                                                                      <  2:197223/250‑1 (MQ=255)                                                                                              tcgaaattCGTCTCCGGTCGTGAGGCCTGACGCGTTACGTAGTCATTGCGCTCTGACAGGTTTTTATAACCGTGAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGACGGCATATGGAATACGCCCCC‑ATCTGCAGTCGGCTTTTGACCAGGTCGGCAAGCGGCTCCTGAACGATGCGGCGTTGATTAGGGCGCGATTGGGGCCTTGGCTCAGGGAAAAGGAGCTGCGCCAt                                                                                                                   <  2:45480/244‑1 (MQ=255)                                                                                              tCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATAC‑TCCCCCAGCTGCATTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTTCCACGGGTCAGGGAAAAAGAGCTGCACCAt                                                                                                                   <  1:350165/250‑1 (MQ=255)                                                                                                aaaTTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATAATTGATCTCTGACAGGTTTTTATAACAGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGCACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATgcg                                                                                                                >  1:358681/1‑251 (MQ=255)                                                                                                         cACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCATTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTAGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAAt                                                                                                        <  1:244325/250‑1 (MQ=255)                                                                                                                    aTGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCa                                                                                            <  1:145186/251‑1 (MQ=255)                                                                                                                                gTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATgcagca                                                                                <  1:183455/251‑1 (MQ=255)                                                                                                                                      aaTCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGTCGCTCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATTAATTGTCAGTAATCATTTTATGCAGCACTTCa                                                                           >  2:203877/1‑250 (MQ=255)                                                                                                                                         cATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAAcc                                                                        <  1:235708/250‑1 (MQ=255)                                                                                                                                          aTTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACGTCAACCGc                                                                      <  1:139811/251‑1 (MQ=255)                                                                                                                                                         gTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCAc                                                       <  1:243419/251‑1 (MQ=255)                                                                                                                                                           ttttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCAcg                                                      <  1:335289/250‑1 (MQ=255)                                                                                                                                                                 aaCCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACgcgcggg                                                <  2:150910/250‑4 (MQ=255)                                                                                                                                                                           agaagaCATCACTGCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAAcc                                     <  1:161256/251‑1 (MQ=255)                                                                                                                                                                              agaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACcttc                                  <  1:15665/251‑1 (MQ=255)                                                                                                                                                                               gaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTAGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTAACAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGAAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACcttct                                 >  2:215819/1‑251 (MQ=255)                                                                                                                                                                                aCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACcttctt                                >  1:90214/1‑251 (MQ=255)                                                                                                                                                                                         cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATCCGCCCCC‑AGCTGAAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAACAGAACTGCACCATGCGCAATGCATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTTCTTAAACCTTCTTCATgcgc                         >  2:10436/1‑249 (MQ=255)                                                                                                                                                                                          aaGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCAt                                                                                                                   <  2:281164/158‑1 (MQ=255)                                                                                                                                                                                          aaGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCAt                                                                                                                   >  1:281164/1‑158 (MQ=255)                                                                                                                                                                                           aGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGaa                     <  2:76875/251‑1 (MQ=255)                                                                                                                                                                                            gCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAg                    <  2:200162/251‑1 (MQ=255)                                                                                                                                                                                                       gCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGgcacgc          >  2:79657/1‑250 (MQ=255)                                                                                                                                                                                                        cATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGgcacgcac        <  1:101916/251‑1 (MQ=255)                                                                                                                                                                                                         aTAAGGTTCCCCGTCGGTCGCCATATGGAATACGACCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAACAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCAcc       >  2:344975/1‑251 (MQ=255)                                                                                                                                                                                                              gTTCCCAGTCGGTCGCCATATGGAATACGCCCCC‑AGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGc  >  1:149186/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                          TTTACGCAGACGACGGCTACGGTTCTTTGCCATTATTTCACCCTCTCGAACATTAAGTCCCATACTCCGTGACCAAGACGATGACCACGTTGTTCAAATTTCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGC  >  NC_000913/3102101‑3102557

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