Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,925,932 (T)7→6 coding (493/657 nt) yobB → C‑N hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,925,9260T.98.8% 367.6 / ‑2.5 84coding (487/657 nt)yobBC‑N hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (47/36);  total (48/36)
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.02e-01

CAACCGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCC  >  NC_000913/1925827‑1926021
                                                                                                   |                                                                                               
caacCGCAAGGTATGGATATGGACCCGACCTGTTCCCTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAt                                                                                                 >  1:78893/1‑100 (MQ=255)
  acCGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttt                                                                                               >  2:878838/1‑100 (MQ=255)
  acCGAAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTTTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttt                                                                                               <  2:40010/100‑1 (MQ=255)
    cGCAAGGTATGGTTATGGACCCGAGCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGAGCTGCTGGCGTCGGCCCGTCGCTTACAAttttt                                                                                             >  2:573475/1‑100 (MQ=255)
    cGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttt                                                                                             <  1:448958/100‑1 (MQ=255)
     gCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAtttttt                                                                                            >  1:526578/1‑100 (MQ=255)
     gCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAtttttt                                                                                            >  2:2825703/1‑100 (MQ=255)
      cAAGGTATGGATATGGCCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttc                                                                                           <  1:561457/100‑2 (MQ=255)
      cAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttc                                                                                           <  2:1595784/100‑2 (MQ=255)
      cAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttc                                                                                           <  1:439464/100‑2 (MQ=255)
        aGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttcac                                                                                         <  1:878839/100‑4 (MQ=255)
        aGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttcac                                                                                         <  2:1798442/100‑4 (MQ=255)
             tggatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAg                                                                                   >  1:2189697/1‑100 (MQ=255)
             tggatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAg                                                                                   <  1:1195145/100‑1 (MQ=255)
               gatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTa                                                                                 >  1:2407444/1‑100 (MQ=255)
                atatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTAc                                                                                <  2:1393953/100‑1 (MQ=255)
                 tatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACt                                                                               <  2:2715838/100‑1 (MQ=255)
                     gaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGAt                                                                           <  1:1611953/100‑1 (MQ=255)
                     gaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGAt                                                                           >  2:1812813/1‑100 (MQ=255)
                       cccGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTg                                                                         >  2:1469778/1‑100 (MQ=255)
                       cccGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTg                                                                         <  1:1096135/100‑1 (MQ=255)
                          gACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCg                                                                      >  2:2570791/1‑100 (MQ=255)
                           aCCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTACTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCg                                                                      <  2:1474227/99‑1 (MQ=255)
                            ccTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGCTTGCCGt                                                                     >  2:2457177/1‑99 (MQ=255)
                            ccTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTg                                                                    <  2:2189684/100‑1 (MQ=255)
                            ccTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTg                                                                    >  1:1168609/1‑100 (MQ=255)
                             cTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGc                                                                   >  1:2771126/1‑100 (MQ=255)
                              tGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCt                                                                  >  2:1107286/1‑100 (MQ=255)
                               gTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTg                                                                 >  1:150097/1‑100 (MQ=255)
                               gTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTg                                                                 >  1:516571/1‑100 (MQ=255)
                                   aCTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATgg                                                             <  2:1073407/100‑1 (MQ=255)
                                   aCTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATgg                                                             <  1:1902304/100‑1 (MQ=255)
                                    cTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGc                                                            >  1:1846141/1‑100 (MQ=255)
                                    cTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGc                                                            >  1:1319783/1‑100 (MQ=255)
                                     tGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCa                                                           >  1:2488346/1‑100 (MQ=255)
                                      gTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTAGAA‑TTTTTTCAGAGCAGTACTCGATTGCCGTGCTGATGGCaa                                                          >  2:2408973/1‑100 (MQ=255)
                                      gTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCaa                                                          >  2:422970/1‑100 (MQ=255)
                                       tttACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCaaa                                                         <  1:677812/100‑1 (MQ=255)
                                         tactacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATg                                                       >  1:2725021/1‑100 (MQ=255)
                                         tactacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATg                                                       <  1:103132/100‑1 (MQ=255)
                                          actacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGc                                                      >  1:1450606/1‑100 (MQ=255)
                                          actacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGc                                                      <  2:903514/100‑1 (MQ=255)
                                           ctacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGcc                                                     >  1:979140/1‑100 (MQ=255)
                                             acAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGTTGGCAAATGCCCg                                                   >  1:705294/1‑100 (MQ=255)
                                             acAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCg                                                   >  1:2031580/1‑100 (MQ=255)
                                              cAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGt                                                  >  1:187301/1‑100 (MQ=255)
                                              cAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGt                                                  <  1:1120788/100‑1 (MQ=255)
                                               aGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTg                                                 <  1:592861/100‑1 (MQ=255)
                                                gggCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCTGTgg                                                <  1:769499/100‑1 (MQ=255)
                                                 ggCACTGTCTGGGGGAACCCGACCTTCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGt                                               <  2:1235060/100‑1 (MQ=255)
                                                 ggCAATGTCTGGGAGACCCCGGCCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTGACATCAGTACTCGGTTGCCGCGTTGATGGCAAATGCCCGTGGt                                               <  2:1235049/100‑1 (MQ=255)
                                                 ggCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGt                                               <  2:98078/100‑1 (MQ=255)
                                                 ggCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGt                                               <  1:1639760/100‑1 (MQ=255)
                                                    aaTGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAAc                                            <  2:1573459/100‑1 (MQ=255)
                                                      tGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAg                                          >  2:362069/1‑100 (MQ=255)
                                                       gTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGt                                         <  1:320781/100‑1 (MQ=255)
                                                        tCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTg                                        >  2:2259140/1‑100 (MQ=255)
                                                        tCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTg                                        >  1:1214662/1‑100 (MQ=255)
                                                         cTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGt                                         >  2:1215063/1‑98 (MQ=255)
                                                         cTGGGAGAACCCGACCTGCTGGAGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGc                                       >  2:1683866/1‑100 (MQ=255)
                                                          tGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCa                                      >  1:1510185/1‑100 (MQ=255)
                                                              agaACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGGTGGCAACTGCCCGTGGTAACAGTGCACTAt                                  >  1:409623/1‑100 (MQ=255)
                                                                aaCCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATgg                                <  1:2835020/100‑1 (MQ=255)
                                                                 aCCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATggg                               >  1:2815331/1‑100 (MQ=255)
                                                                 aCCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATggg                               >  2:1386743/1‑100 (MQ=255)
                                                                 aCCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATggg                               >  1:1779820/1‑100 (MQ=255)
                                                                  cccGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATg                                 >  2:2645208/1‑97 (MQ=255)
                                                                   ccGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGAt                             <  2:1886420/100‑1 (MQ=255)
                                                                     gACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGa                           <  2:328598/100‑1 (MQ=255)
                                                                     gACCTGCGGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACGTCAGTGCTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGa                           <  1:1254793/100‑1 (MQ=255)
                                                                       ccTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAat                         >  1:37880/1‑100 (MQ=255)
                                                                        ctgctgGCGTCGGCCCGTCGCTTACAA‑TTTTTGCACATCCGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAata                        <  2:2031566/100‑1 (MQ=255)
                                                                         tgctgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTTCACTATGGGATGAatat                       >  2:2881480/1‑100 (MQ=255)
                                                                         tgctgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAatat                       >  1:1110870/1‑100 (MQ=255)
                                                                         tgctgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAatat                       <  1:573477/100‑1 (MQ=255)
                                                                          gctgGCGTCGGCCCGTCGCTTACAATTTTTTT‑TTTTTTTTTTTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAatat                       >  2:2071054/1‑99 (MQ=255)
                                                                           ctgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGTTAACAGTGCACTATGGGATGAatattg                     >  2:715041/1‑98 (MQ=255)
                                                                           ctgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATgg                     <  1:1174977/100‑1 (MQ=255)
                                                                           ctgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATgg                     >  1:1633061/1‑100 (MQ=255)
                                                                            tgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGgt                    <  1:1881483/100‑1 (MQ=255)
                                                                            tgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGgt                    <  1:2602182/100‑1 (MQ=255)
                                                                               cGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAAAATGgtcgt                 >  1:1045668/1‑100 (MQ=255)
                                                                                  cGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGtct               <  1:1072819/99‑1 (MQ=255)
                                                                                    gCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCAt            >  1:1525596/1‑100 (MQ=255)
                                                                                    gCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCAt            >  1:1917121/1‑100 (MQ=255)
                                                                                    gCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCAt            >  1:114430/1‑100 (MQ=255)
                                                                                      ccGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCg          >  2:2817360/1‑100 (MQ=255)
                                                                                       cGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGt         <  1:1562138/100‑1 (MQ=255)
                                                                                        gTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGtt        >  2:2218942/1‑100 (MQ=255)
                                                                                           gCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTcgc     >  1:549731/1‑100 (MQ=255)
                                                                                            cTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTcgcg    <  1:2457193/100‑1 (MQ=255)
                                                                                             ttACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCgcgc   >  1:149890/1‑100 (MQ=255)
                                                                                             ttACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCgcgc   <  1:1129610/100‑1 (MQ=255)
                                                                                              tACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGcc  <  1:1466111/100‑1 (MQ=255)
                                                                                              tACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGcc  <  1:1934706/100‑1 (MQ=255)
                                                                                              tACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGcc  >  1:356788/1‑100 (MQ=255)
                                                                                                   |                                                                                               
CAACCGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCC  >  NC_000913/1925827‑1926021

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