Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,995 T→G 17.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,9950TG17.0% 230.1 / 40.5 98intergenic (+16/+23)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (45/35);  new base G (7/10);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAG  >  NC_000913/2364897‑2365094
                                                                                                  |                                                                                                   
cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                    >  2:1391790/1‑100 (MQ=255)
cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                    <  2:1943464/100‑1 (MQ=255)
cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTAGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                    >  2:2245763/1‑100 (MQ=255)
cgTGGGTAAAACCTGAAGATCTGGTACATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                    >  2:2113963/1‑100 (MQ=255)
cgTGGGTAAAACCTGAAGATATGGTGAATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAtt                                                                                                    >  2:905374/1‑100 (MQ=255)
 gTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTa                                                                                                   >  1:1760059/1‑100 (MQ=255)
   gggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACt                                                                                                 >  1:2305076/1‑100 (MQ=255)
    ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                >  1:1285497/1‑100 (MQ=255)
    ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                >  2:487149/1‑100 (MQ=255)
     gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                               >  1:164226/1‑100 (MQ=255)
      tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                              <  1:1158603/100‑1 (MQ=255)
      tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                              >  1:1797508/1‑100 (MQ=255)
      tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                              <  2:113260/100‑1 (MQ=255)
      tAAAACCTGAAGATCTCGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                              <  2:1394507/100‑1 (MQ=255)
        aaaCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTg                                                                                            <  2:336942/100‑1 (MQ=255)
         aaCCTGAAGATCTGGTGCATTATGATCTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGt                                                                                           >  1:1447787/1‑100 (MQ=255)
           ccTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCa                                                                                         <  2:1506135/100‑1 (MQ=255)
            cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAg                                                                                        >  2:72688/1‑100 (MQ=255)
             tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                       <  2:1606248/100‑3 (MQ=255)
             tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                       <  2:1993001/100‑3 (MQ=255)
             tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                       <  1:1323814/100‑3 (MQ=255)
              gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                      <  1:320134/100‑4 (MQ=255)
               aaGATCTGGTGTATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                     <  1:328627/100‑5 (MQ=38)
               aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                     <  1:1091745/100‑1 (MQ=255)
               aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                     >  1:315353/1‑100 (MQ=255)
               aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                     <  1:1292495/100‑1 (MQ=255)
               aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                     <  1:1857753/100‑5 (MQ=39)
                aGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcatt                                                                                    <  2:1334718/100‑6 (MQ=39)
                   tCTGGTGCATTATGATTTGAGTGTACCGACCCGAAAAACGTTCCGGTTGAGAGGTCTTGTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaa                                                                                 <  2:834257/100‑1 (MQ=255)
                   tCTGGTGCATTATGATTTGAATCTCGCCACCCGAAAAACGTTACGTTTGAAATGTCTTCTGTAACCATGATGACAGCCAGTAATCGCTGTCatcatt                                                                                    <  2:710437/97‑6 (MQ=255)
                     tGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAt                                                                               >  2:184342/1‑100 (MQ=255)
                       gTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAttt                                                                             <  2:341228/100‑1 (MQ=255)
                        tGCATTATGATTTGAATGTTGCCACCCGAAAAACGTTACGCCTGAAAGGTCTTCTGTATCAATGATGACAGCCATTCCTAGCTGTCAGATTGaaaa                                                                                >  2:711399/1‑96 (MQ=255)
                         gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                           <  1:621539/100‑1 (MQ=255)
                         gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                           <  1:1733893/100‑1 (MQ=255)
                         gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGAAAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                           <  2:1389358/100‑1 (MQ=255)
                          cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                          <  2:619740/100‑1 (MQ=255)
                           aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                         >  1:1335852/1‑100 (MQ=255)
                            ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                          <  1:1559462/98‑1 (MQ=255)
                            ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                        <  1:72689/100‑1 (MQ=255)
                            ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                        >  2:39827/1‑100 (MQ=255)
                              aTGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                      <  1:997687/100‑1 (MQ=255)
                                 aTTTGAATGTCGCCACCCGAAAAACGTTACTTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                   <  2:437648/100‑1 (MQ=255)
                                 aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                   <  1:835520/100‑1 (MQ=255)
                                 aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                   >  1:1774943/1‑100 (MQ=255)
                                 aTTTGAATGTCGCCACACGACAAACGTTACGTTTGAAAGATCTTCTGTACCAATGCTGACAGCCATAACTGGCTGTCAGATTGAAACTTTAATTGTTAAc                                                                   >  2:809169/1‑100 (MQ=255)
                                  tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACg                                                                  <  2:1240848/100‑1 (MQ=255)
                                      aaTGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                              >  1:341219/1‑100 (MQ=255)
                                      aaTGTCGCCACCCGAAAAAAGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                              >  2:1753488/1‑100 (MQ=255)
                                        tGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                            >  2:416090/1‑100 (MQ=255)
                                        tGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                            <  2:1998167/100‑1 (MQ=255)
                                         gTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTccc                                                           >  1:200586/1‑100 (MQ=255)
                                           cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                         >  2:1409156/1‑100 (MQ=255)
                                           cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                         >  1:1601960/1‑100 (MQ=255)
                                            gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCAt                                                        >  2:638422/1‑100 (MQ=255)
                                            gCCACCCGAAAAACGTTACGTTTGAAAGATCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCAt                                                        >  2:2096499/1‑100 (MQ=255)
                                              cACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATcc                                                      >  1:1114900/1‑100 (MQ=255)
                                               aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                     <  1:2113980/100‑1 (MQ=255)
                                                cctttcaaaCGTAACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                    >  1:930889/7‑100 (MQ=255)
                                                cccGAAAAACGTTACGTTTGACAGGTCTTCTGTAACAAGGATGCCATCCATTACTGGCTGTCAGAGTGAAGATTTAGTGGTTAACGAATTCCCCATCCAg                                                    >  2:1463851/1‑100 (MQ=255)
                                                cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                    >  1:1552487/1‑100 (MQ=255)
                                                cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                    <  2:1285490/100‑1 (MQ=255)
                                                   gAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCTCATCCAg                                                    >  2:1355589/1‑97 (MQ=255)
                                                     aaaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATaaa                                               >  1:2314730/1‑100 (MQ=255)
                                                       aaCGTTCCGTTTGAAAGGTCTTCTGTAACAATGATGACCGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCCGATAAACt                                             <  2:824666/100‑1 (MQ=255)
                                                            tacgttTGAAAGGTCTTCTGTACCAATGATGACAGCCATTACTGGCTGTCAGATTGAAATTTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGc                                         >  2:502871/1‑99 (MQ=255)
                                                            tacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGcc                                        <  2:667302/100‑1 (MQ=255)
                                                              cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                      >  2:1647441/1‑100 (MQ=255)
                                                              cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                      >  1:1419370/1‑100 (MQ=255)
                                                                 ttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttttt                                   >  1:2137066/1‑100 (MQ=255)
                                                                   gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttttac                                 >  2:1494904/1‑98 (MQ=255)
                                                                   gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                 <  1:1536561/100‑1 (MQ=255)
                                                                   gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                 <  2:2305055/100‑1 (MQ=255)
                                                                          cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                          >  2:141537/1‑100 (MQ=255)
                                                                            tctGTAACAAGGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                          <  2:612183/98‑1 (MQ=255)
                                                                             ctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGa                       >  1:1671565/1‑100 (MQ=255)
                                                                               gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                     >  1:1889054/1‑100 (MQ=255)
                                                                               gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                     >  2:2381377/1‑100 (MQ=255)
                                                                                  ttcaatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCACCATACATTTAAACTTTGCATTTTTCCACATGCATGAGTaaa                  >  2:453823/8‑100 (MQ=25)
                                                                                   cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                 >  1:1236689/1‑100 (MQ=255)
                                                                                   cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                 >  2:898770/1‑100 (MQ=255)
                                                                                    caatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATTAACTTTGCCTTTTTCCACATTCATGACTaa                   >  2:1411150/6‑97 (MQ=38)
                                                                                     atgatgACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCa               >  2:1410003/1‑100 (MQ=255)
                                                                                     atgatgACAGCCACTACTGTCTGTCAGATTGATAATTTACTTCTTAACGAATTCCCCCCCCAGATACACTTCGCCTTTTTCCACATGCATGACTAAACCa               <  2:1674822/100‑1 (MQ=255)
                                                                                     aatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCa               >  1:2050104/5‑100 (MQ=39)
                                                                                      atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAt              >  1:54288/4‑100 (MQ=255)
                                                                                       tctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCa               >  2:2177458/3‑98 (MQ=255)
                                                                                       gatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATc             >  2:1796210/1‑100 (MQ=255)
                                                                                        ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc            >  1:860342/2‑100 (MQ=255)
                                                                                        ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc            >  1:2038115/2‑100 (MQ=255)
                                                                                        atgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc            <  1:710385/100‑1 (MQ=255)
                                                                                         tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa           <  2:339878/100‑1 (MQ=255)
                                                                                            cAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGAt        <  1:2245783/100‑1 (MQ=255)
                                                                                                cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt     <  2:2200830/99‑1 (MQ=255)
                                                                                                 aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCa   <  1:1033448/100‑1 (MQ=255)
                                                                                                 aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCa   >  1:194773/1‑100 (MQ=255)
                                                                                                 aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCa   <  1:1033489/100‑1 (MQ=255)
                                                                                                  ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg  <  2:2213578/100‑1 (MQ=255)
                                                                                                  ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg  <  1:1117849/100‑1 (MQ=255)
                                                                                                  ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg  >  1:681525/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAG  >  NC_000913/2364897‑2365094

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: