Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,995 T→G 17.5% 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,9950TG17.5% 352.8 / 70.7 154intergenic (+16/+23)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (71/54);  new base G (11/16);  total (82/72)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 156 total aligned reads displayed.

GCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCA  >  NC_000913/2364896‑2365093
                                                                                                   |                                                                                                  
gcgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAt                                                                                                    >  2:1690082/1‑100 (MQ=255)
 cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                   >  1:1767280/1‑100 (MQ=255)
  gTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTa                                                                                                  <  2:266672/100‑1 (MQ=255)
   tGGGTAAACCCTGAAGATCTTGTGCATTATTATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGCTCTTCTGTAACAATGATGACAGCCATTAc                                                                                                 >  1:136539/1‑100 (MQ=255)
   tGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTAc                                                                                                 >  1:489023/1‑100 (MQ=255)
   tGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTAc                                                                                                 >  2:921796/1‑100 (MQ=255)
     ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                               >  1:964385/1‑100 (MQ=255)
      gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                              <  1:1464195/100‑1 (MQ=255)
       tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTGTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                             >  1:285790/1‑100 (MQ=255)
       tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                             <  2:489029/100‑1 (MQ=255)
       tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                             <  2:914941/100‑1 (MQ=255)
       tAAAACCTGAAGATCCAGTGCTTTATTATTTGAATGAAAAAAAAAAAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                             <  1:3294026/100‑1 (MQ=255)
         aaaCGTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTg                                                                                           <  1:2699948/100‑1 (MQ=255)
         aaaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTg                                                                                           >  2:2643441/1‑100 (MQ=255)
          aaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGt                                                                                          >  1:3045223/1‑100 (MQ=255)
          aaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGt                                                                                          >  1:3266468/1‑100 (MQ=255)
           aCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTc                                                                                         >  2:3330825/1‑100 (MQ=255)
             cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCat                                                                                       <  2:3730158/100‑2 (MQ=255)
             cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCat                                                                                       <  2:874618/100‑2 (MQ=255)
              tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                      <  2:1776914/100‑3 (MQ=255)
              tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                      <  1:1331527/100‑3 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAACAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                     <  1:1885059/100‑4 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                     >  1:3081708/1‑100 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                     >  1:596141/1‑100 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                     <  1:302950/100‑1 (MQ=255)
               gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                     <  1:951973/100‑4 (MQ=255)
                aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                    <  2:3266486/100‑5 (MQ=39)
                  gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcattg                                                                                  <  2:3270565/100‑7 (MQ=38)
                  gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcattg                                                                                  <  2:3547162/100‑7 (MQ=38)
                  gATCGGGTGCATTATGATTTTAATGTCGCCACGCGAACAACGTTACGTTAGAAAGGTCTTCTGTAACAATGCTGACAGCCAGTAATGGCTGTCatc                                                                                      <  1:1876137/96‑3 (MQ=37)
                     ctagTGCATTATGATTTGAATGTCTCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                               <  1:3367460/97‑1 (MQ=255)
                        gTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAttt                                                                            <  2:1137581/100‑1 (MQ=255)
                          gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                          >  2:2183192/1‑100 (MQ=255)
                           cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                         <  2:2533301/100‑1 (MQ=255)
                            aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                        <  1:7054/100‑1 (MQ=255)
                             ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                       <  2:1821773/100‑1 (MQ=255)
                             ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                       <  1:437577/100‑1 (MQ=255)
                                   tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACg                                                                 <  1:1019628/100‑1 (MQ=255)
                                   tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACg                                                                 >  1:3543503/1‑100 (MQ=255)
                                   tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTAGTTAACg                                                                 >  1:380839/1‑100 (MQ=255)
                                       aaTGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                             >  1:1598143/1‑100 (MQ=255)
                                        aTGTCGCCACCCGAAAAACGTTACGTTTGAACGGTCTTCTGTAACAATGATGACAGCCCTTACTTGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTc                                                            >  1:2819535/1‑100 (MQ=255)
                                         tGTCGCCACCCGAATAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCACTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                           <  1:1690073/100‑1 (MQ=255)
                                          gTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTccc                                                          >  2:1779359/1‑100 (MQ=255)
                                           tCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcccc                                                         <  1:1986332/100‑1 (MQ=255)
                                            cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                        >  2:1897584/1‑100 (MQ=255)
                                            cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                        >  1:837462/1‑100 (MQ=255)
                                            cGCCACCCGAAAAACGTTACGGTTGAAAGGTCTTCTGTAACAATGATGACAGCTATGACTGGCTGTAAGATTGAACATTTAGTGGTTATCGAATTaccc                                                         >  1:1315701/1‑99 (MQ=255)
                                             gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCAt                                                       >  2:829522/1‑100 (MQ=255)
                                                 cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                   >  2:1960389/1‑100 (MQ=255)
                                                  ccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGa                                                  >  1:911946/1‑100 (MQ=255)
                                                  ccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGa                                                  >  1:2443921/1‑100 (MQ=255)
                                                  ccGAAAAACGTTACGTTTGAAAGGTCTGCTGTAACAATGATGACAGCCATTACTGGCTGCCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                   >  1:174197/1‑99 (MQ=255)
                                                   cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGAt                                                 <  2:2240320/100‑1 (MQ=255)
                                                    gAAAAACGTTACGTTTGAAAGGTCCTCTATAACAAAGATTACAACCATTACTGCCTCTCAAATTGCAAATTTAATGGTTAACGAATTACCCATCGAGATa                                                >  1:1749518/1‑100 (MQ=255)
                                                    gAAAAAAGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATa                                                <  2:172833/100‑1 (MQ=255)
                                                      aaaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATaaa                                              <  1:654212/100‑1 (MQ=255)
                                                           gttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTg                                         >  2:2629988/1‑100 (MQ=255)
                                                               cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                     >  2:2827801/1‑100 (MQ=255)
                                                               cgttTGAAAGGTCTTCTATAAAAAAGATGACAGCCATTACTGGCTGTCAGATTGAAAACTTAGTGGTTAACGAATTCCCCATCCCGATAAACTTTGCCtt                                     >  1:893052/1‑100 (MQ=255)
                                                                   tGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                 >  1:1965090/1‑100 (MQ=255)
                                                                   tGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                 >  2:1000167/1‑100 (MQ=255)
                                                                    gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCTAGATAAACTTTGCCTTTTTcc                                >  1:2390202/1‑100 (MQ=255)
                                                                    gAAAGGTCTGCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                >  1:2663941/1‑100 (MQ=255)
                                                                      aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcac                              <  1:1169259/100‑1 (MQ=255)
                                                                       aGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                             >  2:14411/1‑100 (MQ=255)
                                                                       aGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                             >  1:3511603/1‑100 (MQ=255)
                                                                           cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                         <  2:2072500/100‑1 (MQ=255)
                                                                           cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                         >  2:1133086/1‑100 (MQ=255)
                                                                           cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                         <  2:1760714/100‑1 (MQ=255)
                                                                             tctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                       >  2:1193232/1‑100 (MQ=255)
                                                                             tctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                       >  2:3696150/1‑100 (MQ=255)
                                                                             tctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                       >  1:3223843/1‑100 (MQ=255)
                                                                             tctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                       <  2:3223861/100‑1 (MQ=255)
                                                                               tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                     <  2:1632306/100‑1 (MQ=255)
                                                                               tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                     <  2:1232841/100‑1 (MQ=255)
                                                                                gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                    >  1:2223837/1‑100 (MQ=255)
                                                                                 tAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTa                   <  2:1005971/100‑1 (MQ=255)
                                                                                  aaCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaa                  >  2:2593828/1‑100 (MQ=255)
                                                                                   aCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTTCACATGCATGACTaaa                 >  1:786462/1‑100 (MQ=255)
                                                                                   aCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                 <  2:2469035/100‑1 (MQ=255)
                                                                                   aCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                 >  1:2803523/1‑100 (MQ=255)
                                                                                    tcaatctgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATGCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                <  2:1572658/94‑1 (MQ=255)
                                                                                    cgatcttgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                <  2:885037/94‑1 (MQ=255)
                                                                                    cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                <  2:2315210/100‑1 (MQ=255)
                                                                                     caatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTTGTAAACGAAGTCGCCATCCAGATAAAATTTTCCTTTTACCACTTGCATGGCCaaaa                >  1:2047308/6‑98 (MQ=16)
                                                                                      atgatgACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCa              <  1:151103/100‑1 (MQ=255)
                                                                                      aatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAcaa              >  1:894570/5‑98 (MQ=39)
                                                                                      aatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAGTCCCCATACAGATAAACTTTGCCTTTTTCCACATGCATCACTAAACCa              >  1:2176154/5‑100 (MQ=38)
                                                                                       tgatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAt             <  2:3813032/100‑1 (MQ=255)
                                                                                       atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAt             >  2:31101/4‑100 (MQ=255)
                                                                                        gatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATc            >  2:2882869/1‑100 (MQ=255)
                                                                                          tgaCAGCCCTTACTGGCTGTCAGATTGAAAATTTAGTAGTTATCTAATTCCCCATCCAGATCACCTTTGCCTTTTTCCACATGCATGACTAAACCATcc           <  1:1339910/99‑1 (MQ=255)
                                                                                           gaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa          >  1:2256227/1‑99 (MQ=255)
                                                                                           gaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAg         >  1:2310832/1‑100 (MQ=255)
                                                                                            aCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGa        <  1:3696127/100‑1 (MQ=255)
                                                                                             cAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGAt       >  2:2985218/1‑100 (MQ=255)
                                                                                               gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa     >  2:6362/1‑100 (MQ=255)
                                                                                               gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa     >  2:943199/1‑100 (MQ=255)
                                                                                               gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATa      >  2:587851/1‑99 (MQ=255)
                                                                                                   |                                                                                                  
GCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCA  >  NC_000913/2364896‑2365093

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: