Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	8,741	G→A	6.5%	A168A (GCG→GCA)	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	8,741	0	G	A	6.5%	163.6 / 2.1	62	A168A (GCG→GCA)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base G (23/34);  new base A (2/2);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

CTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAG  >  NC_000913/8602‑8882                                                                                                                                            |                                                                                                                                              cTGGTATTAGCAACGAGCTTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGaa                                                                                                                                            <  2:704317/143‑1 (MQ=255)          gCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGTAGAACAGCTGGAAAAAGAAGGCACCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCCCAGGCTCGTGCTTGTGCGGCCGCGGGCgtg                                                                                                                                   >  2:61033/1‑143 (MQ=255)          gCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCAGAAGCGGGCgtg                                                                                                                                   >  1:1084819/1‑143 (MQ=255)          gCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTGCTCCTTCGCTCAGGCGCGTGCTTGTGCGGAAGCGGGcccg                                                                                                                                   >  2:873854/1‑140 (MQ=255)          gCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATAAATGATGCAGAACAGCTGGAAAAAGAAGGAATCAACAGAAACATGACAATGCTGAACTCCATCGCTCAGGCTCGTGCTTGTGCAGAAGCGGGCgtg                                                                                                                                   >  1:984845/1‑143 (MQ=255)                gtCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGATCCTGCTGTTCTCCTTCGCTCAGTCTCCTGCTTGTTCGGACGCGGGAGTAATCCTg                                                                                                                             >  1:987797/2‑143 (MQ=255)                   gTATTCTGATCAAACCGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGAtc                                                                                                                          <  1:674901/143‑1 (MQ=255)                    tATTCTGATCAAACTGGCTTCTCCCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCCGCTGTTCTCCTTCGCTCCGGCTCGTGCTTGTGCGGAGGCGCGCGTGTTCCTGAtcc                                                                                                                         >  2:435942/1‑142 (MQ=255)                      ttCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCg                                                                                                                       <  1:612515/143‑1 (MQ=255)                        ccaaacAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGcc                                                                                                                     <  1:502460/138‑1 (MQ=255)                        cTGATCAAACTGGCTTCTAACTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGcc                                                                                                                     <  2:1035020/143‑1 (MQ=255)                          gagcAAACTGGCTTCTACCTGGCAGGGTATACGAGCTACAAAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGt                                                                                                                   <  1:462287/140‑1 (MQ=255)                          gATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGt                                                                                                                   <  2:980267/143‑1 (MQ=255)                            tcgaaCTGGGTGCTACGTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAACGAAGGCATCAACTGTAGCCTGGCCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGttt                                                                                                                 <  2:802817/140‑1 (MQ=255)                               aaCTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTttgtt                                                                                                              <  1:890643/143‑1 (MQ=255)                                 cTGGCTTCTACCTGGCAGGGTATCCGTGCTGCCGAACAGCTGCAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTgg                                                                                                            <  2:230670/143‑1 (MQ=255)                                  atgcgtgtACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGc                                                                                                           <  2:711096/136‑1 (MQ=255)                                   ggCTTCTCCCCGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAGCCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGcc                                                                                                          <  2:200005/143‑1 (MQ=255)                                      ttCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGGAACCTGACCCTGCTGTTCTCCTTCGCCCAGGCTCGTGCTAGTGCGGAAGCGCGCGTGATCCCGAGCTCGCCGTGTCTTGGCCGTa                                                                                                       >  2:510086/1‑143 (MQ=255)                                           ccTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTACCCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCtt                                                                                                  <  2:749824/143‑1 (MQ=255)                                            cTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTCACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGATTGTGCGGAAGCGGGAGTGGTACTGATCTAGCCGTTTGTTGGTCGTATTCTTg                                                                                                 >  1:527149/1‑143 (MQ=255)                                                  gggTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTGCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGt                                                                                           <  1:442692/143‑1 (MQ=255)                                                   ggTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTa                                                                                          <  1:888159/143‑1 (MQ=255)                                                      aTCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACaa                                                                                       >  1:320219/1‑143 (MQ=255)                                                      aTCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTAGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGCCTGGTAGaa                                                                                       >  1:493248/1‑143 (MQ=255)                                                          gTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCg                                                                                   >  1:619802/1‑143 (MQ=255)                                                           tgctgcAGAACAGCTGGAAAAAGAAGGCATCACCTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGa                                                                                  <  1:954014/143‑1 (MQ=255)                                                            gctgcAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCAGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGac                                                                                 <  2:1084819/143‑2 (MQ=255)                                                            gctgcAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTGCTCCTTCGCTCAGGCTCGTGCTTGTGCAGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGaa                                                                                 <  2:984845/143‑1 (MQ=255)                                                                    aCAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGACTACCGATa                                                                         >  1:865304/1‑143 (MQ=255)                                                                         tGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaagaaa                                                                    >  1:702367/1‑143 (MQ=255)                                                                         tGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGAACTACCTAGTCGATTTGTTTTTTCCTGATCTCCCCCTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaagaaa                                                                    <  1:811538/143‑1 (MQ=255)                                                                         tGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCCTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaagaaa                                                                    <  1:344733/143‑1 (MQ=255)                                                                           gAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaagaaaga                                                                  <  1:438886/143‑1 (MQ=255)                                                                                agaagGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACg                                                             <  1:157946/143‑1 (MQ=255)                                                                                     gCATCAACTGTACCCTGACCCTGCGGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTGCCTGAGCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCg                                                        <  1:362541/143‑1 (MQ=255)                                                                                     gCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGCATTCTTGGCTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCg                                                        >  2:1078560/1‑143 (MQ=255)                                                                                        tCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACCAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCa                                                     >  2:1080049/1‑143 (MQ=255)                                                                                         cAACTGTATCCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGATTGTGGTTAATCGTGTGTGTTTTTTAGCTCTCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCag                                                    <  1:465791/143‑1 (MQ=255)                                                                                              gTAACCTGACCCTGCTGTTCTCCTCCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGAt                                               <  2:320219/143‑1 (MQ=255)                                                                                                 aCCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGCATACCGATAAGAGAGAGTACGCGCCGGCAGACGAGCCg                                            >  2:148388/1‑143 (MQ=255)                                                                                                  ccTGACCCGTCTGGTCTCCTTCGCTCTGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCgg                                           <  2:144101/143‑1 (MQ=255)                                                                                                        ccggctgTTCTCACTAGCTCCGACTCTTGCTAGTGCGGAAGCGGGCTTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTgg                                     <  2:527374/140‑1 (MQ=255)                                                                                                          tgctgTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGGATACCGATAAGAAAGAGTACGCTCCGGCGGACGATCCCGGCGAGGtt                                   >  2:362359/1‑143 (MQ=255)                                                                                                           gctgTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGttt                                  >  1:952348/1‑143 (MQ=255)                                                                                                               ttCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTCACTGGTACGAAGCGAAGACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGt                              >  1:512622/1‑143 (MQ=255)                                                                                                                   ccTTCGCTCAGGCTCGTGCTCGTGCGGACGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCCCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCt                          <  1:648842/143‑1 (MQ=255)                                                                                                                   ccTTCGCTCAGGCTCGGGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCt                          <  1:842902/143‑1 (MQ=255)                                                                                                                    cTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTg                         <  1:602266/143‑1 (MQ=255)                                                                                                                      tCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTTCAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGaa                       <  1:557576/143‑1 (MQ=255)                                                                                                                           gaGGCTCGTGCTTGTGAGGCAGCGGGCGTGTCCGTTCTCTCGCCGTTTGTTGCCCGTATGCTTGACTGGTACAAAGCGAATAGCGATAAGAAAGAGTACGCTCCTGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTa                  <  2:96480/142‑1 (MQ=255)                                                                                                                            aGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTAc                 <  1:61033/143‑1 (MQ=255)                                                                                                                            aGGCTCGTGCGGGTGCGGAAGCGGGCGTGGTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAAGACCGATAAGAAAGAGGACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTAc                 <  2:163908/143‑1 (MQ=255)                                                                                                                              gCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCa               <  1:321604/143‑1 (MQ=255)                                                                                                                               cTCGTGCGTGGGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAg              <  2:261447/143‑1 (MQ=255)                                                                                                                                tcgtgctcgtgcGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAACGCGAATACCGATAAGAAAGAGTACGCTCCAGCAGAAAATCCGGGCGTGGGTTCTGTATCCGAAATATAcaaaa             >  2:27372/1‑139 (MQ=255)                                                                                                                                 cgtgcgtgGGCGGAAGCGGGCGTCTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGta            <  1:166637/143‑1 (MQ=255)                                                                                                                                    gCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGGGGTTTCTGTATCTGAAATCTACCAGtacta         >  1:904030/1‑143 (MQ=255)                                                                                                                                       tgtgCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACaa      >  1:327247/1‑143 (MQ=255)                                                                                                                                        gtgCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATCAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTGCaaa     >  2:216566/1‑143 (MQ=255)                                                                                                                                        gtgCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATACGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACaaa     >  1:567457/1‑143 (MQ=255)                                                                                                                                           cGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACAGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTATGTAGCTGAAATGTACCACCACTACAAagag  >  2:221816/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                              CTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAG  >  NC_000913/8602‑8882

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑