Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,995 T→G 16.0% intergenic (+16/+23) nudI → / ← ais nucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,364,9950TG16.0% 474.2 / 84.8 194intergenic (+16/+23)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (87/75);  new base G (13/18);  total (100/94)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 201 total aligned reads displayed.

GCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAG  >  NC_000913/2364896‑2365094
                                                                                                   |                                                                                                   
gcgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAt                                                                                                     <  1:1007524/100‑1 (MQ=255)
 cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                    >  1:1193057/1‑100 (MQ=255)
  gTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCTACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                    <  1:1072750/99‑1 (MQ=255)
   tGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTAc                                                                                                  <  1:2217269/100‑1 (MQ=255)
     ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                >  2:812858/1‑100 (MQ=255)
     ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                >  1:3745283/1‑100 (MQ=255)
     ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                >  1:2818841/1‑100 (MQ=255)
     ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAACCATTACTg                                                                                                >  1:3461063/1‑100 (MQ=255)
      gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                               <  2:4339587/100‑1 (MQ=255)
       tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                              >  2:898967/1‑100 (MQ=255)
       tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                              <  2:1842454/100‑1 (MQ=255)
        aaaaCCTGAAGAGCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                <  1:1264907/97‑1 (MQ=255)
         aaaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAATGTCTTCTGTAACAATGATGACAGCCATTACTGGCTg                                                                                            <  2:3104772/100‑1 (MQ=255)
         aaaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTg                                                                                            <  2:3706734/100‑1 (MQ=255)
           aCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATTATGACAGCCAGTAATGGCTGTc                                                                                          <  1:2548848/100‑1 (MQ=255)
            ccTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCa                                                                                         <  2:2024285/100‑1 (MQ=255)
             cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAg                                                                                        >  1:248127/1‑100 (MQ=255)
             cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAg                                                                                        >  1:4026655/1‑100 (MQ=255)
              tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                       <  2:2269183/100‑3 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                      >  1:1479280/1‑100 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                      <  2:750574/100‑4 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                      <  1:2008304/100‑4 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                      <  2:2381366/100‑4 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAATACTGGCTGTCAGAt                                                                                      <  1:1379118/100‑1 (MQ=255)
                aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                     <  2:1468080/100‑5 (MQ=39)
                aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                     <  1:3677846/100‑5 (MQ=39)
                aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                     <  2:191311/100‑5 (MQ=39)
                 aGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTg                                                                                    >  2:3855137/1‑100 (MQ=255)
                    tCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcattgtt                                                                                 <  2:779124/100‑9 (MQ=38)
                     cTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGTCATTACTGGCTGTCAGATTGaaaa                                                                                >  2:3669973/1‑100 (MQ=255)
                      tGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAt                                                                               <  2:3194355/100‑1 (MQ=255)
                       ggTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAtt                                                                              <  2:671939/100‑1 (MQ=255)
                       ggTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAtt                                                                              <  2:1633190/100‑1 (MQ=255)
                        gTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAttt                                                                             <  1:2270775/100‑1 (MQ=255)
                         tGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTa                                                                            >  1:754119/1‑100 (MQ=255)
                          gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                           <  1:4118740/100‑1 (MQ=255)
                          gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                           >  2:1298780/1‑100 (MQ=255)
                           cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                          <  2:1747431/100‑1 (MQ=255)
                             ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                        <  1:2657306/100‑1 (MQ=255)
                                   tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACg                                                                  <  2:4089412/100‑1 (MQ=255)
                                   tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACg                                                                  <  1:2185400/100‑1 (MQ=255)
                                   tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCAATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACg                                                                  >  1:70112/1‑100 (MQ=255)
                                      gAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAt                                                               >  1:2412304/1‑100 (MQ=255)
                                       aaTGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                              >  2:3055251/1‑100 (MQ=255)
                                            cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                         >  1:3375230/1‑100 (MQ=255)
                                            cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                         >  1:2902693/1‑100 (MQ=255)
                                                aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                     >  1:330575/1‑100 (MQ=255)
                                                aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                     <  1:3022890/100‑1 (MQ=255)
                                                aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                     >  1:2964353/1‑100 (MQ=255)
                                                aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                     >  1:2521457/1‑100 (MQ=255)
                                                aCCCGAAAAACGTTACGGTTGAAAGGTCTTCTGTAACAAAGATGACAGCCATTACTGGCTGTCTGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                     >  1:503913/1‑100 (MQ=255)
                                                 cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                    <  2:2479081/100‑1 (MQ=255)
                                                 cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTAGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                     <  1:1905617/99‑1 (MQ=255)
                                                 cacGAAAAACGTTACGTTTGAAAGGTTTTCTGTAACACTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                    <  2:3760869/98‑1 (MQ=255)
                                                   cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGAt                                                  >  1:679516/1‑100 (MQ=255)
                                                   cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGAt                                                  <  1:471621/100‑1 (MQ=255)
                                                   cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGAt                                                  <  2:2123675/100‑1 (MQ=255)
                                                      aaaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATaaa                                               <  2:267382/100‑1 (MQ=255)
                                                          cgttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAACTTTAGTGGTTAACGAATTCCCCATCCAGATAAACtt                                            >  1:4295393/1‑99 (MQ=255)
                                                           gttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTg                                          >  2:3914431/1‑100 (MQ=255)
                                                            ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGc                                         <  2:4183614/100‑1 (MQ=255)
                                                               cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                      >  2:1422103/1‑100 (MQ=255)
                                                                gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttt                                     >  2:559153/1‑100 (MQ=255)
                                                                 tttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttt                                    >  1:614473/1‑100 (MQ=255)
                                                                    gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                 <  2:1128401/100‑1 (MQ=255)
                                                                    gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                 >  2:479657/1‑100 (MQ=255)
                                                                    gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                 <  1:3623665/100‑1 (MQ=255)
                                                                     aaaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCca                                <  1:1158227/100‑1 (MQ=255)
                                                                     aaaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCca                                >  1:3999452/1‑100 (MQ=255)
                                                                      aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcac                               >  2:2846602/1‑100 (MQ=255)
                                                                       aGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                              <  2:3127702/100‑1 (MQ=255)
                                                                         gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatg                            <  2:4056005/100‑1 (MQ=255)
                                                                           cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                          >  1:1873888/1‑100 (MQ=255)
                                                                             tctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                        >  2:3354021/1‑100 (MQ=255)
                                                                               tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                      <  1:4249008/100‑1 (MQ=255)
                                                                                 tAACAATGATGACAGCGATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTa                    <  1:2578289/100‑1 (MQ=255)
                                                                                 tAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTa                    <  1:1301284/100‑1 (MQ=255)
                                                                                    cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTTCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                 >  1:3295227/1‑100 (MQ=255)
                                                                                    cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                 <  2:3586051/100‑1 (MQ=255)
                                                                                     aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAcc                >  1:3202331/1‑100 (MQ=255)
                                                                                      aatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCGCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCa               >  1:4284358/5‑100 (MQ=38)
                                                                                       tgatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAt              >  1:4241691/1‑100 (MQ=255)
                                                                                       tgatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAt              <  2:3396119/100‑1 (MQ=255)
                                                                                       atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAt              >  2:1343166/4‑100 (MQ=255)
                                                                                         ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc            >  2:1484613/2‑100 (MQ=255)
                                                                                         atgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc            >  1:544504/1‑100 (MQ=255)
                                                                                         atgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc            <  1:2318859/100‑1 (MQ=255)
                                                                                          tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc            >  1:394442/1‑99 (MQ=255)
                                                                                          tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa           >  2:2289084/1‑100 (MQ=255)
                                                                                           gaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTGTTGCTTTTGCATGCATAAACCa               >  1:4179904/1‑95 (MQ=255)
                                                                                           gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAcagatcgg          >  2:3406519/1‑93 (MQ=255)
                                                                                           gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAg          >  2:1867716/1‑100 (MQ=255)
                                                                                            aCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGa         >  1:121121/1‑100 (MQ=255)
                                                                                            aCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGa         >  2:1615878/1‑100 (MQ=255)
                                                                                             cAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGAt        >  1:3340131/1‑100 (MQ=255)
                                                                                               gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa      >  2:3314967/1‑100 (MQ=255)
                                                                                                 cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATc    <  1:1020646/100‑1 (MQ=255)
                                                                                                  aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAATGTATTCTTCATTCAGTTATACTTTGTCTTTTTGCACATGCATGACTAAACCATTCTGCTAAt     >  1:1115012/1‑98 (MQ=255)
                                                                                                   ttACTGGCTGTTAGATTGAAAATTTAGTGGTTAAGGAATTCCTCATCCAGATAAACTTGGCCTTTTTCCACATGCATGACTAAAACATCCAGATAATCAg  <  1:2394874/100‑1 (MQ=255)
                                                                                                   ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg  <  2:3525794/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAG  >  NC_000913/2364896‑2365094

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: