Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,210 A→G 72.5% 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,2100AG72.5% 149.1 / 63.3 98K287E (AAA→GAA) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (11/16);  new base G (34/36);  total (45/53)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 103 total aligned reads displayed.

GTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGC  >  NC_000913/2553113‑2553308
                                                                                                 |                                                                                                  
gTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGATTCGACAACCAGAAAGCACATTCGgaa                                                                                                  >  1:920984/1‑100 (MQ=255)
  tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                                >  1:725236/1‑100 (MQ=255)
  tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                                <  1:1679533/100‑1 (MQ=255)
   aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                                <  2:2499340/99‑1 (MQ=255)
   aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAg                                                                                               <  1:24081/100‑1 (MQ=255)
    ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGc                                                                                              <  2:1036063/100‑1 (MQ=255)
    ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTAGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                              >  1:829623/1‑100 (MQ=255)
     gCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                              <  2:825638/99‑1 (MQ=255)
       aCGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACACCCAGAAAGCACATTCGGAAAAGCGAt                                                                                           <  2:217901/100‑1 (MQ=255)
         ggcggcGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATaa                                                                                         <  1:365024/100‑1 (MQ=255)
         ggcggcGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATaa                                                                                         <  2:725231/100‑1 (MQ=255)
             gcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTa                                                                                     >  1:2613213/1‑100 (MQ=255)
              cgatctgTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTAt                                                                                    <  1:25089/94‑1 (MQ=255)
              cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTAt                                                                                    <  2:1405952/100‑1 (MQ=255)
              cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTAt                                                                                    >  2:94030/1‑100 (MQ=255)
               gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATGTCCACTGGTTCGACAACCAGAAAGCACATTCGGATAAGCGATAAGTTATg                                                                                   <  2:1085787/100‑1 (MQ=255)
                tttCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGa                                                                                  <  1:2200803/100‑1 (MQ=255)
                   cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATTAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                               <  1:1472229/100‑1 (MQ=255)
                    gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTAtagatcg                                                                              >  2:25089/1‑94 (MQ=255)
                     tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAccc                                                                             <  1:2086928/100‑1 (MQ=255)
                     tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                                             >  2:487458/1‑100 (MQ=255)
                     tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                                             >  2:783538/1‑100 (MQ=255)
                      cAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                                             >  2:979747/1‑99 (MQ=255)
                        gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                          >  1:2449986/1‑100 (MQ=255)
                        gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                          >  1:634804/1‑100 (MQ=255)
                        gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGAc                                                                          <  2:226066/100‑1 (MQ=255)
                         aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAATCGATAAGTTATGAAACCCGACg                                                                         <  2:609542/100‑1 (MQ=255)
                            aTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcac                                                                      >  2:2000711/1‑100 (MQ=255)
                             tCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcaca                                                                     <  1:1945733/100‑1 (MQ=255)
                              cGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCacac                                                                    >  2:1880929/1‑100 (MQ=255)
                              cGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGAacac                                                                    >  1:1154339/1‑100 (MQ=255)
                               gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                   >  2:1232191/1‑100 (MQ=255)
                                  acaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAgg                                                                <  2:2005024/100‑1 (MQ=255)
                                   caGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGAAAAGTTATGAAACACGACGCACACCAGGt                                                               >  1:491237/1‑100 (MQ=255)
                                    aGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGtt                                                              <  2:1431487/100‑1 (MQ=255)
                                    aGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGtt                                                              <  2:10160/100‑1 (MQ=255)
                                    aGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                               <  1:433248/99‑1 (MQ=255)
                                        aTTGCTGAAGGCGTGATTAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAAc                                                          <  1:1283528/100‑1 (MQ=255)
                                        aTTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAAc                                                          <  1:1881710/100‑1 (MQ=255)
                                           gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                       >  1:1959297/1‑100 (MQ=255)
                                            cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                      <  1:130053/100‑1 (MQ=255)
                                             tGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGttt                                                     >  1:1097415/1‑100 (MQ=255)
                                                 ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                                 <  2:920975/100‑1 (MQ=255)
                                                  gCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGATATAAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                                >  1:5474/1‑100 (MQ=255)
                                                   cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                               <  2:186338/100‑1 (MQ=255)
                                                   cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTAGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                               >  1:1971750/1‑100 (MQ=255)
                                                     tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCAATAACTTATGAACCCCGACGCCCAACAGGTTAAACAGTTTATGCTCa                                              >  2:2516491/1‑99 (MQ=255)
                                                     tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCaa                                             <  2:2301125/100‑1 (MQ=255)
                                                      gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCa                                              >  2:1524962/1‑98 (MQ=255)
                                                      gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                            >  2:981020/1‑100 (MQ=255)
                                                      gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAc                                            >  2:982427/1‑100 (MQ=255)
                                                       aTCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACACGACGAACACCAGGTAAAACAGTTTCTGCTCAAcc                                           >  1:856888/1‑100 (MQ=255)
                                                        tCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCt                                          >  2:683500/1‑100 (MQ=255)
                                                        tCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCt                                          >  1:1879483/1‑100 (MQ=255)
                                                         cAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCtt                                         >  2:2620655/1‑100 (MQ=255)
                                                          aGTTATTTCCACTGCTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTc                                        <  2:1809257/100‑1 (MQ=255)
                                                           gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                       <  1:1603569/100‑1 (MQ=255)
                                                            ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAg                                      <  1:2331553/100‑1 (MQ=255)
                                                            ttATTTCCACTGGTTCGAAAACTATAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAg                                      <  1:332269/100‑1 (MQ=255)
                                                             tATTTCCACTTGTTAGACAACAAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                     <  1:1897266/100‑1 (MQ=255)
                                                             tATTTCCACTGGTTCGATAACCAGAAAGTACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                     <  1:2363370/100‑1 (MQ=255)
                                                              aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                    <  1:949859/100‑1 (MQ=255)
                                                               tttCCACTGGTGCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTGAAACAGTTTCTGCTCAACCTTCa                                       <  2:1862040/96‑1 (MQ=255)
                                                                ttCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                                  >  2:2240132/1‑100 (MQ=255)
                                                                 tCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCGAAACCTTCAGGATa                                  >  1:2611118/1‑99 (MQ=255)
                                                                  ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                                >  2:1339206/1‑100 (MQ=255)
                                                                  ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                                >  2:1931561/1‑100 (MQ=255)
                                                                  ccACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                                >  2:1251955/1‑100 (MQ=255)
                                                                   cACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGa                               >  2:1599636/1‑100 (MQ=255)
                                                                    aCTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAt                              >  2:925021/1‑100 (MQ=255)
                                                                      tGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCGCACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                            <  2:1489987/100‑1 (MQ=255)
                                                                      tGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                            <  1:1191011/100‑1 (MQ=255)
                                                                       ggTTCGACAACCAGAAATCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAtttn                           <  1:1124928/100‑2 (MQ=255)
                                                                       ggTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                           >  1:1069379/1‑100 (MQ=255)
                                                                        gTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGt                          >  1:1637223/1‑100 (MQ=255)
                                                                         ttCGACAACCAGACAGCACATTCGAAAACTAGATAACTTATGACACCCGACGTACTCCAGTTTAAACTGTTTCTTCTTTTCTTTCAGGATACGATTTGTc                         <  2:588065/100‑1 (MQ=255)
                                                                         ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                         <  1:1920373/100‑1 (MQ=255)
                                                                          tCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTca                        <  2:1440689/100‑1 (MQ=255)
                                                                           cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                       <  1:1597555/100‑1 (MQ=255)
                                                                           cGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTACACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcag                       <  2:1412162/100‑1 (MQ=255)
                                                                            gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                      >  2:2424942/1‑100 (MQ=255)
                                                                            gACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                      >  2:1446949/1‑100 (MQ=255)
                                                                              cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                    >  2:563641/1‑100 (MQ=255)
                                                                              cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                    >  2:2414406/1‑100 (MQ=255)
                                                                              cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                    >  2:1218271/1‑100 (MQ=255)
                                                                               aaCCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCagcagc                   <  2:464412/100‑1 (MQ=255)
                                                                               aaCCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCagcagc                   <  1:142992/100‑1 (MQ=255)
                                                                                 ccAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTg                 >  1:2566752/1‑100 (MQ=255)
                                                                                 ccAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTg                 >  2:21821/1‑100 (MQ=255)
                                                                                   aGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc               <  2:806082/100‑1 (MQ=255)
                                                                                   aGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc               <  2:1098693/100‑1 (MQ=255)
                                                                                      aaGCACATTCGTAAAATCGATAAGTTATGTAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc               <  1:1145579/97‑1 (MQ=255)
                                                                                      aaGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgc            <  1:1259532/100‑1 (MQ=255)
                                                                                      aaGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgc            <  1:2000730/100‑1 (MQ=255)
                                                                                        gCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgccg          >  2:2453196/1‑100 (MQ=255)
                                                                                          acaTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTc        <  2:1637207/100‑1 (MQ=255)
                                                                                          acaTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTc        <  1:2377050/100‑1 (MQ=255)
                                                                                           caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg       >  2:1362641/1‑100 (MQ=255)
                                                                                                gAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGgc  >  2:572091/1‑100 (MQ=255)
                                                                                                 aaaaaGCGATAAGTTATGAAACCCTACGTTTTTCAGGTTTAATAGTTTCTGCTCAACCTTAAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATgg   <  1:719267/98‑1 (MQ=255)
                                                                                                 |                                                                                                  
GTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGC  >  NC_000913/2553113‑2553308

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