Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,032 (G)6→7 coding (1002/1047 nt) wzzE → Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,0261.G98.9% 323.1 / ‑3.4 89coding (996/1047 nt)wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (43/45);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 103 total aligned reads displayed.

TCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAG  >  NC_000913/3969927‑3970125
                                                                                                    |                                                                                                    
tCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTc                                                                                                       >  2:2992883/1‑100 (MQ=255)
 ccGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑g                                                                                                    >  2:171372/1‑100 (MQ=255)
 ccGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑g                                                                                                    >  2:673595/1‑100 (MQ=255)
  cgcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGCC‑‑gg                                                                                                   <  2:2649928/100‑1 (MQ=255)
     gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑ggggg                                                                                                <  2:1610887/100‑1 (MQ=255)
       tttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑ggggggg                                                                                              >  2:2813111/1‑99 (MQ=255)
        ttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggc                                                                                             <  1:1493282/100‑3 (MQ=255)
         tCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggct                                                                                            <  2:1006748/100‑4 (MQ=255)
         tCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggct                                                                                            >  2:2396201/1‑97 (MQ=255)
         tCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggct                                                                                            <  2:2044085/100‑4 (MQ=255)
         tCAGACCGATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggct                                                                                            >  1:1629356/1‑97 (MQ=255)
          cAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggctg                                                                                           <  2:1817623/100‑5 (MQ=255)
          cAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑gggggggctg                                                                                           >  1:2337096/1‑96 (MQ=255)
            gACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGAt                                                                                          >  1:267604/1‑100 (MQ=255)
             aCCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATc                                                                                         <  1:1585743/100‑1 (MQ=255)
                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCgggg                                                                                     <  2:1136528/100‑1 (MQ=255)
                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCgggg                                                                                     >  1:572469/1‑100 (MQ=255)
                  tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGc                                                                                    <  2:2283201/100‑1 (MQ=255)
                   cGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCt                                                                                   <  1:2710073/100‑1 (MQ=255)
                   cGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCt                                                                                   <  1:2992913/100‑1 (MQ=255)
                    gCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGc                                                                                    <  2:2544590/98‑1 (MQ=255)
                      tATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGgt                                                                                >  1:325095/1‑100 (MQ=255)
                       aTTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGgtg                                                                               <  2:1950479/100‑1 (MQ=255)
                        tttGCGATCGCCGGAAGAACCGGTAGACCGCGATAGCCCACGTCCGGCCTTCCTCATGATTATGTGGGGCCTCGTCGGGGGGGCTGATCGGGGCTGgtg                                                                               <  2:1388599/99‑1 (MQ=255)
                          tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCg                                                                            >  1:2041675/1‑100 (MQ=255)
                           gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGc                                                                           >  1:1533083/1‑100 (MQ=255)
                           gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGc                                                                           <  2:2172808/100‑1 (MQ=255)
                           gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGc                                                                           <  1:375442/100‑1 (MQ=255)
                                cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTaa                                                                      >  1:328458/1‑100 (MQ=255)
                                  ccGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAAcc                                                                    <  1:2055194/100‑1 (MQ=255)
                                  ccGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAAcc                                                                    >  1:1460608/1‑100 (MQ=255)
                                    ggAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccg                                                                  >  1:2491173/1‑100 (MQ=255)
                                    ggAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccg                                                                  >  1:2514291/1‑100 (MQ=255)
                                       agaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTCACccgccg                                                               <  2:1490756/100‑1 (MQ=255)
                                       agaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                               >  2:823330/1‑100 (MQ=255)
                                        gaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGt                                                              >  2:1749541/1‑100 (MQ=255)
                                         aaCCGGTAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                               <  2:1584007/98‑1 (MQ=255)
                                         aaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGtt                                                             <  1:394410/100‑1 (MQ=255)
                                          aCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                            >  1:2439293/1‑100 (MQ=255)
                                           ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                           <  2:1596683/100‑1 (MQ=255)
                                           ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                           >  1:2682452/1‑100 (MQ=255)
                                            cGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                          <  1:1749560/100‑1 (MQ=255)
                                            cGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                          >  1:236187/1‑100 (MQ=255)
                                             ggTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                         >  1:2071361/1‑100 (MQ=255)
                                             ggTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                         <  1:2376430/100‑1 (MQ=255)
                                              gTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCTGCCCTTGCTCg                                                        <  2:1629338/100‑1 (MQ=255)
                                              gTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCg                                                        <  1:1460629/100‑1 (MQ=255)
                                              gTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCg                                                        <  2:1740696/100‑1 (MQ=255)
                                              gTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                       >  2:438191/1‑100 (MQ=255)
                                                aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaa                                                      >  2:1110001/1‑100 (MQ=255)
                                                aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaa                                                      >  1:791317/1‑100 (MQ=255)
                                                 aaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa                                                     <  1:171373/100‑1 (MQ=255)
                                                 aaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa                                                     >  1:1138463/1‑100 (MQ=255)
                                                  aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAt                                                    <  2:1104643/100‑1 (MQ=255)
                                                  aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAt                                                    <  1:2391766/100‑1 (MQ=255)
                                                      cgATGGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCa                                                <  2:1460594/100‑1 (MQ=255)
                                                      cgATAGCCCACGTCGTGCCTTCATGATGATTCTGTGGCGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCa                                                <  2:2075434/100‑1 (MQ=255)
                                                        aTAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAac                                              <  1:937378/100‑1 (MQ=255)
                                                           gCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACAct                                           >  1:2813243/1‑100 (MQ=255)
                                                           gCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACAct                                           <  2:2405410/100‑1 (MQ=255)
                                                            cccACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACAct                                           >  2:2855350/1‑100 (MQ=255)
                                                              cACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACActgct                                        >  1:2824971/1‑100 (MQ=255)
                                                                  tcgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCgg                                    <  1:1206466/100‑1 (MQ=255)
                                                                  tcgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCgg                                    <  1:2663431/100‑1 (MQ=255)
                                                                  tcgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCgg                                    <  1:2692369/100‑1 (MQ=255)
                                                                    gtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACActgctatttcg                                  <  1:93711/100‑7 (MQ=255)
                                                                    gtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTg                                  >  2:76444/1‑100 (MQ=255)
                                                                    gtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTg                                  <  1:188716/100‑1 (MQ=255)
                                                                     tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGa                                 >  1:2655396/1‑100 (MQ=255)
                                                                     tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGa                                 >  2:2633906/1‑100 (MQ=255)
                                                                       ccTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgc                               <  1:1542131/100‑1 (MQ=255)
                                                                       ccTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgc                               >  2:316413/1‑100 (MQ=255)
                                                                        cTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgcg                              <  1:684578/100‑1 (MQ=255)
                                                                        cTTCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgcg                              <  2:267603/100‑1 (MQ=255)
                                                                         ttCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgcgc                             <  2:2414221/100‑1 (MQ=255)
                                                                          tCCTGATGATTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCa                            >  2:674787/1‑100 (MQ=255)
                                                                             tgatgaTTATGTGGGGCATTGTCGGGGGGGCTGGTCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGAGAGCGCAAAg                         >  2:2021085/1‑100 (MQ=255)
                                                                              gatgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAgg                        >  1:628218/1‑100 (MQ=255)
                                                                              gatgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAgg                        <  2:1533067/100‑1 (MQ=255)
                                                                              gatgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAgg                        <  2:2491150/100‑1 (MQ=255)
                                                                               atgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgc                       <  2:1269004/100‑1 (MQ=255)
                                                                               atgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgc                       >  2:461763/1‑100 (MQ=255)
                                                                               atgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGACTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgc                       >  2:1655961/1‑100 (MQ=255)
                                                                                tgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgcg                      >  2:545231/1‑100 (MQ=255)
                                                                                tgatTATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgcg                      >  1:1828059/1‑100 (MQ=255)
                                                                                    tATGTGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcg                  <  2:1289640/100‑1 (MQ=255)
                                                                                       gtgGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg               >  2:2222108/1‑100 (MQ=255)
                                                                                       gtgGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg               >  1:2295920/1‑100 (MQ=255)
                                                                                        tgGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccgc              <  2:2071337/100‑1 (MQ=255)
                                                                                         ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt             >  2:2907643/1‑100 (MQ=255)
                                                                                         ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt             >  2:2907586/1‑100 (MQ=255)
                                                                                         ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt             <  2:3072022/100‑1 (MQ=255)
                                                                                         ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt             >  1:69786/1‑100 (MQ=255)
                                                                                         ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTCCGGTGAGCGCCAAGGCGCTCGCCGCt             >  2:282217/1‑100 (MQ=255)
                                                                                            gCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt          >  2:1769510/1‑100 (MQ=255)
                                                                                            gCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt          >  1:1739064/1‑100 (MQ=255)
                                                                                              aTTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTc        <  2:1413080/100‑1 (MQ=255)
                                                                                               ttGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg       <  2:3019721/100‑1 (MQ=255)
                                                                                               ttGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg       >  1:1707976/1‑100 (MQ=255)
                                                                                                   c‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAagag  <  2:758036/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
TCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAG  >  NC_000913/3969927‑3970125

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