Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 889,923 A→G intergenic (‑104/‑166) ybjL ← / → ybjM putative transport protein YbjL/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913889,9230AG100.0% 251.9 / NA 68intergenic (‑104/‑166)ybjL/ybjMputative transport protein YbjL/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/33);  total (35/33)

TCCCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGT  >  NC_000913/889786‑890058
                                                                                                                                         |                                                                                                                                       
tCCCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCg                                                                                                                                       >  2:521134/1‑140 (MQ=255)
 cccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGa                                                                                                                                      <  1:414931/140‑1 (MQ=255)
  ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                                     >  2:130619/1‑140 (MQ=255)
  ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                                     >  1:372992/1‑140 (MQ=255)
    aTTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATag                                                                                                                                   >  2:1171081/1‑140 (MQ=255)
          cAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagt                                                                                                                             <  1:521134/140‑1 (MQ=255)
            aTTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagtaa                                                                                                                           >  1:526855/1‑140 (MQ=255)
               cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                        <  1:654101/140‑1 (MQ=255)
                  cGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAatat                                                                                                                     >  1:1250744/1‑140 (MQ=255)
                         tataTTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTAc                                                                                                              <  1:1171081/140‑1 (MQ=255)
                             ttCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgt                                                                                                          <  2:372992/140‑1 (MQ=255)
                               cACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgtgt                                                                                                        >  1:476205/1‑140 (MQ=255)
                                aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                       >  2:988258/1‑140 (MQ=255)
                                aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                       <  2:1068116/140‑1 (MQ=255)
                                  ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                                     >  1:804774/1‑140 (MQ=255)
                                  ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                                     >  2:320610/1‑140 (MQ=255)
                                  ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                                     >  2:1025853/1‑140 (MQ=255)
                                     ggCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcac                                                                                                  >  1:318412/1‑140 (MQ=255)
                                     ggCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcac                                                                                                  >  1:155318/1‑140 (MQ=255)
                                      gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAGGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                                 >  2:1206834/1‑140 (MQ=255)
                                      gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                                 <  2:1250744/140‑1 (MQ=255)
                                      gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                                 >  2:1167549/1‑140 (MQ=255)
                                       cTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACAt                                                                                                <  2:318412/140‑1 (MQ=255)
                                              cttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTAt                                                                                         >  1:50868/1‑140 (MQ=255)
                                               ttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATc                                                                                        <  1:320610/140‑1 (MQ=255)
                                                tGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCa                                                                                       >  1:891005/1‑140 (MQ=255)
                                                 gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAGGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                      <  1:1206834/140‑1 (MQ=255)
                                                 gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                      <  1:130619/140‑1 (MQ=255)
                                                 gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                      <  2:334463/140‑1 (MQ=255)
                                                    tACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACgg                                                                                   >  1:491480/1‑140 (MQ=255)
                                                      cTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCt                                                                                 >  1:361474/1‑140 (MQ=255)
                                                                 ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                                      >  2:296749/1‑140 (MQ=255)
                                                                 ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                                      >  1:145373/1‑140 (MQ=255)
                                                                  tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                     <  1:988258/140‑1 (MQ=255)
                                                                  tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                     <  1:779929/140‑1 (MQ=255)
                                                                  tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                     <  2:476205/140‑1 (MQ=255)
                                                                  tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                     <  2:891005/140‑1 (MQ=255)
                                                                  tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                     <  2:155318/140‑1 (MQ=255)
                                                                  tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                     <  2:170587/140‑1 (MQ=255)
                                                                   gaaagaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                                    <  2:50868/140‑1 (MQ=255)
                                                                     aagaaaTGGTAATTTACGTTAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTg                                                                  <  1:302017/140‑1 (MQ=255)
                                                                      agaaaTGGTGATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                                 >  2:453903/1‑140 (MQ=255)
                                                                      agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                                 >  2:1297870/1‑140 (MQ=255)
                                                                       gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                                <  2:491480/140‑1 (MQ=255)
                                                                       gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                                <  2:878838/140‑1 (MQ=255)
                                                                       gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                                <  1:1167549/140‑1 (MQ=255)
                                                                       gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                                <  2:231923/140‑1 (MQ=255)
                                                                           tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtat                                                            >  2:477697/1‑140 (MQ=255)
                                                                           tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtat                                                            >  2:551255/1‑140 (MQ=255)
                                                                           tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtat                                                            >  2:255372/1‑140 (MQ=255)
                                                                                 tttACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                                <  1:633757/130‑1 (MQ=255)
                                                                                 tttACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                                >  2:633757/1‑130 (MQ=255)
                                                                                        tAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAAt                                               >  2:407900/1‑140 (MQ=255)
                                                                                           tGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCg                                            >  1:1140808/1‑140 (MQ=255)
                                                                                           tGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCg                                            >  1:623049/1‑140 (MQ=255)
                                                                                             ttttttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                                          <  1:453903/140‑1 (MQ=255)
                                                                                             ttttttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                                          <  1:1297870/140‑1 (MQ=255)
                                                                                                 ttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCgg                                      <  1:551255/140‑1 (MQ=255)
                                                                                                  tACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATAGCATGGAATGCGAAGCGGCt                                    <  1:889429/140‑1 (MQ=255)
                                                                                                             aGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACt                          <  2:361474/140‑1 (MQ=255)
                                                                                                                       ttcattcaGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCcaca                <  2:1210480/140‑1 (MQ=255)
                                                                                                                                cgcACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTg                         >  2:320361/1‑122 (MQ=255)
                                                                                                                                cgcACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTg                         <  1:320361/122‑1 (MQ=255)
                                                                                                                                  cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCca                  >  2:997159/1‑127 (MQ=255)
                                                                                                                                  cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCca                  <  1:997159/127‑1 (MQ=255)
                                                                                                                                  cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGaaaaa     >  2:664361/1‑140 (MQ=255)
                                                                                                                                  cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGaaaaa     >  2:1000856/1‑140 (MQ=255)
                                                                                                                                     cTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAAtgt  >  1:228157/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TCCCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGT  >  NC_000913/889786‑890058

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: