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‑‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  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‑‑aaaaa‑a                                                                                              >  2:2408343/1‑100 (MQ=255)
   aagaAGATCGACATACCTTGCTGGTACGCTTAAATGATACTTCATACAGTCTTTATAAAAAGAATTAAACAGTATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                                 <  1:1549955/96‑1 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑a                                                                                            <  1:2005664/100‑1 (MQ=255)
    tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑a                                                                                            <  2:1924750/100‑1 (MQ=255)
     gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑aa                                                                                           >  2:1208721/1‑99 (MQ=255)
      aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑aac                                                                                          <  2:1922595/100‑3 (MQ=255)
      aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑aa                                                                                           >  1:487643/1‑98 (MQ=255)
       aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑aaca                                                                                         <  2:1157019/100‑4 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑aacag                                                                                        <  1:1858139/100‑5 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑aacag                                                                                        >  2:2197542/1‑96 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa‑a‑aacag                                                                                        <  2:516520/100‑5 (MQ=255)
        gATCGACTTCCCTGGCTGGAACGCTTAAATGATACTTCATACAGTCTTTATAAAAAGAATGAAACAGGATATATTGTCTGTTCATCATT‑‑aaaaa‑a‑aacag                                                                                        >  1:1554745/1‑96 (MQ=255)
               ttCCTTGGCTGGAACTCTTACATGATGCTTCATAAAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAA‑AGACAGTATCGTg                                                                                  <  1:1362193/100‑1 (MQ=255)
               ttCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTg                                                                                  <  2:1788504/100‑1 (MQ=255)
                tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGa                                                                                 >  1:434549/1‑100 (MQ=255)
                 cccTGGCTGGAACGCTTAAATGATGCGTCATACAGTCTTTATAAAAAGAATGAAACATGATTTATTTTCTGTTAATCATTA‑AAAAA‑A‑ACAGTATCGTg                                                                                  <  1:106631/98‑1 (MQ=255)
                  ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAA‑A‑ACAGTATCGTGat                                                                                >  2:1727201/1‑100 (MQ=255)
                   ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGatat                                                                              >  1:2496410/1‑100 (MQ=255)
                    tggctagAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAAAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGa                                                                                 <  1:205413/96‑1 (MQ=255)
                      gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATAtttt                                                                           >  2:321731/1‑100 (MQ=255)
                       ctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTa                                                                          >  1:1187848/1‑100 (MQ=255)
                         ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACg                                                                        >  2:1627096/1‑100 (MQ=255)
                              gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAgg                                                                   <  2:2031870/100‑1 (MQ=255)
                              gCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAgg                                                                   <  1:1441398/100‑1 (MQ=255)
                               cttgactgatgCTTCATACAATCATGATAAAAATAATGAAACAGAATTTATTGTCTTTTTATCATTA‑AAAAA‑A‑ACCGTATCAAGATGTTTTAGCCAAgg                                                                   >  1:288929/7‑99 (MQ=255)
                               cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAggg                                                                  >  2:2308480/1‑100 (MQ=255)
                                ttAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCaa                                                                     <  1:2409809/96‑1 (MQ=255)
                                  aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAg                                                               <  1:2316969/100‑1 (MQ=255)
                                  aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAg                                                               >  1:2549699/1‑100 (MQ=255)
                                     tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            <  1:2490556/100‑1 (MQ=255)
                                     tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGCATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            >  1:1595311/1‑100 (MQ=255)
                                       atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                          >  1:2558161/1‑100 (MQ=255)
                                        tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                         >  2:1348826/1‑100 (MQ=255)
                                         gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATTGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  2:1282521/1‑100 (MQ=255)
                                         gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  1:830349/1‑100 (MQ=255)
                                         gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  2:91109/1‑100 (MQ=255)
                                          cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                       >  2:2167909/1‑100 (MQ=255)
                                           ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGt                                                      >  1:228121/1‑100 (MQ=255)
                                            tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                     <  1:1294863/100‑1 (MQ=255)
                                              agaCAGGCTTTATAAAAAGAAGTAAACAGGATTTATTGTGTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACTCGAGGGTAGCCCCCATGTTCa                                                   <  1:840577/98‑1 (MQ=255)
                                              aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   >  2:1506286/1‑100 (MQ=255)
                                              aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTACAACAAAA‑AAAAA‑A‑ACAGGATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   <  1:1387153/100‑1 (MQ=255)
                                              aTACAGTATTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                     <  2:1086356/99‑1 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTCCGCAAGGGTAGCACCCAt                                                        >  1:528628/1‑94 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                  >  2:1518731/1‑100 (MQ=255)
                                               tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                  >  2:852304/1‑100 (MQ=255)
                                                aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTTTTCATCAATA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                                 <  1:1091786/100‑1 (MQ=255)
                                                  aGTCTTTATAAAAAGAATGAAGCAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                               <  2:683857/100‑1 (MQ=255)
                                                  aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                               <  2:1763608/100‑1 (MQ=255)
                                                   gTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTGTTCAAAAATA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                              <  1:2085997/100‑1 (MQ=255)
                                                   gTCTTTATAAAAAGAAAGAAACAAGAGTCATTGTCTTTTCATCATTA‑AAAAA‑A‑ACAGAATCAATATATTTTACGCAAGGGTAGACCCCagattgtt                                                  >  1:2158585/1‑89 (MQ=255)
                                                    tCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                             <  2:1595320/100‑1 (MQ=255)
                                                      tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                           <  2:562261/100‑1 (MQ=255)
                                                        tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                         >  2:713701/1‑100 (MQ=255)
                                                        tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        <  1:443859/100‑1 (MQ=255)
                                                         ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                        <  2:952981/100‑1 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  2:1130894/100‑1 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  2:1970958/100‑1 (MQ=255)
                                                          taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  1:2308466/100‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTGTGGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  1:1886703/99‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      <  2:685525/100‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      <  2:907884/100‑1 (MQ=255)
                                                           aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCACCATGTTCATTTCCTCTGGTTa                                      >  1:465377/1‑100 (MQ=255)
                                                            aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                     >  2:1839412/1‑100 (MQ=255)
                                                            aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                     >  1:1264796/1‑100 (MQ=255)
                                                             aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                    >  2:918182/1‑100 (MQ=255)
                                                             aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                    <  1:511914/100‑1 (MQ=255)
                                                                gAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                                 >  1:474788/1‑100 (MQ=255)
                                                                gAATGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATgg                                 <  1:1518722/100‑1 (MQ=255)
                                                                gAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                                <  2:1475564/100‑1 (MQ=255)
                                                                  aTGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                               >  2:2075109/1‑100 (MQ=255)
                                                                  aTGAAACAGGATTTATTGTCTGTTAAAAAAAA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                               <  1:642631/100‑1 (MQ=255)
                                                                   tGAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                              >  1:1795172/1‑100 (MQ=255)
                                                                    gAAACAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTTCGCAAGGGTAGCCCCCATGTTCATTTCACCTGGTTAGATg                                  >  1:1152813/1‑95 (MQ=255)
                                                                      aaCAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                            >  1:749361/1‑99 (MQ=255)
                                                                        cAGGATTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         <  1:1506277/100‑1 (MQ=255)
                                                                        cAGGATTTATTGTCTGTTCATCATT‑‑AAAAAGA‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         >  1:2014006/1‑100 (MQ=255)
                                                                         aGGATTTATTGTCTGATCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                        <  2:2558180/100‑1 (MQ=255)
                                                                            aTTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  2:2029002/100‑1 (MQ=255)
                                                                            aTTTATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATGTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  1:1655097/100‑1 (MQ=255)
                                                                               tATTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt                  <  2:431746/100‑1 (MQ=255)
                                                                                aTTGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt                 >  2:308731/1‑100 (MQ=255)
                                                                                  tGTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCTTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct               <  2:2451372/100‑1 (MQ=255)
                                                                                   gTCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc              >  1:2295047/1‑100 (MQ=255)
                                                                                    tCTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcg             >  1:2168488/1‑100 (MQ=255)
                                                                                     cTGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc            >  2:1252790/1‑100 (MQ=255)
                                                                                      tGTTCATCATTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg           >  2:1979966/1‑100 (MQ=255)
                                                                                        ttCATCATTC‑AAAAA‑A‑ACAGTATCGTGATATTTTACCCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgcg         <  2:116563/100‑1 (MQ=255)
                                                                                        ttCATCATT‑‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTTCTCTGGTTAGATGGCGACTATGAACCTATTCTcgcgcg         <  1:1597975/99‑1 (MQ=255)
                                                                                           atcatTA‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt      <  2:2326699/100‑1 (MQ=255)
                                                                                              atT‑‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAg  >  1:1753703/1‑100 (MQ=255)
                                                                                                 |                                                                                                    
CTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAA‑A‑ACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  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: ‑