Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,494,883:1 (A)7→8 100% coding (261/564 nt) idnK → D‑gluconate kinase, thermosensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,494,8761.A93.9% 496.5 / 21.6 147coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (4/5);  new base A (80/58);  total (84/63)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCG  >  NC_000913/4494777‑4494966
                                                                                                    |                                                                                               
tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAtt                                                                                                  <  1:1621519/100‑1 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aa                                                                                            >  1:1393658/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaa                                                                                           <  2:3616762/100‑1 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaa                                                                                           >  1:1367908/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaa                                                                                           >  1:455420/1‑100 (MQ=255)
      atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑a                                                                                        <  2:2606573/100‑1 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGCAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa                                                                                          >  1:4212451/1‑98 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACCGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aa                                                                                       >  1:4198350/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aa                                                                                       >  2:3039795/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aa                                                                                       <  2:2059121/100‑1 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aa                                                                                       >  2:1566301/1‑100 (MQ=255)
        gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aaa                                                                                      >  1:3434000/1‑99 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aaac                                                                                     <  1:1673103/100‑3 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aaacag                                                                                   >  2:923648/1‑96 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a‑aaaa‑aaacag                                                                                   >  2:1392349/1‑96 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑A‑AAAA‑AA‑CA‑Gt                                                                                 >  1:2736245/1‑100 (MQ=255)
               gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATc                                                                              >  2:1714501/1‑100 (MQ=255)
               gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATc                                                                              >  2:4205217/1‑100 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCg                                                                             >  2:3460666/1‑100 (MQ=255)
                 cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGt                                                                            <  2:2816107/100‑1 (MQ=255)
                          ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTa                                                                   >  2:4177039/1‑100 (MQ=255)
                           tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTAc                                                                  <  2:628471/100‑1 (MQ=255)
                            ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACg                                                                 <  2:2029939/100‑1 (MQ=255)
                            ggAACGATTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTTTGAAAAAAAAAA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTAc                                                                  <  1:2197134/99‑1 (MQ=255)
                             gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGTTTTATTGTTTGTTCATCATT‑‑ATAAAA‑AA‑CA‑GTATAGTGTtttttt                                                                    >  1:122545/1‑96 (MQ=255)
                             gAACGCTTAAATGATGCTTCATACAGTCTTCATCAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGc                                                                >  1:2903089/1‑100 (MQ=255)
                                 gCTTAAATGATGCTTCATACAGTCTTTATAGAAGGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAgg                                                            >  1:2898323/1‑100 (MQ=255)
                                 gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAgg                                                            >  2:1568528/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCaa                                                              >  1:322005/1‑97 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAggg                                                           >  1:1255953/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAggg                                                           <  2:3673744/100‑1 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTa                                                         >  2:3394903/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAGCAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAg                                                        >  1:2358778/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTa                                                         >  1:751678/1‑99 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAg                                                        >  2:4316150/1‑100 (MQ=255)
                                      aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGGAAGGGTAGc                                                       >  2:3159157/1‑100 (MQ=255)
                                      aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGc                                                       >  2:225756/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGcc                                                      >  1:1179591/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAA‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGcc                                                      >  2:1934500/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGccc                                                     >  2:731808/1‑100 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGccc                                                     <  2:3434022/100‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGcccc                                                    >  2:1903113/1‑100 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTAGCGTGATATTTTACGCAAGGGTAGcccc                                                    <  2:1916082/100‑1 (MQ=255)
                                          atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGccccc                                                   >  1:4353838/1‑100 (MQ=255)
                                          atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGccccc                                                   >  2:1108965/1‑100 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                  <  2:1845171/100‑1 (MQ=255)
                                             cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                <  1:941622/100‑1 (MQ=255)
                                              ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGt                                               <  2:1531347/100‑1 (MQ=255)
                                                cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                             <  1:1732931/100‑1 (MQ=255)
                                                  tagaGTCTTTATAAAAAGAATGAAACAGGATTTATTGTTTGTTCATTATTA‑A‑AAAA‑AA‑CA‑GTATTGTGATATTTTACGTAAGGGTAGCGCCCATGTTCAt                                           >  1:3670098/4‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                           >  1:2386816/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                           >  2:2924572/1‑100 (MQ=255)
                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                          >  1:960784/1‑100 (MQ=255)
                                                    cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                         <  2:1005684/100‑1 (MQ=255)
                                                    cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                         >  2:3085512/1‑100 (MQ=255)
                                                    cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                         <  1:3178824/100‑1 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                       <  2:1002188/100‑1 (MQ=255)
                                                        ntttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTGTCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                     >  1:2815751/2‑100 (MQ=255)
                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                     <  1:503894/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                  >  1:551682/1‑100 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                  <  2:3947970/100‑1 (MQ=255)
                                                            aaaaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                 <  1:2117240/98‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                <  2:2815768/100‑1 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                 >  2:3981925/1‑99 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                               >  1:666551/1‑100 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                               >  1:2021400/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                              <  1:2010369/100‑1 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                              >  1:2524216/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCACCCATGTTCATTTCCTCTGGTTAg                              >  1:311538/1‑100 (MQ=255)
                                                                aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                             <  1:4044384/100‑1 (MQ=255)
                                                                    aaTGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                         <  2:868452/100‑1 (MQ=255)
                                                                     attaaaCAGGATTTATTGTCTGTGAATCATT‑‑ACAAAA‑AA‑CA‑GTATCGTGATATTTTACGCAATGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                        <  1:1066990/97‑1 (MQ=255)
                                                                     aTGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                        >  1:3766007/1‑100 (MQ=255)
                                                                      tGAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                       <  2:286158/100‑1 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                      <  2:4195969/100‑1 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                      <  2:1062079/100‑1 (MQ=255)
                                                                         aaCAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                    <  1:3039777/100‑1 (MQ=255)
                                                                          aCAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                   <  1:821544/100‑1 (MQ=255)
                                                                           cAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                  <  2:2343854/100‑1 (MQ=255)
                                                                           cAGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                  <  1:3251463/100‑1 (MQ=255)
                                                                            aGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                 <  2:3375252/100‑1 (MQ=255)
                                                                            aGGATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAACCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                 >  1:1109651/1‑100 (MQ=255)
                                                                            aGGATTTATTGTCTGTTCATCATT‑‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                <  2:219828/100‑1 (MQ=255)
                                                                             ggATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGan                <  1:3416895/100‑2 (MQ=255)
                                                                             ggATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                <  2:3529756/100‑1 (MQ=255)
                                                                              gATTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaaa               >  1:3489345/1‑100 (MQ=255)
                                                                              gATTTATTGTCTGTTCATCATT‑‑A‑AAAA‑AA‑AATGTATCGTGGTATTTTACGCATGGGTAGCCACCATGTTCCTTTCTTTTTGTTAGTTGTAGACTATGaga               >  1:1034391/1‑98 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc              <  1:3593398/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc              <  2:3726584/100‑1 (MQ=255)
                                                                                 ttATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTa            <  2:2268474/100‑1 (MQ=255)
                                                                                 ttATTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTa            <  2:3003116/100‑1 (MQ=255)
                                                                                 tgaTTGTTTGTTCATCATT‑‑A‑AAAA‑AAGCA‑GTACCGTGATATTTTACGCAAGGGTTGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc              <  1:1494874/96‑1 (MQ=255)
                                                                                   aTTGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt           >  2:2107861/1‑99 (MQ=255)
                                                                                     tGTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct        <  2:3161517/100‑1 (MQ=255)
                                                                                      gTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc       <  2:3766033/100‑1 (MQ=255)
                                                                                      gTCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc       >  2:4096769/1‑100 (MQ=255)
                                                                                       tCTGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTATATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg      >  1:3898754/1‑100 (MQ=255)
                                                                                         tGTTCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg    >  2:3939239/1‑100 (MQ=255)
                                                                                           ttCATCATTA‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg  >  2:1250557/1‑100 (MQ=255)
                                                                                           ttCATCAAT‑‑A‑AAAAGTA‑CA‑GGATCGTTATAGTTTTCGCAAGGGGAGCCTCCATGCTCGTTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc   <  1:1535139/99‑1 (MQ=255)
                                                                                                    |                                                                                               
TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑A‑AAAA‑AA‑CA‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCG  >  NC_000913/4494777‑4494966

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