Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,210 A→G 36.2% K287E (AAA→GAA)  amiA → N‑acetylmuramoyl‑l‑alanine amidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,2100AG36.2% 118.7 / 128.1 127K287E (AAA→GAA) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (41/40);  new base G (26/20);  total (67/60)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 132 total aligned reads displayed.

TGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGCG  >  NC_000913/2553112‑2553309
                                                                                                  |                                                                                                   
tGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGaa                                                                                                    <  1:1940256/100‑1 (MQ=255)
 gTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGaaa                                                                                                   <  2:975514/100‑1 (MQ=255)
  ttAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaa                                                                                                  >  1:1255113/1‑100 (MQ=255)
    aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGAATGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAg                                                                                                >  1:2125065/1‑100 (MQ=255)
     ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                               >  1:876366/1‑100 (MQ=255)
     ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                               <  2:876373/100‑1 (MQ=255)
     ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGc                                                                                               >  2:2046014/1‑100 (MQ=255)
       cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                             >  1:414938/1‑100 (MQ=255)
       cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGa                                                                                             >  2:2805748/1‑100 (MQ=255)
        aCGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                                                            <  1:2072652/100‑1 (MQ=255)
        aCGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt                                                                                            <  2:968629/100‑1 (MQ=255)
            cggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGt                                                                                        <  1:242996/100‑1 (MQ=255)
               cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTAt                                                                                     <  2:104202/100‑1 (MQ=255)
               cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAAAAAATAAAAAAAAGCGATAAGTTAt                                                                                     <  1:2593910/100‑1 (MQ=255)
                gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATg                                                                                    >  2:2977005/1‑100 (MQ=255)
                gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATg                                                                                    >  2:50636/1‑100 (MQ=255)
                    cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACGTGTTCGACAACCAGAAAAAAAAAAAAAAAAAGCGATAAGTTATGAAAc                                                                                <  1:1131987/100‑1 (MQ=255)
                     gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAcc                                                                               >  2:2048054/1‑100 (MQ=255)
                     gTCAGAAAATCGCCACAGCGATTGCTCAAGGCGTGATCAGGTATTTCCACTGGTTCGACAACAAGAAAGCACATTCGAAAAAGCGATAGGTTATGaa                                                                                  >  1:3255079/1‑97 (MQ=255)
                       cAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCg                                                                             >  2:1178239/1‑100 (MQ=255)
                       cAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                                              <  2:2758466/99‑1 (MQ=255)
                          aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg                                                                          <  2:2859916/100‑1 (MQ=255)
                                gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                    >  2:2118336/1‑100 (MQ=255)
                                gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                    <  2:2717819/100‑1 (MQ=255)
                                 ccACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCa                                                                   >  2:1576925/1‑100 (MQ=255)
                                  cacaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAg                                                                  >  1:3395654/1‑100 (MQ=255)
                                    caGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                                <  1:105293/100‑1 (MQ=255)
                                    caGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                                <  2:1148233/100‑1 (MQ=255)
                                      ttccgatcTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAAAAAAAAAAAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                              <  1:2664904/93‑1 (MQ=255)
                                      gCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                              <  2:1394127/100‑1 (MQ=255)
                                      gCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                              >  1:924770/1‑100 (MQ=255)
                                      gCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                              <  1:2549282/100‑1 (MQ=255)
                                         aTTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAAc                                                           <  1:1908770/100‑1 (MQ=255)
                                          ttGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACa                                                          <  2:963584/100‑1 (MQ=255)
                                             cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                       >  2:736381/1‑100 (MQ=255)
                                             cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                       <  1:303487/100‑1 (MQ=255)
                                             cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                       >  1:648695/1‑100 (MQ=255)
                                             cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                       >  1:2317694/1‑100 (MQ=255)
                                             cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                       >  1:3016317/1‑100 (MQ=255)
                                               gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGGACACCAGGTTAAACAGTTTc                                                     <  2:3251788/100‑1 (MQ=255)
                                               gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                     >  1:839482/1‑100 (MQ=255)
                                               gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaagatcggaa                                                     >  2:2664924/1‑92 (MQ=255)
                                               gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                     <  2:2425757/100‑1 (MQ=255)
                                               gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                     >  1:948162/1‑100 (MQ=255)
                                                aaGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCt                                                    <  2:3317453/100‑1 (MQ=255)
                                                 aGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTg                                                   >  1:3165024/1‑100 (MQ=255)
                                                    cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                                >  1:3275203/1‑100 (MQ=255)
                                                    cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                                >  2:915243/1‑100 (MQ=255)
                                                    cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                                >  2:3435240/1‑100 (MQ=255)
                                                     gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                               >  1:1206876/1‑100 (MQ=255)
                                                      tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCaa                                              >  1:2531806/1‑100 (MQ=255)
                                                      tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCaa                                              <  1:673565/100‑1 (MQ=255)
                                                       gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                             <  2:2691953/100‑1 (MQ=255)
                                                         tCAGTTAGTTCTACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCTACGCACACCAGGTTAAACAGTTTCTGCTCAACCt                                           >  1:1817244/1‑100 (MQ=255)
                                                             ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAg                                       <  1:2980210/100‑1 (MQ=255)
                                                              tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTAAgg                                      >  2:841393/1‑100 (MQ=255)
                                                               aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                     <  2:1076966/100‑1 (MQ=255)
                                                                  tCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATAc                                  >  2:3419139/1‑100 (MQ=255)
                                                                   ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                                 >  1:186446/1‑100 (MQ=255)
                                                                   ccACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                                 <  1:935309/100‑1 (MQ=255)
                                                                      cTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAtt                              >  2:2016893/1‑100 (MQ=255)
                                                                       tGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                             >  2:499139/1‑100 (MQ=255)
                                                                        ggTTCGACAATCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTACACAGTTTCTGCTCAACCTTCAGGATACGATTTg                            <  1:2273594/100‑1 (MQ=255)
                                                                        ggTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                            <  2:2531149/100‑1 (MQ=255)
                                                                         gTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGt                           <  2:2476843/100‑1 (MQ=255)
                                                                          ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                          <  1:1576910/100‑1 (MQ=255)
                                                                          ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                          <  1:1360054/100‑1 (MQ=255)
                                                                          ttCGACAACCAGAAAGCACATTCGAAAACGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                          >  1:2699042/1‑100 (MQ=255)
                                                                          ttCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                          <  2:1484850/100‑1 (MQ=255)
                                                                          ttCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                          <  1:250877/100‑1 (MQ=255)
                                                                          ttCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                          <  1:412355/100‑1 (MQ=255)
                                                                           tCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTca                         <  1:1020839/100‑1 (MQ=255)
                                                                             gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                       >  1:1486041/1‑100 (MQ=255)
                                                                             gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                       >  1:1023198/1‑100 (MQ=255)
                                                                             gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                       <  2:1271561/100‑1 (MQ=255)
                                                                             gACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                       >  1:1748057/1‑100 (MQ=255)
                                                                             gACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                       >  2:1770156/1‑100 (MQ=255)
                                                                             gACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                       <  2:3189879/100‑1 (MQ=255)
                                                                               cAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                     >  1:3099681/1‑100 (MQ=255)
                                                                                aaCCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCagcagc                    <  1:3141799/100‑1 (MQ=255)
                                                                                 aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt                   >  1:571600/1‑100 (MQ=255)
                                                                                    aGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc                <  2:2167937/100‑1 (MQ=255)
                                                                                    aGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc                <  1:1802796/100‑1 (MQ=255)
                                                                                    aGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc                <  2:2956608/100‑1 (MQ=255)
                                                                                    aGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc                >  2:1070206/1‑100 (MQ=255)
                                                                                    aGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc                <  2:3212007/100‑1 (MQ=255)
                                                                                       aaGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgc             <  1:2756836/100‑1 (MQ=255)
                                                                                           acaTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTc         <  1:1445896/100‑1 (MQ=255)
                                                                                           acaTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTc         <  2:2221642/100‑1 (MQ=255)
                                                                                            caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGATGACAGCCGTCg        >  2:1148669/1‑100 (MQ=255)
                                                                                            caTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg        >  1:1906840/1‑100 (MQ=255)
                                                                                            caTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg        >  1:169076/1‑100 (MQ=255)
                                                                                             aTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGa       >  1:829181/1‑100 (MQ=255)
                                                                                              ttCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAt      >  1:2050728/1‑100 (MQ=255)
                                                                                              ttCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAt      >  2:2684287/1‑100 (MQ=255)
                                                                                               tCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATg     >  2:1282419/1‑100 (MQ=255)
                                                                                                 gAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGCTTTGTCAGCAGCTGACAGCCGTCGATGgc   >  2:3171072/1‑100 (MQ=255)
                                                                                                 gAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGgc   >  1:943700/1‑100 (MQ=255)
                                                                                                 gAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGgc   >  1:1222562/1‑100 (MQ=255)
                                                                                                 gAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGgc   >  1:480762/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGCG  >  NC_000913/2553112‑2553309

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: