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.A96.9% 463.6 / 6.0 131coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base A (71/56);  total (72/59)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494777‑4494970
                                                                                                     |                                                                                                
tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCAtt                                                                                                   >  2:2923881/1‑100 (MQ=255)
  tCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aa                                                                                               <  1:1707188/100‑1 (MQ=255)
   cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                              >  2:2579055/1‑100 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                             <  1:2366798/100‑1 (MQ=255)
    tgatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaa                                                                                             <  1:519521/100‑1 (MQ=255)
     gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                            <  1:906932/100‑1 (MQ=255)
     gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                            >  2:1752776/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                         >  2:1227921/1‑100 (MQ=255)
       tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                         >  2:1630588/1‑100 (MQ=255)
        gaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaa                                                                                        >  1:315406/1‑99 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTTTAAAAAGAATGAAACAGGTTTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                         >  2:402625/1‑98 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  1:2815483/100‑3 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  1:215171/100‑3 (MQ=255)
         aaGATCGACTTCCCTGGCTGGAACGCTGAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                       <  2:618823/100‑3 (MQ=255)
          aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaaca                                                                                      >  2:507704/1‑97 (MQ=255)
           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                     >  1:3028781/1‑96 (MQ=255)
            aTCGACTTCCCTGGCTGGAACGCTTAAATGATTCTTCATACAGTCTTCATAAAAAGAATGAAACAGGATTTATTGTCTGTTAATCATT‑‑aaaaaaaacag                                                                                     >  2:1259028/1‑95 (MQ=255)
             tCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTa                                                                                    >  1:1582210/1‑100 (MQ=255)
                aCTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCg                                                                                 <  2:158825/100‑1 (MQ=255)
                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                              <  1:370196/100‑1 (MQ=255)
                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                              >  2:1456334/1‑100 (MQ=255)
                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                              <  1:2579067/100‑1 (MQ=255)
                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATATATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGa                                                                              >  2:2884306/1‑100 (MQ=255)
                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACACTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTACATCATTAAAAAAAACAGTATCGt                                                                                >  2:934906/1‑98 (MQ=255)
                    cccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGat                                                                             >  2:2937892/1‑100 (MQ=255)
                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGatat                                                                           >  2:1862196/1‑100 (MQ=255)
                       tggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtt                                                                          <  1:2936099/100‑1 (MQ=255)
                        ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAttt                                                                         >  2:2892468/1‑100 (MQ=255)
                         gctggAACGCTTAAATGATGCTTCATACAGTCTTTATATAAAGAATGAAACAGGATTTGTTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtttt                                                                        <  2:1783143/100‑1 (MQ=255)
                         gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATAtttt                                                                        <  1:475924/100‑1 (MQ=255)
                            ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAc                                                                      >  2:2251052/1‑99 (MQ=255)
                            ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     <  1:1827166/100‑1 (MQ=255)
                            ggAACGCTTAAATGATGCTTCAGACAGTCTTTATAAAAAGAATGAAACAGGATGTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACg                                                                     >  2:1905630/1‑100 (MQ=255)
                              aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCa                                                                   >  2:826845/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  2:1386751/1‑100 (MQ=255)
                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAggg                                                               >  2:645635/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                             >  2:3000264/1‑100 (MQ=255)
                                    tAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTa                                                             >  2:3177250/1‑100 (MQ=255)
                                     aaaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAg                                                            >  1:1275208/1‑100 (MQ=255)
                                      aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGc                                                           >  1:1144550/1‑100 (MQ=255)
                                       atgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                          <  2:1124698/100‑1 (MQ=255)
                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                         <  1:2174542/100‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        <  2:2188576/100‑1 (MQ=255)
                                         gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGcccc                                                        <  1:2367557/100‑1 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAATAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:816314/100‑1 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  1:2165353/100‑1 (MQ=255)
                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:2972536/100‑1 (MQ=255)
                                            gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  1:577781/1‑100 (MQ=255)
                                            gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     <  1:2934593/100‑1 (MQ=255)
                                            gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  1:2584138/1‑100 (MQ=255)
                                              ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGt                                                   >  2:407379/1‑100 (MQ=255)
                                               tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                  >  2:903737/1‑100 (MQ=255)
                                               tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                  <  1:903739/100‑1 (MQ=255)
                                                cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  1:2910076/97‑1 (MQ=255)
                                                 aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                >  2:943108/1‑100 (MQ=255)
                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATAAAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                <  2:2869787/100‑1 (MQ=255)
                                                  tACAGTATTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATG‑‑AAAAAGAACAGTATGGTGATATTTTACTCAAGGGTAGCCCCCATGTGCAt                                               >  2:269146/1‑100 (MQ=255)
                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                              <  2:646857/100‑1 (MQ=255)
                                                    cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                             >  1:2874990/1‑100 (MQ=255)
                                                    cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                             <  1:2653747/100‑1 (MQ=255)
                                                     aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                            <  2:496584/100‑1 (MQ=255)
                                                     aGTCTTTATAAAAAGAATGAAACAGGATTAATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGAGATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                            <  2:1974581/100‑1 (MQ=255)
                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                           >  1:1111727/1‑100 (MQ=255)
                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                         >  2:367876/1‑100 (MQ=255)
                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                        <  1:1230254/100‑1 (MQ=255)
                                                          ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                       <  1:993463/100‑1 (MQ=255)
                                                          ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGATAGCCCCCATGTTCATTTCCTCTg                                       >  2:1393734/1‑100 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAGGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  1:1819507/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  2:315403/100‑1 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTAAGCAAGGTTAGCCCCCATGTTCATTTCCTCTgg                                      >  2:2192037/1‑100 (MQ=255)
                                                           tataAAAAGAATGAAACAGGATTTATTGCCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      >  1:1852285/1‑100 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                     >  2:700126/1‑100 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                     <  2:2415891/100‑1 (MQ=255)
                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCACACATGTTCATTTCATCTGGt                                     >  2:1941226/1‑100 (MQ=255)
                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    <  2:2584126/100‑1 (MQ=255)
                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                   >  2:2694827/1‑100 (MQ=255)
                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACTCAAGGGTAACCCCCATGTTCATTTCCTCTGGTTAg                                  >  2:924337/1‑100 (MQ=255)
                                                                  aGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATg                               <  1:3144513/100‑1 (MQ=255)
                                                                    aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             >  2:456300/1‑100 (MQ=255)
                                                                     aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                            >  1:698942/1‑100 (MQ=255)
                                                                      tGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                           >  2:9617/1‑100 (MQ=255)
                                                                       gAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                          <  1:2002292/100‑1 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGTCGACt                         >  2:3064131/1‑100 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  2:412012/1‑100 (MQ=255)
                                                                        aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  2:1351964/1‑100 (MQ=255)
                                                                         aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                        >  1:1163947/1‑100 (MQ=255)
                                                                             ggATTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     >  1:2016062/1‑99 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  1:2923897/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  <  1:1534041/100‑1 (MQ=255)
                                                                               aTTTATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                  >  2:1166535/1‑100 (MQ=255)
                                                                                  tATTGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt               <  2:1582206/100‑1 (MQ=255)
                                                                                    ttGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtc             >  1:741247/1‑100 (MQ=255)
                                                                                     tGTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct            >  1:71241/1‑100 (MQ=255)
                                                                                      gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           >  1:1920662/1‑100 (MQ=255)
                                                                                      gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           >  1:964173/1‑100 (MQ=255)
                                                                                      gTCTGTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAACCACCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           >  2:584404/1‑100 (MQ=255)
                                                                                          gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  2:1614807/1‑100 (MQ=255)
                                                                                          gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  1:118791/1‑100 (MQ=255)
                                                                                          gTTCATCATTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       >  2:2133078/1‑100 (MQ=255)
                                                                                               tcatTAAAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  2:2329865/100‑1 (MQ=255)
                                                                                                     |                                                                                                
TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494777‑4494970

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