Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,979,102 G→A 17.1% A168T (GCG→ACG)  wzyE → putative ECA polysaccharide chain elongation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,979,1020GA17.1% 318.3 / 53.2 142A168T (GCG→ACG) wzyEputative ECA polysaccharide chain elongation protein
Reads supporting (aligned to +/- strand):  ref base G (65/50);  new base A (15/9);  total (82/60)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

Only 100 of 143 total aligned reads displayed.

ACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGATTGTCG  >  NC_000913/3979004‑3979201
                                                                                                  |                                                                                                   
aCGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGc                                                                                                    >  1:1804194/1‑100 (MQ=255)
    cTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATg                                                                                                >  2:1591962/1‑100 (MQ=255)
         tGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCtggt                                                                                           >  2:947417/1‑100 (MQ=255)
           cTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCtggtgg                                                                                         <  2:1337756/100‑1 (MQ=255)
           cTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGACGATGCtggtgg                                                                                         >  2:689448/1‑100 (MQ=255)
           cTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGACGATGCtggtgg                                                                                         >  1:1590493/1‑100 (MQ=255)
             gTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTc                                                                                       <  1:925996/100‑1 (MQ=255)
             gTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTc                                                                                       >  2:569429/1‑100 (MQ=255)
             gTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTc                                                                                       <  1:1023292/100‑1 (MQ=255)
             gTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTc                                                                                       >  1:1101547/1‑100 (MQ=255)
                ccGGCTTAACTCCTACAGACAAAACTTTTCCACTGAAGTCTCCGGCGTAACGATAAAACGCTTCATTTACTTTTTCCTCCCGGCGATGCTGACGCTCTAc                                                                                    >  2:924435/1‑100 (MQ=255)
                    cTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACtt                                                                                  >  2:1864648/1‑98 (MQ=255)
                    cTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTc                                                                                >  1:1377668/1‑100 (MQ=255)
                      tAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTg                                                                              >  1:883099/1‑100 (MQ=255)
                      tAACTCCTACAGTAAGATCTTTTCAAGTGAAGTCTCCGACGTCGCGTTAAAACGCTTCTCTTACTTTTTCATCCCGGCGATGCTGGAGCTCTACATTc                                                                                >  2:53933/1‑98 (MQ=255)
                       aaCTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTgc                                                                             >  2:2222590/1‑100 (MQ=255)
                         cTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTgcgc                                                                           >  1:1011700/1‑100 (MQ=255)
                         cTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTgcg                                                                            >  1:1889218/1‑99 (MQ=255)
                           ccTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCa                                                                         >  2:1786464/1‑100 (MQ=255)
                           ccTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGCGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCa                                                                         >  2:818055/1‑100 (MQ=255)
                            cTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAg                                                                        <  1:811537/100‑1 (MQ=255)
                            cTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAg                                                                        >  2:1902451/1‑100 (MQ=255)
                              aCAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCTCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGa                                                                      >  1:1755132/1‑100 (MQ=255)
                               cAGTCAGATCTTTTCCAGTGAAGTCGCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGAc                                                                     <  1:741032/100‑1 (MQ=255)
                                 gTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAg                                                                   >  1:2208485/1‑100 (MQ=255)
                                 gTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAg                                                                   <  2:2105593/100‑1 (MQ=255)
                                  tCAGATCTTTTCCAGTGAAGTCTCCGGCGTTGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTTCGCCAGGAc                                                                     <  1:2056635/97‑1 (MQ=255)
                                  tCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGc                                                                  <  2:2103443/100‑1 (MQ=255)
                                  tCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGc                                                                  <  1:319626/100‑1 (MQ=255)
                                   cAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGc                                                                  >  2:598652/1‑99 (MQ=255)
                                     gATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaaa                                                               >  2:1677443/1‑100 (MQ=255)
                                     caTCTTTTCCAGTGAAGTATACAACATAAAGTTAAAACGCTTCTTTTACTTTTTAATCCCAGCGATGATGGTGGAATACTTTCTGCGCCAGGTAAGCaa                                                                >  2:575189/2‑99 (MQ=255)
                                      aTCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAg                                                              <  1:859735/100‑1 (MQ=255)
                                      aTCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAg                                                              >  1:2241097/1‑100 (MQ=255)
                                        cTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCg                                                            <  2:1861856/100‑1 (MQ=255)
                                         ttttCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGt                                                           >  2:570609/1‑100 (MQ=255)
                                          tttCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTg                                                          >  1:718894/1‑100 (MQ=255)
                                            tCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGc                                                        <  2:1045118/100‑1 (MQ=255)
                                            tCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGc                                                        <  1:1323447/100‑1 (MQ=255)
                                            tCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGc                                                        <  2:2401283/100‑1 (MQ=255)
                                             ccAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCt                                                       <  1:2080466/100‑1 (MQ=255)
                                             ccAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCt                                                       >  2:232285/1‑100 (MQ=255)
                                                gTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGtt                                                    <  1:191818/100‑1 (MQ=255)
                                                gTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGtt                                                    <  2:689650/100‑1 (MQ=255)
                                                    aGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTc                                                >  1:1090692/1‑100 (MQ=255)
                                                    aGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTc                                                >  1:780777/1‑100 (MQ=255)
                                                    aGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTc                                                >  1:1735370/1‑100 (MQ=255)
                                                     gTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTcc                                               >  1:111313/1‑100 (MQ=255)
                                                      tctcCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCt                                              <  2:732953/100‑1 (MQ=255)
                                                       ctcCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTc                                             >  2:761958/1‑100 (MQ=255)
                                                       ctcCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTc                                             >  1:1875809/1‑100 (MQ=255)
                                                       ctcCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTc                                             <  2:424511/100‑1 (MQ=255)
                                                           ggtgtggcgttaaATCGCTTTTTTTACTTTTTCATCCCGGCGTTGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTCCTCGTCA                                         <  1:193488‑M1/96‑1 (MQ=255)
                                                           ggcgtggcgtTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCA                                         <  1:1677463‑M1/100‑1 (MQ=255)
                                                           ggcgtggcgtTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCA                                         <  1:192564‑M1/100‑1 (MQ=255)
                                                              gtggcgtTAAAACGCTTCTTTTACTTTTTCATCCCTTCGATGCTGGTGGTCTACTTTCTGCGCCAGGACCGCACAGCGTGGCTGTTTTTCCTCGTCACCA                                      >  2:1490751‑M1/1‑100 (MQ=255)
                                                              gtggcgtTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCA                                      >  1:1838121‑M1/1‑100 (MQ=255)
                                                               tggcgtTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAc                                     <  1:1528692/100‑1 (MQ=255)
                                                                 gcgtTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACgg                                   >  1:1496635/1‑100 (MQ=255)
                                                                  cgtTAAAACGCTTCTTTTACTTTATCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGt                                  >  2:547198/1‑100 (MQ=255)
                                                                   gtTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTc                                 <  2:1213224/100‑1 (MQ=255)
                                                                      aaaaCGCTTCTTTTACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGcc                              >  2:649312/1‑100 (MQ=255)
                                                                       aaaCGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCt                             >  1:673605/1‑100 (MQ=255)
                                                                       aaaCGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCt                             >  2:1306593/1‑100 (MQ=255)
                                                                        aaCGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCt                             <  2:816642/99‑1 (MQ=255)
                                                                         aCGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCttt                           <  2:1601402/100‑1 (MQ=255)
                                                                          cGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTg                          <  1:1041363/100‑1 (MQ=255)
                                                                            cttcttTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGc                        <  2:1562954/100‑1 (MQ=255)
                                                                            cttcttTTACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGc                        <  1:649316/100‑1 (MQ=255)
                                                                            cttctcTTACTTTCTCAGCCCGGCGATCCTGTTGCACTACATTCTGCGCCAGTACAGTAAAGCGTGGCTGTTTTTCCTCGTCAACACGGTCGCCTTTgac                        <  2:630759/100‑3 (MQ=255)
                                                                             ttcttTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCt                       <  1:2062049/100‑1 (MQ=255)
                                                                              tCTTTTACTTTTTTATTTCTTCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTt                      <  2:1084304‑M2/100‑2 (MQ=255)
                                                                              tCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTt                      <  1:2176837‑M2/100‑2 (MQ=255)
                                                                               cTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTtg                     >  1:1418488‑M2/1‑98 (MQ=255)
                                                                               cTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTtg                     >  2:60033‑M2/1‑98 (MQ=255)
                                                                               cTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTtg                     >  1:588563‑M2/1‑98 (MQ=255)
                                                                                ttttACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCATCAAGGTCGCCTTTGGCTt                      <  1:448476‑M2/98‑2 (MQ=255)
                                                                                 tttACTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTtgct                   >  2:2370352‑M2/1‑96 (MQ=255)
                                                                                   tACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTtgctga                 >  1:832686‑M2/1‑94 (MQ=255)
                                                                                   tACTTTTTCATCCCGACGATGCTGGTGCTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTtgctg                  >  2:1794012‑M2/1‑94 (MQ=255)
                                                                                    aCTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTtgctgac                <  2:1196101‑M2/100‑8 (MQ=255)
                                                                                     cTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGCCAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACt               >  1:736493/1‑100 (MQ=255)
                                                                                     cTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACt               >  2:1126664/1‑100 (MQ=255)
                                                                                     cTTTTTCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACt               >  1:600479/1‑100 (MQ=255)
                                                                                      tttttCATCCCGGCGATGCTGGTGGTCTATTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACtt              <  2:484573/100‑1 (MQ=255)
                                                                                      tttttCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAagccaaaggcgaccgtgctgacga              <  1:23702‑M1/100‑25 (MQ=255)
                                                                                       ttttCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTta             <  1:499843/100‑1 (MQ=255)
                                                                                       ttttCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTta             <  2:698550/100‑1 (MQ=255)
                                                                                        tttCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAagccaaaggcgaccgtgctgacgagg            <  2:105496‑M1/100‑27 (MQ=255)
                                                                                         ttCATCCCGGCGATGCTGTTGGTCTACTTTCTGTGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAagccaaaggcgaccgtgctgacgag             <  2:1099910‑M1/98‑26 (MQ=255)
                                                                                         ttCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTtata           >  1:977910/1‑100 (MQ=255)
                                                                                          tCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTatat          >  1:539108/1‑100 (MQ=255)
                                                                                          tCATCCCGACGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAagccaaaggcgaccgtgctgacgaggaa          <  1:1762547‑M1/100‑29 (MQ=255)
                                                                                           cATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATg         <  2:1541127/100‑1 (MQ=255)
                                                                                            aTCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGa        >  2:2115849/1‑100 (MQ=255)
                                                                                             taccGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGAt       >  1:2224876/3‑100 (MQ=255)
                                                                                             tCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGAt       >  1:22621/1‑100 (MQ=255)
                                                                                             tCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGAt       >  1:2079640/1‑100 (MQ=255)
                                                                                              cccGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGAtt      >  2:161675/1‑100 (MQ=255)
                                                                                                  gCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGATTGTcg  <  1:1074455/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
ACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCGCCTTTGGCTTGCTGACTTATATGATTGTCG  >  NC_000913/3979004‑3979201

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: