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.A92.4% 249.2 / 13.3 79coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base A (37/35);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  NC_000913/4494780‑4494973
                                                                                                 |                                                                                                    
ctaatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                                  <  2:1771789/97‑1 (MQ=255)
  gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                                <  1:1960797/100‑1 (MQ=255)
  gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                                >  1:2389757/1‑100 (MQ=255)
   atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCCTCATT‑‑aaaaaa                                                                                              >  2:2408343/1‑100 (MQ=255)
   aagaAGATCGACATACCTTGCTGGTACGCTTAAATGATACTTCATACAGTCTTTATAAAAAGAATTAAACAGTATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                                 <  1:1549955/96‑1 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                            <  1:2005664/100‑1 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                            <  2:1924750/100‑1 (MQ=255)
     gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaa                                                                                           >  2:1208721/1‑99 (MQ=255)
      aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                          <  2:1922595/100‑3 (MQ=255)
      aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaa                                                                                           >  1:487643/1‑98 (MQ=255)
       aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaaca                                                                                         <  2:1157019/100‑4 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                        <  1:1858139/100‑5 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                        >  2:2197542/1‑96 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                        <  2:516520/100‑5 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATACTTCATACAGTCTTTATAAAAAGAATGAAACAGGATATATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                        >  1:1554745/1‑96 (MQ=255)
               ttCCTTGGCTGGAACTCTTACATGATGCTTCATAAAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAGACAGTATCGTg                                                                                  <  1:1362193/100‑1 (MQ=255)
               ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTg                                                                                  <  2:1788504/100‑1 (MQ=255)
                tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                                 >  1:434549/1‑100 (MQ=255)
                 cccTGGCTGGAACGCTTAAATGATGCGTCATACAGTCTTTATAAAAAGAATGAAACATGATTTATTTTCTGTTAATCATTAAAAAAAACAGTATCGTg                                                                                  <  1:106631/98‑1 (MQ=255)
                  ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGat                                                                                >  2:1727201/1‑100 (MQ=255)
                   ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGatat                                                                              >  1:2496410/1‑100 (MQ=255)
                    tggctagAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAAAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                                 <  1:205413/96‑1 (MQ=255)
                      gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtttt                                                                           >  2:321731/1‑100 (MQ=255)
                       ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTa                                                                          >  1:1187848/1‑100 (MQ=255)
                         ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                        >  2:1627096/1‑100 (MQ=255)
                              gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAgg                                                                   <  2:2031870/100‑1 (MQ=255)
                              gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAgg                                                                   <  1:1441398/100‑1 (MQ=255)
                               cttgactgatgCTTCATACAATCATGATAAAAATAATGAAACAGAATTTATTGTCTTTTTATCATTAAAAAAAACCGTATCAAGATGTTTTAGCCAAgg                                                                   >  1:288929/7‑99 (MQ=255)
                               cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                                  >  2:2308480/1‑100 (MQ=255)
                                ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCaa                                                                     <  1:2409809/96‑1 (MQ=255)
                                  aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                               <  1:2316969/100‑1 (MQ=255)
                                  aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                               >  1:2549699/1‑100 (MQ=255)
                                     tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            <  1:2490556/100‑1 (MQ=255)
                                     tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGCATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            >  1:1595311/1‑100 (MQ=255)
                                       atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                          >  1:2558161/1‑100 (MQ=255)
                                        tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                         >  2:1348826/1‑100 (MQ=255)
                                         gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATTGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  2:1282521/1‑100 (MQ=255)
                                         gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  1:830349/1‑100 (MQ=255)
                                         gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  2:91109/1‑100 (MQ=255)
                                          cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                       >  2:2167909/1‑100 (MQ=255)
                                           ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGt                                                      >  1:228121/1‑100 (MQ=255)
                                            tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                     <  1:1294863/100‑1 (MQ=255)
                                              agaCAGGCTTTATAAAAAGAAGTAAACAGGATTTATTGTGTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACTCGAGGGTAGCCCCCATGTTCa                                                   <  1:840577/98‑1 (MQ=255)
                                              aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   >  2:1506286/1‑100 (MQ=255)
                                              aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTACAACAAAAAAAAAAACAGGATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   <  1:1387153/100‑1 (MQ=255)
                                              aTACAGTATTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                     <  2:1086356/99‑1 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTCCGCAAGGGTAGCACCCAt                                                        >  1:528628/1‑94 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                  >  2:1518731/1‑100 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                  >  2:852304/1‑100 (MQ=255)
                                                aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTTTTCATCAATAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                                 <  1:1091786/100‑1 (MQ=255)
                                                  aGTCTTTATAAAAAGAATGAAGCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                               <  2:683857/100‑1 (MQ=255)
                                                  aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                               <  2:1763608/100‑1 (MQ=255)
                                                   gTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTGTTCAAAAATAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                              <  1:2085997/100‑1 (MQ=255)
                                                   gTCTTTATAAAAAGAAAGAAACAAGAGTCATTGTCTTTTCATCATTAAAAAAAACAGAATCAATATATTTTACGCAAGGGTAGACCCCagattgtt                                                  >  1:2158585/1‑89 (MQ=255)
                                                    tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                             <  2:1595320/100‑1 (MQ=255)
                                                      tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                           <  2:562261/100‑1 (MQ=255)
                                                        tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                         >  2:713701/1‑100 (MQ=255)
                                                        tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        <  1:443859/100‑1 (MQ=255)
                                                         ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        <  2:952981/100‑1 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  2:1130894/100‑1 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  2:1970958/100‑1 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  1:2308466/100‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTGTGGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  1:1886703/99‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      <  2:685525/100‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      <  2:907884/100‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTa                                      >  1:465377/1‑100 (MQ=255)
                                                            aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                     >  2:1839412/1‑100 (MQ=255)
                                                            aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                     >  1:1264796/1‑100 (MQ=255)
                                                             aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                    >  2:918182/1‑100 (MQ=255)
                                                             aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                    <  1:511914/100‑1 (MQ=255)
                                                                gAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                                 >  1:474788/1‑100 (MQ=255)
                                                                gAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                                 <  1:1518722/100‑1 (MQ=255)
                                                                gAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                                <  2:1475564/100‑1 (MQ=255)
                                                                  aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                               >  2:2075109/1‑100 (MQ=255)
                                                                  aTGAAACAGGATTTATTGTCTGTTAAAAAAAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                               <  1:642631/100‑1 (MQ=255)
                                                                   tGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                              >  1:1795172/1‑100 (MQ=255)
                                                                    gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTTCGCAAGGGTAGCCCCCATGTTCATTTCACCTGGTTAGATg                                  >  1:1152813/1‑95 (MQ=255)
                                                                      aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                            >  1:749361/1‑99 (MQ=255)
                                                                        cAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         <  1:1506277/100‑1 (MQ=255)
                                                                        cAGGATTTATTGTCTGTTCATCATT‑‑AAAAAGAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         >  1:2014006/1‑100 (MQ=255)
                                                                         aGGATTTATTGTCTGATCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                        <  2:2558180/100‑1 (MQ=255)
                                                                            aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  2:2029002/100‑1 (MQ=255)
                                                                            aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATGTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  1:1655097/100‑1 (MQ=255)
                                                                               tATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt                  <  2:431746/100‑1 (MQ=255)
                                                                                aTTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt                 >  2:308731/1‑100 (MQ=255)
                                                                                  tGTCTGTTCATCATTAAAAAAAACAGTATCTTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct               <  2:2451372/100‑1 (MQ=255)
                                                                                   gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc              >  1:2295047/1‑100 (MQ=255)
                                                                                    tCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg             >  1:2168488/1‑100 (MQ=255)
                                                                                     cTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc            >  2:1252790/1‑100 (MQ=255)
                                                                                      tGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg           >  2:1979966/1‑100 (MQ=255)
                                                                                        ttCATCATTCAAAAAAACAGTATCGTGATATTTTACCCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg         <  2:116563/100‑1 (MQ=255)
                                                                                        ttCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTTCTCTGGTTAGATGGCGACTATGAACCTATTCTcgcgcg         <  1:1597975/99‑1 (MQ=255)
                                                                                           atcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt      <  2:2326699/100‑1 (MQ=255)
                                                                                              atT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAg  >  1:1753703/1‑100 (MQ=255)
                                                                                                 |                                                                                                    
CTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  NC_000913/4494780‑4494973

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: