Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,891 T→G 8.2% intergenic (+38/‑111) garD → / → prlF galactarate dehydratase/antitoxin PrlF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,276,8910TG 8.2% 140.2 / 3.3 63intergenic (+38/‑111)garD/prlFgalactarate dehydratase/antitoxin PrlF
Reads supporting (aligned to +/- strand):  ref base T (21/33);  new base G (2/4);  total (26/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 = 8.34e-01

TCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACT  >  NC_000913/3276750‑3277033
                                                                                                                                              |                                                                                                                                               
ccggtcctATTCGCGACGTCCGCAGCGGGACGAAGACACCCTTCTCGGATCAAGGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTg                                                                                                                                                <  2:334685/136‑1 (MQ=255)
     tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACGGGTGACACGATTTCTATTGGAACGTCCCTCGCAAAACATGGCCTCAGCGCCATg                                                                                                                                           >  2:149280/1‑143 (MQ=255)
        aTTCTCCACGTCGCCAGCGGGAAGAAGAAATCCTTCTCGGATCAATGGGTGCTGCATAACCAGCTTGTAGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCCTAAGGCCATGttt                                                                                                                                        <  1:468309/143‑1 (MQ=255)
         tgctcGGCGTCGCCAGCGGGAAGGTGAAGACCTGCTCGGATCAATGGGGGGTGCATAACCAGCTGGCGGTGTGTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGtttt                                                                                                                                       <  2:242361/141‑1 (MQ=255)
         gcctcGACGACGCCGGGGGGAAGAAGAAAACCTTCGCGGCTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAGGCCATGtttt                                                                                                                                       <  2:96417/141‑1 (MQ=255)
          tctcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGAACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGttttt                                                                                                                                      >  1:461114/1‑143 (MQ=255)
           ctcGACGTCGCCAGCGGGAAGCAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTa                                                                                                                                     <  1:532333/143‑1 (MQ=255)
             gtgcggcGCCCGCGGGCAGTACAACACCTGCTCGGATCCATGGGCGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTAtt                                                                                                                                   <  2:534334/137‑1 (MQ=255)
                   cGCCAGCGGGAAGACTACACCCTTCTCGCATCACGGGGGGCTGCATAACCAGCTGGCGGTGTTCAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTaa                                                                                                                             <  2:551867/143‑1 (MQ=255)
                    gCCAGCGGGAAGAAGAAGACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                            <  2:490183/143‑1 (MQ=255)
                     ccAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTTGAACGTCGCTCGCAAAACATGGCGTTAGTGCGCGGTTTTTATTGTATAAAg                                                                                                                           >  1:193477/1‑143 (MQ=255)
                                  gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAt                                                                                                              <  1:10957/143‑1 (MQ=255)
                                         ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTaataat                                                                                                      >  1:881846/1‑143 (MQ=255)
                                               gATCAATGGGGCCTCCATAACCAGCTGGCGGTCTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                                <  2:526274/143‑1 (MQ=255)
                                                  cAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgca                                                                                             <  1:597618/143‑1 (MQ=255)
                                                    aTGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgcatt                                                                                           >  1:890591/1‑143 (MQ=255)
                                                      gggggCTGCCTACCCAGCCGGCGGTGTGCGCCCCGGCGCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                         <  1:710668/143‑1 (MQ=255)
                                                       ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAcc                                                                                        >  1:885141/1‑143 (MQ=255)
                                                       ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAcc                                                                                        >  2:805499/1‑143 (MQ=255)
                                                           cTGCATAGCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACAGGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAg                                                                                    <  1:13503/143‑1 (MQ=255)
                                                             gCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGtt                                                                                  <  1:644395/143‑1 (MQ=255)
                                                              cATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCCTTACCTTAGTTc                                                                                 >  1:258437/1‑143 (MQ=255)
                                                               aTAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCa                                                                                >  1:307374/1‑143 (MQ=255)
                                                                    cAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCtt                                                                           <  2:34613/143‑1 (MQ=255)
                                                                    cAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCAGTACCTTAGTTCAAGCtt                                                                           >  2:220396/1‑143 (MQ=255)
                                                                    cAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGAATTCGATTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCAAGTTTTTATTGTGTAAAGCCCCCACGTCCATTAATAATGCATTTGCGTTACATTAGTGCATGCtt                                                                           >  2:232956/1‑143 (MQ=255)
                                                                     aGCTGGCGGTGTTTAGCCCGGGACCGGTGACCGTATTTCTTTTGGAACGTCCCTCGCAAACCATGGCCTTCGTGCCATGTTTTTATTGTTTACAGCCCCCCCGTCCATTAATAATGCATTTGCATTACCTTAGTTCAGGCTta                                                                          <  2:608882/143‑1 (MQ=255)
                                                                           cGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATCATGCATTTGCATTACCTTAGTTCAAGCTTATAAttt                                                                    <  1:593002/143‑1 (MQ=255)
                                                                            ggTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCGCGCTTATCCTTTg                                                                   >  1:229707/1‑143 (MQ=255)
                                                                            ggTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGATTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGATCAAGCTTATAAttta                                                                   >  2:895534/1‑142 (MQ=255)
                                                                              tgtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAg                                                                 <  1:670152/143‑1 (MQ=255)
                                                                               gtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                                <  1:11220/143‑1 (MQ=255)
                                                                                tttAACCCGGCACTGGTGACCTGATTTCTTTTGGAGCGTCCCTCGCAAAACAGGGCCTTAGTGCCATGTTTTTATTGTTTAAAGTCCCCACGTCCATTAATAAGGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCa                                                               <  2:351421/143‑1 (MQ=255)
                                                                                tttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCa                                                               >  1:370474/1‑143 (MQ=255)
                                                                                  tAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGa                                                             <  2:658787/143‑1 (MQ=255)
                                                                                     cccGGCCCCGGTGCCCTGATTTCTTTTGGAACGTCCCTCGCAAAACCTGGGCTTAGGGCCATGTTTTTATTGTTTCAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaaa                                                          <  2:226986/143‑1 (MQ=255)
                                                                                     cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTGAGTGCCATGTTTGTATTGTTTAAAGCCCCCACGTCCATTAATACTGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGCGCAGaaaa                                                          >  2:552027/1‑143 (MQ=255)
                                                                                      ccGGCACCGGTGACCTGATCTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                         <  1:775551/143‑1 (MQ=255)
                                                                                          cACCGGTGACCTGATTTCTTTGGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGc                                                     <  2:380544/143‑1 (MQ=255)
                                                                                             cGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTCAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTGCAAGCTTATAATTTGAGCAGAAAACAGGCTga                                                  >  2:459407/1‑142 (MQ=255)
                                                                                                 gACCTGATTTCTTTTGGAACGTCCCTCGCAAAGCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCACCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                              <  1:185206/143‑1 (MQ=255)
                                                                                                 gACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCAAGTTTATACTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                              <  1:767847/143‑1 (MQ=255)
                                                                                                  aCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAg                                             <  1:359825/143‑1 (MQ=255)
                                                                                                  aCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAg                                             <  1:767144/143‑1 (MQ=255)
                                                                                                  aCCGGAGTTCTTGGGGAACGTCCCGGGCAGAACAGGGCCGTAGTGCCATGTGTTTATTGTTGAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGGTCAAGCTTAGAATTTGAGCAGAAAACAGGCTGTAAAAg                                             <  2:202229/143‑1 (MQ=255)
                                                                                                    cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGGGCCATGTTTTTATTGTTTAAAGCCCGCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAGAAGGa                                           >  1:309474/1‑143 (MQ=255)
                                                                                                      gATCTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                         <  1:149280/143‑1 (MQ=255)
                                                                                                         ttCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCACGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGACTTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTg                                      <  1:869944/143‑1 (MQ=255)
                                                                                                          tCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAACGCCCCCCCGTCCATTCATAATCAATTTGCATTACCTTAGTTCCAGCTTATAATTTGTGACCAAAACACGCTGTAAAAGGACAGTGa                                     >  2:673886/1‑143 (MQ=255)
                                                                                                           cTTTTGGAACGTCCCTCGTATAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGCAAACAGGCTGTAAAAGGACAGTGaa                                    <  2:706455/143‑1 (MQ=255)
                                                                                                           cTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATACTGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                    >  2:653106/1‑143 (MQ=255)
                                                                                                            ttttGGAACGTCCCTCGCAAAACATGGACTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAAt                                   <  1:128541/143‑1 (MQ=255)
                                                                                                                    cGTCCCTCGCACAACATGACCTTAGTGCCGTGTTTTTATTGTTTTAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCg                           <  2:119199/143‑1 (MQ=255)
                                                                                                                            gCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTATTTCAAGCTTATAATTTGAGCAGAACACAGGCTGTCCAAGGCCAGTGAATCCTGCCAGATAATgca                   >  1:419368/1‑142 (MQ=255)
                                                                                                                            gCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGACTACAAGTCCATTAATAATGAATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGAAAGTGAATCATGCACGATAATgct                   >  1:872964/1‑143 (MQ=255)
                                                                                                                              aaaaCATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcg                 >  1:725590/1‑143 (MQ=255)
                                                                                                                              aaaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcg                 >  2:594737/1‑143 (MQ=255)
                                                                                                                               agaCATCGCCTTAGTGACATGTTTTTATTGTTTAAAGCCCCAAACAAATAACTAAAGGCATTTGAATGACCTTATTGCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgc                <  2:193477/141‑1 (MQ=255)
                                                                                                                               aaaCATGGCACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCAAGACCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgcacgc                >  2:673481/1‑143 (MQ=255)
                                                                                                                               aaaCATGGCACTAAGGCCATGTTTTTATAGTTTAAAGCCCACACGTCGATTAATAATGCAATTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGTAAACAGGCTGAGAAAGGACAGTGAATTATGACCGCTAATgatggc                >  1:356361/1‑138 (MQ=255)
                                                                                                                                           tAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCAATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGt     >  1:578677/1‑143 (MQ=255)
                                                                                                                                           tAGTGCCATGTTTTTATTGTTTAAAGCCCCCAAGTCCGTAAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGt     <  1:621530/143‑1 (MQ=255)
                                                                                                                                            aGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTa    <  1:805499/143‑1 (MQ=255)
                                                                                                                                              tGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCACAAAACAGGCTGTAAACGGACAGTCAATCATGCCCCCTACTGCTCGCTCTCACgatctaat  >  2:586291/1‑136 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCCTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACT  >  NC_000913/3276750‑3277033

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: