Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,319,581 A→G 16.9% V99A (GTG→GCG) ‡ rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,319,5810AG16.9% 193.0 / 32.0 89V99A (GTG→GCG) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base A (42/32);  new base G (5/10);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGG  >  NC_000913/2319487‑2319680
                                                                                              |                                                                                                   
gcgcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATcc                                                                                                <  2:301299/100‑1 (MQ=255)
  gcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGc                                                                                              >  1:1932118/1‑100 (MQ=255)
   cATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGTTGTTACTCATTGCGGAACAATCGGAGTCGCCAAACAGCGGTTCAAACGCAGGCACATCCGcc                                                                                             >  2:892648/1‑100 (MQ=255)
    aTAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCt                                                                                            <  1:644056/100‑1 (MQ=255)
     tactccACGCCAACGCCTCCATAGAACCTCGCGAGGAGTTACTCATTGCGCAACAACCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGcc                                                                                             <  2:2178552/94‑1 (MQ=255)
       aaCCACGCCAGTGACTCCAGAGAGCCTTGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGc                                                                                          <  2:1935216/99‑1 (MQ=255)
         ccACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGtt                                                                                       <  1:1977610/100‑1 (MQ=255)
         ccACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGtt                                                                                       >  2:2403458/1‑100 (MQ=255)
         ccACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACCCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGtt                                                                                       <  2:1143195/100‑1 (MQ=255)
          cACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTGAAACGCAGGCGGATGTGCCTGCGttt                                                                                      <  1:1556691/100‑1 (MQ=38)
          cACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                      <  2:1755845/100‑1 (MQ=38)
          cACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                      <  1:322101/100‑1 (MQ=38)
          cACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                      <  2:441687/100‑1 (MQ=38)
          cACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                      <  1:645106/100‑1 (MQ=38)
          cACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                      <  1:66863/100‑1 (MQ=38)
           aCGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtc                                                                                     <  1:1863141/100‑1 (MQ=255)
           aCGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtc                                                                                     <  1:1860878/100‑1 (MQ=255)
            cGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCTTTtct                                                                                    >  2:2356587/1‑100 (MQ=255)
            cGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtct                                                                                    <  2:1852566/100‑1 (MQ=255)
             gCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGtttgaa                                                                                   <  2:97184/100‑4 (MQ=38)
              ccAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGtttgaac                                                                                  <  1:1015342/100‑5 (MQ=37)
              ccAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtctc                                                                                   >  1:2105782/1‑99 (MQ=255)
              ccAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                                                                                  <  2:2383253/100‑1 (MQ=255)
                  tGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTcaa                                                                              >  2:1551813/1‑98 (MQ=255)
                  tGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCa                                                                              >  1:1862566/1‑100 (MQ=255)
                     cTCCAGAGACCCTTGCCAGGTGTTACTCATTGCGGAACAGTCGGAGTCGGCAAACAGCGTTTCAAACGCAGGTACATCCGCCTGTGTTTCTCGTCCAcg                                                                            >  2:171706/1‑99 (MQ=255)
                     cTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc                                                                           >  2:541361/1‑100 (MQ=255)
                     cTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgc                                                                           >  1:1509031/1‑100 (MQ=255)
                       ccAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgg                                                                         >  1:1651416/1‑100 (MQ=255)
                         agagAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgggg                                                                       >  1:2552266/1‑100 (MQ=255)
                          gagaACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGa                                                                      <  1:1528571/100‑1 (MQ=255)
                              aCCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGc                                                                  <  2:1861061/100‑1 (MQ=255)
                               ccTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGca                                                                 <  2:665287/100‑1 (MQ=255)
                               ccTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGca                                                                 >  1:1618080/1‑100 (MQ=255)
                                cTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGcac                                                                >  1:2001515/1‑100 (MQ=255)
                                cTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGcac                                                                >  1:415437/1‑100 (MQ=255)
                                  cGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCacac                                                              <  2:939307/100‑1 (MQ=255)
                                  cGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCacac                                                              <  2:1932104/100‑1 (MQ=255)
                                   gCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCacac                                                              <  2:126957/99‑1 (MQ=255)
                                   gCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACAcc                                                             >  1:1500895/1‑100 (MQ=255)
                                    ccAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCg                                                            >  1:270313/1‑100 (MQ=255)
                                       ggTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAACGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGttt                                                         >  2:1579578/1‑100 (MQ=255)
                                          gtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGCAAGCACACCGTTTTCt                                                      >  2:2275019/1‑100 (MQ=255)
                                          gtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCt                                                      >  2:2164965/1‑100 (MQ=255)
                                              cTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCg                                                  <  1:533844/100‑1 (MQ=255)
                                              cTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCg                                                  >  1:2010318/1‑100 (MQ=255)
                                               tCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGa                                                 <  1:1540238/100‑1 (MQ=255)
                                                cATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAt                                                >  1:271764/1‑100 (MQ=255)
                                                    gaggAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGAAAAGCACACCGTTTTCTGCCGATAAGc                                            >  1:1726476/3‑100 (MQ=255)
                                                    gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGc                                            >  1:1074152/1‑100 (MQ=255)
                                                    gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGc                                            >  1:1280691/1‑100 (MQ=255)
                                                    gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGc                                            >  1:463326/1‑100 (MQ=255)
                                                     cGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCg                                           <  2:566530/100‑1 (MQ=255)
                                                     cGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCg                                           <  1:319325/100‑1 (MQ=255)
                                                               ggAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCCCGAGGTGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCg                                 <  2:853469/100‑1 (MQ=255)
                                                               ggAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCg                                 <  1:2231550/100‑1 (MQ=255)
                                                               ggAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCg                                 >  1:459528/1‑100 (MQ=255)
                                                                  gTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATg                              <  1:743173/100‑1 (MQ=255)
                                                                   tCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAGGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATg                              <  2:1144164/99‑1 (MQ=255)
                                                                    cGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTa                            >  1:22186/1‑100 (MQ=255)
                                                                    cGGCAAACAGCGGTTCAAACGCAGGCACATACGCCTGCGTTTCTCGTCAACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTa                            >  2:332739/1‑100 (MQ=255)
                                                                     ggCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTAc                           >  1:1221929/1‑100 (MQ=255)
                                                                     ggCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTAc                           >  2:5657/1‑100 (MQ=255)
                                                                     ggCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTAc                           >  1:814487/1‑100 (MQ=255)
                                                                       cAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACtt                         >  2:154176/1‑100 (MQ=255)
                                                                        aaaCAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTc                        <  2:1594382/100‑1 (MQ=255)
                                                                        aaaCAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTc                        >  1:48463/1‑100 (MQ=255)
                                                                          agaGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAg                      <  2:1554309/98‑1 (MQ=255)
                                                                           cAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGc                     >  2:588981/1‑100 (MQ=255)
                                                                             gCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGctc                   <  2:463322/100‑1 (MQ=255)
                                                                              gacgagaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGctct                  >  1:2116691/7‑100 (MQ=37)
                                                                               ggTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCtt                 <  1:2185204/100‑1 (MQ=255)
                                                                               ggTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCtt                 >  1:687134/1‑100 (MQ=255)
                                                                                gTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCttt                >  1:1441130/1‑100 (MQ=255)
                                                                                cttCAAACGCAGGCACATACGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGTCGATAAGCGATTTTCGGCGATGTACTTCAGCTCtt                 >  2:1244445/2‑99 (MQ=255)
                                                                                cgagaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCCGCTCttt                >  2:332967/5‑100 (MQ=37)
                                                                                 tctaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGTGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtc               <  2:1795686/97‑1 (MQ=38)
                                                                                  tCAAACGCAGGCACATCCGCCTGCTTTTCTCGTCCACGCTTGGCAAGCACACCGTTTTCTGCCGATAAGCTATTTTCGGCGATGTACTTCAGCTCTTtca              <  2:1817691/100‑1 (MQ=255)
                                                                                    aaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAGAAGCGATTTTCGGCGATGTACTTCAGCTCTTtcaga            >  1:1795700/1‑98 (MQ=38)
                                                                                    aaaCGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtcatc            <  2:2411689/100‑1 (MQ=255)
                                                                                    aaaCGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtcatc            <  1:2023894/100‑1 (MQ=255)
                                                                                     aaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACTGTTTTCTGCCGATAAGCGATTTTAGGCGATGTACTTCAGCTCTTtcacca           >  1:2301858/1‑100 (MQ=255)
                                                                                     aaCGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtcatca           >  1:1468842/1‑100 (MQ=255)
                                                                                      aCGCAGGCGGATGTGACTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcac          >  2:481272/1‑100 (MQ=255)
                                                                                       cGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcaca         >  1:1368114/1‑100 (MQ=255)
                                                                                       cGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGGTTTCTGCCGATAAGCGATTTTCGGCGAGGTACTTCAGCTCTTTcatcaca         >  2:630462/1‑100 (MQ=255)
                                                                                         cAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcacatc       <  2:825775/100‑1 (MQ=255)
                                                                                         cAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcacatc       <  1:661520/100‑1 (MQ=255)
                                                                                           ggCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCag     >  2:795490/1‑100 (MQ=255)
                                                                                             cacaTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCagag   >  1:1269697/1‑100 (MQ=255)
                                                                                             cacaTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCagag   >  1:1163359/1‑100 (MQ=255)
                                                                                              acaTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAgg  >  1:1399083/1‑100 (MQ=255)
                                                                                              acaTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAgg  >  1:1400541/1‑100 (MQ=255)
                                                                                              |                                                                                                   
GCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGG  >  NC_000913/2319487‑2319680

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: