Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,998 C→A 16.5% intergenic (+19/+20) 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,9980CA16.5% 407.9 / 68.4 158intergenic (+19/+20)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (73/59);  new base A (11/15);  total (84/74)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

Only 100 of 162 total aligned reads displayed.

TGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTT  >  NC_000913/2364899‑2365097
                                                                                                   |                                                                                                   
tGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTAc                                                                                                     >  1:489023/1‑100 (MQ=255)
tGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTAc                                                                                                     >  2:921796/1‑100 (MQ=255)
   gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                                  <  1:1464195/100‑1 (MQ=255)
    tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                 <  2:914941/100‑1 (MQ=255)
    tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                 <  2:489029/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:1173829/100‑2 (MQ=255)
          cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCat                                                                                           <  2:1057138/100‑2 (MQ=255)
          cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCat                                                                                           <  2:3129948/100‑2 (MQ=255)
          cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCat                                                                                           <  2:3730158/100‑2 (MQ=255)
           tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                          <  2:1776914/100‑3 (MQ=255)
           tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                          <  1:1331527/100‑3 (MQ=255)
            gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                         >  1:3081708/1‑100 (MQ=255)
            gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                         >  1:596141/1‑100 (MQ=255)
            gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                         <  1:951973/100‑4 (MQ=255)
             aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                        <  2:3266486/100‑5 (MQ=39)
               gATCGGGTGCATTATGATTTTAATGTCGCCACGCGAACAACGTTACGTTAGAAAGGTCTTCTGTAACAATGCTGACAGCCAGTAATGGCTGTCatc                                                                                          <  1:1876137/96‑3 (MQ=37)
                  ctagTGCATTATGATTTGAATGTCTCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                   <  1:3367460/97‑1 (MQ=255)
                  cTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                   >  1:1743565/1‑100 (MQ=255)
                  cTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaa                                                                                    <  1:2518637/99‑1 (MQ=255)
                   tGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcattgttac                                                                                  <  2:3288384/100‑11 (MQ=38)
                     gTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAttt                                                                                <  2:1137581/100‑1 (MQ=255)
                       gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                              <  2:3749817/100‑1 (MQ=255)
                       gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                              >  2:2183192/1‑100 (MQ=255)
                       gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                              >  1:944090/1‑100 (MQ=255)
                        cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                             <  2:2533301/100‑1 (MQ=255)
                        cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                             <  2:3262616/100‑1 (MQ=255)
                         aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                            <  1:7054/100‑1 (MQ=255)
                          ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                            >  1:2017039/1‑99 (MQ=255)
                               aTTTGAATGTCGCCATCCGAAAAACGTTAAGTGTTAAAGGGGTTCTGTAACAATGATGACAGGAATTACTGGGTGTCAGATTGAAAATTTAGTGGTTa                                                                        <  1:711253/98‑1 (MQ=255)
                               aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTaa                                                                       >  1:1517022/1‑99 (MQ=255)
                                tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACg                                                                     <  1:1019628/100‑1 (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                                                             >  2:3647034/1‑100 (MQ=255)
                                         cGCCACCCGAAAAACGTTACGGTTGAAAGGTCTTCTGTAACAATGATGACAGCTATGACTGGCTGTAAGATTGAACATTTAGTGGTTATCGAATTaccc                                                             >  1:1315701/1‑99 (MQ=255)
                                          gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCAt                                                           >  2:829522/1‑100 (MQ=255)
                                          gCCACCCGAAAAAAGTTACGTTTGAAAGGTCTTCTGTACCAATGATGACAGCCATTACTGGCTGTCAGATTGCAAATTTAGTGGTTCACGAATTCGCCa                                                            >  1:2359835/1‑99 (MQ=255)
                                              cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                       >  2:1960389/1‑100 (MQ=255)
                                                cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGAt                                                     <  2:2240320/100‑1 (MQ=255)
                                                 gAAAAACGTTACGTTTGAAAGGTCCTCTATAACAAAGATTACAACCATTACTGCCTCTCAAATTGCAAATTTAATGGTTAACGAATTACCCATCGAGATa                                                    >  1:1749518/1‑100 (MQ=255)
                                                 gAAAAAAGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATa                                                    <  2:172833/100‑1 (MQ=255)
                                                    aaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                 >  1:2924574/1‑100 (MQ=255)
                                                      acgttacgttTGAAAGGTCTTCTGTAACAATGATGACAGACATTACTGGCTGTCCGATTGAAAATTTAGTGGTTAACGAATTCACCATCCAGTTACACtt                                               >  1:364436/1‑100 (MQ=255)
                                                        gttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTg                                             >  2:2629988/1‑100 (MQ=255)
                                                          tacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGcc                                           <  2:865189/100‑1 (MQ=255)
                                                            cgttTGAAAGGTCTTCTATAAAAAAGATGACAGCCATTACTGGCTGTCAGATTGAAAACTTAGTGGTTAACGAATTCCCCATCCCGATAAACTTTGCCtt                                         >  1:893052/1‑100 (MQ=255)
                                                                tGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                     >  2:1000167/1‑100 (MQ=255)
                                                                tGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                     >  1:1965090/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                                 >  1:3511603/1‑100 (MQ=255)
                                                                    aGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                                 >  2:14411/1‑100 (MQ=255)
                                                                       tcttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgc                              <  2:605246/100‑1 (MQ=255)
                                                                       tcttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgc                              <  1:985437/100‑1 (MQ=255)
                                                                        cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                             <  2:1760714/100‑1 (MQ=255)
                                                                        cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                             <  2:2072500/100‑1 (MQ=255)
                                                                          tctcTAACACTGCTGACAGCAATTACTGACTGTCAGATTGAAAATTTAGTGGTTAACGCATTCCCCATACAGATAAATTTTGCCTTTTTCCAaatt                               >  1:2227215/1‑95 (MQ=255)
                                                                          tctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                           >  1:3223843/1‑100 (MQ=255)
                                                                            tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                         <  2:1632306/100‑1 (MQ=255)
                                                                            tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                         <  2:1232841/100‑1 (MQ=255)
                                                                             gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                        >  2:1924954/1‑100 (MQ=255)
                                                                             gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                        >  1:89698/1‑100 (MQ=255)
                                                                             gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                        >  1:2223837/1‑100 (MQ=255)
                                                                               tttcaatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaa                      >  1:1931774/9‑100 (MQ=38)
                                                                               aaCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaa                      >  2:2593828/1‑100 (MQ=255)
                                                                                aCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaaa                     <  2:2469035/100‑1 (MQ=255)
                                                                                 tcaatctgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATGCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                    <  2:1572658/94‑1 (MQ=255)
                                                                                 cgatcttgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                    <  2:885037/94‑1 (MQ=255)
                                                                                 cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                    <  2:2315210/100‑1 (MQ=255)
                                                                                 cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                    <  1:14410/100‑1 (MQ=255)
                                                                                  caatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTTGTAAACGAAGTCGCCATCCAGATAAAATTTTCCTTTTACCACTTGCATGGCCaaaa                    >  1:2047308/6‑98 (MQ=16)
                                                                                     gatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCa                  >  1:599209/1‑98 (MQ=255)
                                                                                     gatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAta                >  1:3813007/1‑99 (MQ=255)
                                                                                     gatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATc                >  2:2882869/1‑100 (MQ=255)
                                                                                      atgaCAGCTATTACTGGCTGTGAAATAGAAAATTTAGTGGTTGACGAATTCCCCATCCATATAAACTTTGCCTTTTTCCACGTGCATGACTAAACCATcc               <  1:344250/100‑1 (MQ=255)
                                                                                       tgaCAGCCCTTACTGGCTGTCAGATTGAAAATTTAGTAGTTATCTAATTCCCCATCCAGATCACCTTTGCCTTTTTCCACATGCATGACTAAACCATcc               <  1:1339910/99‑1 (MQ=255)
                                                                                       tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa              >  2:281174/1‑100 (MQ=255)
                                                                                       tgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAcagatcg              >  1:885030/1‑94 (MQ=255)
                                                                                        gaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa              >  1:2256227/1‑99 (MQ=255)
                                                                                        gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAg             >  2:1822583/1‑100 (MQ=255)
                                                                                          cAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGAt           >  2:2985218/1‑100 (MQ=255)
                                                                                            gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa         >  2:943199/1‑100 (MQ=255)
                                                                                            gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATa          >  2:587851/1‑99 (MQ=255)
                                                                                            gCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa         >  2:1329018/1‑100 (MQ=255)
                                                                                               aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCa      >  1:2813041/1‑100 (MQ=255)
                                                                                                  aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt   >  2:3020194/1‑100 (MQ=255)
                                                                                                  aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt   <  1:2947907/100‑1 (MQ=255)
                                                                                                  aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt   >  1:2103342/1‑100 (MQ=255)
                                                                                                  aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt   >  1:837319/1‑100 (MQ=255)
                                                                                                  aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt   <  2:1793125/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTT  >  NC_000913/2364899‑2365097

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