Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,553,210 A→G 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,2100AG99.0% 346.6 / ‑4.4 97K287E (AAA→GAA) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (37/58);  total (38/59)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

TAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGA  >  NC_000913/2553115‑2553304
                                                                                               |                                                                                              
tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                            >  1:251605/1‑100 (MQ=255)
tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                            >  1:893675/1‑100 (MQ=255)
tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                            <  2:644142/100‑1 (MQ=255)
 aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAg                                                                                           >  2:990188/1‑100 (MQ=255)
 aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAg                                                                                           >  2:2484067/1‑100 (MQ=255)
  ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                          <  1:212369/100‑1 (MQ=255)
  ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                          <  2:1390661/100‑1 (MQ=255)
  ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGc                                                                                          >  2:1575562/1‑100 (MQ=255)
    cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGa                                                                                        <  2:1677109/100‑1 (MQ=255)
      cggcggCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATa                                                                                      <  2:2451668/100‑1 (MQ=255)
            cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTAt                                                                                <  2:532412/100‑1 (MQ=255)
             gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATg                                                                               <  2:612002/100‑1 (MQ=255)
               ttCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGaa                                                                             <  1:586714/100‑1 (MQ=255)
                 cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                           >  1:1247128/1‑100 (MQ=255)
                 cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                           <  1:2484047/100‑1 (MQ=255)
                 cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                           <  1:2491806/100‑1 (MQ=255)
                 cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                           <  1:2551251/100‑1 (MQ=255)
                  gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcc                                                                          >  1:1151245/1‑100 (MQ=255)
                   tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAACAGCGATAAGTTATGAAAccc                                                                         <  1:1725839/100‑1 (MQ=255)
                   tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcgc                                                                         >  1:1710492/1‑98 (MQ=255)
                   tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAccc                                                                         <  1:2695956/100‑1 (MQ=255)
                    cAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCg                                                                        <  2:987542/100‑1 (MQ=255)
                    cAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCg                                                                        <  1:1916336/100‑1 (MQ=255)
                     aGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGa                                                                       >  1:1606892/1‑100 (MQ=255)
                      gagaaTCGCCACAGCGATTGCGGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGAGAAGTTATGAAACCCGAc                                                                      >  1:835228/4‑100 (MQ=255)
                      gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                      >  1:306241/1‑100 (MQ=255)
                      gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                      <  2:682106/100‑1 (MQ=255)
                      gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                      <  2:1594045/100‑1 (MQ=255)
                      gAAAATCGCCACAGCGAGTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                      <  2:114460/100‑1 (MQ=255)
                       aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACg                                                                     <  2:1738136/100‑1 (MQ=255)
                           tCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcaca                                                                 <  2:642131/100‑1 (MQ=255)
                            cGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcaca                                                                 <  2:939221/99‑1 (MQ=255)
                               cacaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAg                                                             <  1:35688/100‑1 (MQ=255)
                                 cttccgatcttGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                           <  1:1087917/90‑1 (MQ=255)
                                 caGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                           <  2:100737/100‑1 (MQ=255)
                                 caGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                           <  1:2583446/100‑1 (MQ=255)
                                 caGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                           >  1:1975808/1‑100 (MQ=255)
                                  aGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGtt                                                          <  2:1247136/100‑1 (MQ=255)
                                   gCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGGTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                         >  2:1997036/1‑100 (MQ=255)
                                     gATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaaa                                                       >  2:553618/1‑100 (MQ=255)
                                           tGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGtagatcggaag                                                 >  2:1087925/1‑90 (MQ=255)
                                            gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                <  1:2440094/100‑1 (MQ=255)
                                            gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                <  1:2590539/100‑1 (MQ=255)
                                            gAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTc                                                >  1:1202771/1‑100 (MQ=255)
                                            caaGGCGTGATCATTTAGTTCCACGGGTTCTACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGttt                                                 <  2:518217/98‑1 (MQ=255)
                                               ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTTTGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                             <  2:1990487/100‑1 (MQ=255)
                                               ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                             <  1:1100510/100‑1 (MQ=255)
                                                gCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                            <  1:1821716/100‑1 (MQ=255)
                                                 cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                           <  2:251607/100‑1 (MQ=255)
                                                  gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                          <  1:8937/100‑1 (MQ=255)
                                                  gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                          >  1:958121/1‑100 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                        >  1:1830446/1‑100 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                        >  2:2674079/1‑100 (MQ=255)
                                                    gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                        >  1:1271298/1‑100 (MQ=255)
                                                     aTCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAcn                                       <  1:2478835/100‑2 (MQ=255)
                                                      tCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCt                                      <  2:1482069/100‑1 (MQ=255)
                                                      tCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCt                                      >  2:1660508/1‑100 (MQ=255)
                                                       cAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCtt                                     >  1:144785/1‑100 (MQ=255)
                                                       cAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCtt                                     >  2:1820758/1‑100 (MQ=255)
                                                        aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTc                                    <  1:9046/100‑1 (MQ=255)
                                                        aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTc                                    <  1:856135/100‑1 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                   >  2:2791801/1‑100 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                   <  1:233203/100‑1 (MQ=255)
                                                         gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                   >  2:2064633/1‑100 (MQ=255)
                                                           tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                 >  2:2089080/1‑100 (MQ=255)
                                                           tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                 >  2:2086930/1‑100 (MQ=255)
                                                            aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                <  2:887090/100‑1 (MQ=255)
                                                                ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                            >  1:2485036/1‑100 (MQ=255)
                                                                ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                            >  1:170247/1‑100 (MQ=255)
                                                                 cACTGGTTCGACAACCAGAAAGCACATTCGCAATAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGa                           <  1:574871/100‑1 (MQ=255)
                                                                  aCTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAt                          >  1:2343878/1‑100 (MQ=255)
                                                                  aCTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAt                          <  2:1357244/100‑1 (MQ=255)
                                                                  aCTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAt                          <  1:1848201/100‑1 (MQ=255)
                                                                   cTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAtt                         <  2:269609/100‑1 (MQ=255)
                                                                    tGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                        <  1:1269202/100‑1 (MQ=255)
                                                                    tGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                        <  1:1554158/100‑1 (MQ=255)
                                                                    tGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                        >  1:2159849/1‑100 (MQ=255)
                                                                     ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                       >  1:615494/1‑100 (MQ=255)
                                                                     ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                       >  2:2206314/1‑100 (MQ=255)
                                                                     ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                       <  1:2273970/100‑1 (MQ=255)
                                                                       ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                     <  1:990181/100‑1 (MQ=255)
                                                                       ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                     <  1:492729/100‑1 (MQ=255)
                                                                        tCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCGACCTTCAGGATACGATTTGTca                    <  1:1637925/100‑1 (MQ=255)
                                                                        tCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTca                    <  2:1348538/100‑1 (MQ=255)
                                                                         cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                   >  1:194871/1‑100 (MQ=255)
                                                                         cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcaa                   >  2:1637937/1‑99 (MQ=255)
                                                                         cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACACACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                   >  1:90881/1‑100 (MQ=255)
                                                                          gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                  <  1:2719914/100‑1 (MQ=255)
                                                                          gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                  <  2:2212294/100‑1 (MQ=255)
                                                                            cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                <  2:1202779/100‑1 (MQ=255)
                                                                            cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                <  2:448589/100‑1 (MQ=255)
                                                                            cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                <  1:2206298/100‑1 (MQ=255)
                                                                              aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt              >  1:505660/1‑100 (MQ=255)
                                                                              aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt              >  1:1918484/1‑100 (MQ=255)
                                                                                 aGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTggc           <  2:1918498/100‑3 (MQ=255)
                                                                                    aaGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgc        <  1:2378634/100‑1 (MQ=255)
                                                                                       cacaTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGt     >  1:1049331/1‑100 (MQ=255)
                                                                                        acaTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTc    <  1:1842456/100‑1 (MQ=255)
                                                                                          aTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGa  >  2:2704260/1‑100 (MQ=255)
                                                                                          aTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGa  <  1:1660496/100‑1 (MQ=255)
                                                                                               |                                                                                              
TAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGA  >  NC_000913/2553115‑2553304

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