Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,998 C→A 15.7% 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,9980CA15.7% 270.4 / 36.9 102intergenic (+19/+20)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (48/38);  new base A (7/9);  total (55/47)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

Only 100 of 103 total aligned reads displayed.

GGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGT  >  NC_000913/2364900‑2365096
                                                                                                  |                                                                                                  
 ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                  >  2:487149/1‑100 (MQ=255)
 ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                  >  1:1285497/1‑100 (MQ=255)
  gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                                 >  1:164226/1‑100 (MQ=255)
  gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCAACCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                                 >  1:311328/1‑100 (MQ=255)
   tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                >  1:1797508/1‑100 (MQ=255)
   tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                <  1:1158603/100‑1 (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                                                                                         <  1:1323814/100‑3 (MQ=255)
          tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                         <  2:1606248/100‑3 (MQ=255)
          tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                         <  2:1993001/100‑3 (MQ=255)
           gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatca                                                                                        <  1:320134/100‑4 (MQ=255)
            aaGATCTGGTGTATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcat                                                                                       <  1:328627/100‑5 (MQ=38)
            aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                       <  1:1292495/100‑1 (MQ=255)
            aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                       >  1:315353/1‑100 (MQ=255)
            aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                       <  1:1091745/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:1733893/100‑1 (MQ=255)
                      gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                             <  1:621539/100‑1 (MQ=255)
                      gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGAAAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                             <  2:1389358/100‑1 (MQ=255)
                       cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                            <  2:619740/100‑1 (MQ=255)
                        aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                           <  2:417514/100‑1 (MQ=255)
                        aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                           >  1:1335852/1‑100 (MQ=255)
                         ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                            <  1:1559462/98‑1 (MQ=255)
                         ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                          >  2:39827/1‑100 (MQ=255)
                         ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                          <  1:72689/100‑1 (MQ=255)
                         ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                          <  1:1638308/100‑1 (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)
                                     tGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                              <  2:1998167/100‑1 (MQ=255)
                                     tGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcc                                                              >  2:416090/1‑100 (MQ=255)
                                      gTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTccc                                                             >  1:200586/1‑100 (MQ=255)
                                        cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                           >  2:1409156/1‑100 (MQ=255)
                                        cGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                           >  1:1601960/1‑100 (MQ=255)
                                         gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                           >  2:1796258/1‑99 (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                                                      <  2:1285490/100‑1 (MQ=255)
                                             cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                      >  1:1552487/1‑100 (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                                        >  1:1419370/1‑100 (MQ=255)
                                                           cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                        >  2:1647441/1‑100 (MQ=255)
                                                              ttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttttt                                     >  1:2137066/1‑100 (MQ=255)
                                                                gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttttac                                   >  2:1494904/1‑98 (MQ=255)
                                                                gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                   <  2:2305055/100‑1 (MQ=255)
                                                                gAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                   <  1:1536561/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                   >  2:898770/1‑100 (MQ=255)
                                                                                cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                   >  1:1236689/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:1033489/100‑1 (MQ=255)
                                                                                              aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCa     >  1:194773/1‑100 (MQ=255)
                                                                                               ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg    <  1:1117849/100‑1 (MQ=255)
                                                                                               ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg    >  1:681525/1‑100 (MQ=255)
                                                                                               ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg    <  2:2213578/100‑1 (MQ=255)
                                                                                                 aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAggc  >  1:1861813/1‑99 (MQ=255)
                                                                                                 aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAgg   >  1:4527/1‑99 (MQ=255)
                                                                                                 aCTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt  >  1:1313698/1‑100 (MQ=255)
                                                                                                  cTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg    >  2:792053/1‑97 (MQ=255)
                                                                                                  |                                                                                                  
GGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGT  >  NC_000913/2364900‑2365096

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: