Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,768 T→C 11.2% intergenic (+34/+14) yeaD → / ← yeaE putative aldose 1‑epimerase YeaD/methylglyoxal reductase YeaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,864,7680TC11.2% 146.3 / 8.7 65intergenic (+34/+14)yeaD/yeaEputative aldose 1‑epimerase YeaD/methylglyoxal reductase YeaE
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (3/4);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-02

AAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGG  >  NC_000913/1864630‑1864906
                                                                                                                                          |                                                                                                                                          
acaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAACGTGCCCAATGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                        <  2:562775/141‑1 (MQ=255)
aaaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCgg                                                                                                                                        <  1:58478/143‑1 (MQ=255)
aaaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATCTACGTTAATGTTGTGTGCCGGGGGCAATGCATCCgg                                                                                                                                        >  2:65035/1‑143 (MQ=255)
aaaaCATTTGTTTGTGGAGAAACGGCCTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGATGCGAAACGTTAATTTACTTTAATGTTGTGTGCCGGGTGCAATGCATCCgg                                                                                                                                        <  2:491048/143‑1 (MQ=255)
aaaaCATTTGGTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAGCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCgg                                                                                                                                        <  2:149315/143‑1 (MQ=255)
  aaCATTTGTTTGTGTAGAAACGGCTGACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGca                                                                                                                                      <  2:176971/143‑1 (MQ=255)
    cATTTGTGTGTGTAGACACGTGTGACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcaca                                                                                                                                    <  2:85819/143‑1 (MQ=255)
     aTTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGCAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacac                                                                                                                                   >  1:661786/1‑143 (MQ=255)
       ttgtttgtGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaa                                                                                                                                 <  1:647872/143‑1 (MQ=255)
         gtttgtGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTTTGCCGGGCGCAAAGCATCCGACACacaaca                                                                                                                               >  2:809600/1‑143 (MQ=255)
              tgtAGATACCGCTGACGCCTGAGTCATGCAAAAAATAACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCCACCGGCACAAAACAAcaca                                                                                                                          >  2:838472/1‑143 (MQ=255)
              tgtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACACAACAttaac                                                                                                                          <  1:881065/143‑5 (MQ=38)
              ggtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATcaca                                                                                                                          <  1:162650/142‑1 (MQ=255)
               gtAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCAAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCGACCGCCACACAACATCacac                                                                                                                         >  1:469873/1‑143 (MQ=255)
                   caaCGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTAGCGCCATCCACTCGCGGTCCGAAACGTTAATTTACCTTAAGGTTCTGCGGCGGGTTCAATGCCTGCGGCACAACACATGCCACCcta                                                                                                                     >  2:900891/2‑143 (MQ=255)
                     aCGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGACACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATc                                                                                                                   <  1:841628/143‑1 (MQ=255)
                       gggttACGCTTCAGAAACGGAAAAAGTGACCAAAGAGCAACCTGCATATTTGGCGCAAGTCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGACGGGTGCAATGCATCCGGCACACAACATCACACCATATCCa                                                                                                                 <  2:823266/140‑1 (MQ=255)
                         cTTACGCTTCAGAAACGCAAAAAGTGACCAAGGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgc                                                                                                               <  2:720987/143‑1 (MQ=255)
                          ttACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                                                              >  2:470204/1‑143 (MQ=255)
                          ttACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                                                              >  1:452659/1‑143 (MQ=255)
                          ttACGCGTCCGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGTCATGCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                                                              <  2:415027/143‑1 (MQ=255)
                          ttACGCATCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                                                              <  1:22919/143‑1 (MQ=255)
                                ttCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTCATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGgtt                                                                                                        >  2:301480/1‑143 (MQ=255)
                                ttCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGAAATCCATTCGCGTTGCGAAACGTAAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCCGttc                                                                                                        >  1:767979/1‑142 (MQ=255)
                                       aCGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTACCGGGTGCAATGCAACCGGCACACAACATCACACCATAACCACCGCAGTCTCTCCtgc                                                                                                 >  2:782101/1‑141 (MQ=255)
                                         ttaaaaaTTGACCGAGTAGAAACCTGCACATCTGGCGCAATCCTTTCGCGTTTCTAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                                                               <  2:935628/141‑1 (MQ=255)
                                         gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                                                               <  1:456148/143‑1 (MQ=255)
                                           aaaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCAAACACAATCTCGCCTGGTCCAGCGCCGTTTTTCCttttcct                                                                                             >  1:730466/1‑140 (MQ=255)
                                               aGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCg                                                                                         <  1:560522/143‑1 (MQ=255)
                                                 tGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                                                                       >  2:27250/1‑143 (MQ=255)
                                                    ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGatat                                                                                    <  1:658903/143‑1 (MQ=255)
                                                      aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGatattc                                                                                  >  1:204666/1‑141 (MQ=255)
                                                            aaaCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATc                                                                            >  1:180167/1‑143 (MQ=255)
                                                            aaaCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTCCTGGTGCCGGAACTGccccatc                                                                            >  2:209557/1‑138 (MQ=255)
                                                                 tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGCATTTCCTTATAGTGCCGGATATGCCTTATCCAGCa                                                                       >  2:169557/1‑143 (MQ=255)
                                                                 tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAACACATTCAGCGTAGTTTTTCCTTTTGGAGCCGGATATGCCTTATCCAGCa                                                                       >  1:128410/1‑143 (MQ=255)
                                                                  gCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGGTTTTCCTTTTGGTGCCGGCTATGCCTTATCCAGCAt                                                                      >  1:40213/1‑143 (MQ=255)
                                                                    acaTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                        >  2:298535/1‑75 (MQ=37)
                                                                    acaTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                        <  1:298535/75‑1 (MQ=37)
                                                                    acaGCTGGCGCAATCCATTCGCGTTGCGGACCGTTAATCTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAg                                                                    <  1:390615/143‑1 (MQ=255)
                                                                           gcgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCCTCGCTCAtgc                                                             >  1:157403/1‑141 (MQ=255)
                                                                             gcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATCGCTAATTTCg                                                           >  1:797914/1‑143 (MQ=255)
                                                                             gcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCg                                                           <  2:109853/143‑1 (MQ=255)
                                                                               aaTCCATTCGCGTTGCAAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt                                                         <  1:808751/143‑1 (MQ=255)
                                                                                  gcactCGCGGTGCGCAGCGTTAGGCTACGTCAACGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCCTCTCCGGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa                                                      <  2:357034/139‑1 (MQ=255)
                                                                                  ccgtgggcgTTGCGAAAGGTTCATGTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCAGACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa                                                      <  2:361985/137‑1 (MQ=255)
                                                                                  ccaggcgcgTTGCGAAACGTTAATTTGCGTTAATGTTGTGTGCGGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa                                                      <  2:6681/138‑1 (MQ=255)
                                                                                  ccATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa                                                      <  1:628476/143‑1 (MQ=255)
                                                                                      tCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCAGCCGGCACACAACATCACACCTTATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGaaa                                                  <  2:58131/143‑1 (MQ=255)
                                                                                        gcgTTGCGAAACGTTACAGTACTTTTATGTTGTGATGCGGGTGCAATGCATCCGGCACACCACATGACACGCTATCCAGCGCAGTTTTTCATTCTGGTCCCGGATCAGTCTTATTCAGCAGCGCCAATTGCGCTGAAGAAAGt                                                >  2:876572/1‑143 (MQ=255)
                                                                                          gTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTc                                              <  1:599649/143‑1 (MQ=255)
                                                                                               gAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTCCTTCCGCCGCAGAACGTTCCACCt                                         >  1:744862/1‑143 (MQ=255)
                                                                                                       aaTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGcaca                                 <  1:337221/143‑1 (MQ=255)
                                                                                                              gTTAATGTTGTGTGCCGGGTGCAATGCCTCCGGCACACAACATCACACCATATCCGGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCGAGCATAGCTCATGCCGCTGAAGAACCTGCAACCCCAAACACATCCGcgc                          >  2:52570/1‑141 (MQ=255)
                                                                                                                 aaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCATACTATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCAtttt                       <  2:744865/143‑1 (MQ=255)
                                                                                                                     ttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTCGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCGGAAGAAAGTTCGCCCTCCAGCCCACCCGCATCttgtcg                   >  2:290367/1‑141 (MQ=255)
                                                                                                                      tgtgtgCCGGGTGCAATGCATCCGGCACACCACATCACACCATATCCAGCGCAGTTTCTCCTTTTGGTCCCGAATCTGGCTGAGGCAGCACAGCTAATTCCGCTGAAGAACGTTCGACCTCCACCGCCGCCGCCTCTTGTCgg                  >  2:510543/1‑143 (MQ=255)
                                                                                                                      tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATTGCTCATTCCGCTGAAGAAAGTGCGATCTTAAACACAGCCGCATTTTGTTgg                  >  2:217527/1‑143 (MQ=255)
                                                                                                                            cccggTGCACTGCATCCGGCACACAACATCACACCATCTCCAGCGCAGGTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGATATgg            <  2:382344/140‑1 (MQ=255)
                                                                                                                            ccGGGTGCAATGCATCCAGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCGCAGCCGCATTTTGTTGGACATgg            >  1:23251/1‑143 (MQ=255)
                                                                                                                            ccGGATGCATTGCACCCGGCACACAACATTACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACATCAAGCACAGCCGCATTTTGTTGGACATgg            >  1:720654/1‑143 (MQ=255)
                                                                                                                            ccGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATgg            >  1:487890/1‑143 (MQ=255)
                                                                                                                             ctgctGCAATGCACCCGGCACACAACAGCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATggg           <  2:433842/139‑1 (MQ=255)
                                                                                                                                   cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACATCAAGCACAGCCGCATTTTGTTAGACATGGGCAAtag     >  2:25279/1‑141 (MQ=255)
                                                                                                                                      tGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGTCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTgg  <  1:470204/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
AAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGG  >  NC_000913/1864630‑1864906

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