Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,403,542 A→G S342P (TCT→CCT)  abgA ← p‑aminobenzoyl‑glutamate hydrolase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,403,5420AG100.0% 1037.9 / NA 251S342P (TCT→CCT) abgAp‑aminobenzoyl‑glutamate hydrolase subunit A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (113/138);  total (113/138)

Only 100 of 267 total aligned reads displayed.

GGAACCCGCAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGAACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCAAAAACATAT  >  NC_000913/1403444‑1403639
                                                                                                  |                                                                                                 
ggAACCCGCAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGAc                                                                                                 >  1:2141214/1‑101 (MQ=255)
ggAACCCGCAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGAc                                                                                                 >  2:2141215/1‑101 (MQ=255)
  aaCCCGCAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTg                                                                                               >  2:2141218/1‑101 (MQ=255)
   aCCCGCAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTgg                                                                                              >  2:2141220/1‑101 (MQ=255)
        cAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTa                                                                                         >  1:2141228/1‑101 (MQ=255)
        cAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTa                                                                                         >  2:2141227/1‑101 (MQ=255)
         aGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAg                                                                                        >  1:2141232/1‑101 (MQ=255)
          ggCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGc                                                                                       >  2:2141233/1‑101 (MQ=255)
           gcgcTTCAACACGTTCAATGGCCTGATTGCCCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGCCCCATTGCGGCGAAGGAGGACTGGCGGTAGCt                                                                                      <  2:2140675/101‑1 (MQ=255)
           gcgcTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCt                                                                                      <  2:2140813/101‑1 (MQ=255)
            cgcTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTg                                                                                     <  2:2140692/101‑1 (MQ=255)
             gcTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTg                                                                                     <  1:2140562/100‑1 (MQ=255)
             gcTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGc                                                                                    >  2:2141238/1‑101 (MQ=255)
             gcTTCAACACGTTCAATGGCCTGATTAACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGc                                                                                    >  1:2141237/1‑101 (MQ=255)
              cTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCa                                                                                   <  2:2140718/101‑1 (MQ=255)
               ttCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCAc                                                                                  >  1:2141242/1‑101 (MQ=255)
                tCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCAcc                                                                                 >  1:2141244/1‑101 (MQ=255)
                  aaCACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTa                                                                                         <  2:2140743/91‑1 (MQ=255)
                  aaCACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCa                                                                               >  1:2141248/1‑101 (MQ=255)
                  aaCACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCa                                                                               <  1:2140725/101‑1 (MQ=255)
                   acacGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCAc                                                                                  >  1:2141251/1‑97 (MQ=255)
                   acacGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCAt                                                                              <  2:2140752/101‑1 (MQ=255)
                   acacGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCAt                                                                              <  1:2140804/101‑1 (MQ=255)
                   acacGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCAt                                                                              <  1:2140731/101‑1 (MQ=255)
                    cacGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATc                                                                             >  2:2141252/1‑101 (MQ=255)
                      cGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCAccc                                                                                <  2:2140779/96‑1 (MQ=255)
                       gTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGa                                                                          >  1:2141254/1‑101 (MQ=255)
                        ttCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGAc                                                                         >  1:2141258/1‑101 (MQ=255)
                        tgcAATGGCCTGATTGCCCCCCGCGCCCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCt                                                                                      <  1:2140767/86‑1 (MQ=255)
                         tCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACg                                                                        <  1:2140795/101‑1 (MQ=255)
                          cAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGa                                                                       <  2:2140980/101‑1 (MQ=255)
                           aaTGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAg                                                                      <  1:2140786/101‑1 (MQ=255)
                           aaTGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAg                                                                      <  2:2140952/101‑1 (MQ=255)
                            aTGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGt                                                                     >  2:2141261/1‑101 (MQ=255)
                              ggCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCa                                                                            >  2:2141267/1‑92 (MQ=255)
                              ggCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGttt                                                                   >  1:2141264/1‑101 (MQ=255)
                              ggCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGttt                                                                   >  1:2141266/1‑101 (MQ=255)
                                 cTGATTGCCCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCaa                                                                <  2:2140831/101‑1 (MQ=255)
                                  tGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAAc                                                               <  1:2140857/101‑1 (MQ=255)
                                   taTTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACg                                                              >  1:2141273/2‑101 (MQ=255)
                                   gATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACg                                                              >  2:2141272/1‑101 (MQ=255)
                                       gccccccGCGACCTGAGCCGCCTGACTTTGCAACCATGCGCCCCATTGCGGCGAAGGAGGACTgg                                                                                              <  2:2140866/63‑1 (MQ=255)
                                       gACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCa                                                                 <  1:2140771/94‑1 (MQ=255)
                                         cccccGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGaca                                                        >  2:2141278/1‑101 (MQ=255)
                                          ccccGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGacac                                                       >  1:2141279/1‑101 (MQ=255)
                                          ccccGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGacac                                                       >  2:2141281/1‑101 (MQ=255)
                                            ccGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCa                                                     >  2:2141285/1‑101 (MQ=255)
                                               cgACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATAc                                                  >  1:2141287/1‑101 (MQ=255)
                                                 aCCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACAt                                                <  1:2140874/101‑1 (MQ=255)
                                                  ccTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATg                                               <  1:2140926/101‑1 (MQ=255)
                                                   cTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGGCTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATgg                                              >  1:2141294/1‑101 (MQ=255)
                                                   cTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATgg                                              >  2:2141292/1‑101 (MQ=255)
                                                    tGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGt                                             <  1:2140962/101‑1 (MQ=255)
                                                     gAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGaca                                                        <  2:2141017/89‑1 (MQ=255)
                                                     gAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGtt                                            <  2:2140803/101‑1 (MQ=255)
                                                       gccgccTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTtgc                                          >  2:2141295/1‑101 (MQ=255)
                                                       gccgccTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTtgc                                          <  1:2140859/101‑1 (MQ=255)
                                                        ccgccTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTtgct                                         <  1:2140916/101‑1 (MQ=255)
                                                        ccgccTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTtgct                                         <  1:2140943/101‑1 (MQ=255)
                                                         cgccTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTtgctg                                        <  1:2140991/101‑1 (MQ=255)
                                                            cTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTgcgc                                     >  1:2141298/1‑101 (MQ=255)
                                                              gACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGccc                                   >  2:2141303/1‑101 (MQ=255)
                                                              gACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGccc                                   >  1:2141302/1‑101 (MQ=255)
                                                               aCTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCt                                  <  1:2140888/101‑1 (MQ=255)
                                                                cTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTg                                 <  2:2140955/101‑1 (MQ=255)
                                                                cTTTGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTg                                 <  2:2140932/101‑1 (MQ=255)
                                                                 tttGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGa                                <  2:2141055/101‑1 (MQ=255)
                                                                 tttGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGa                                >  1:2141306/1‑101 (MQ=255)
                                                                   tGCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAAt                              <  1:2140852/101‑1 (MQ=255)
                                                                    gCAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATc                             <  1:2140826/101‑1 (MQ=255)
                                                                     cAACCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCg                            <  2:2140982/101‑1 (MQ=255)
                                                                      aaCCATGCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGc                           >  1:2141311/1‑101 (MQ=255)
                                                                            gCGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCttgttg                     <  2:2140862/101‑1 (MQ=255)
                                                                             cGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGACTCGCTtgttgt                    <  1:2140818/101‑1 (MQ=255)
                                                                             cGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTtgttgt                    <  1:2140834/101‑1 (MQ=255)
                                                                             cGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTtgttgt                    <  1:2140997/101‑1 (MQ=255)
                                                                             cGACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTtgttgt                    >  2:2141313/1‑101 (MQ=255)
                                                                              gACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCttgt                       <  2:2140893/97‑1 (MQ=255)
                                                                              gACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTtgtg                   <  1:2140973/101‑1 (MQ=255)
                                                                              gACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTtgtg                   >  1:2141315/1‑101 (MQ=255)
                                                                              gACCCATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTtgtg                   >  2:2141317/1‑101 (MQ=255)
                                                                                 ccATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCa                                                                 <  1:2140759/52‑1 (MQ=255)
                                                                                 ccATTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCAc                <  2:2140914/101‑1 (MQ=255)
                                                                                   aTTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACgg              <  1:2140883/101‑1 (MQ=255)
                                                                                   aTTGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACgg              >  2:2141321/1‑101 (MQ=255)
                                                                                    ttGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGt             <  2:2140925/101‑1 (MQ=255)
                                                                                     tGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTc            >  1:2141323/1‑101 (MQ=255)
                                                                                     tGCGGCGAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTc            <  2:2141114/101‑1 (MQ=255)
                                                                                       cggcgAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCaa          <  1:2140902/101‑1 (MQ=255)
                                                                                        ggcgAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGttt                                                                   <  1:2140365/43‑1 (MQ=255)
                                                                                         gcgAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCaaaa        <  1:2141011/101‑1 (MQ=255)
                                                                                          cgAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCaaaaa       <  1:2140701/101‑1 (MQ=255)
                                                                                           gAAGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCAAAAAc      >  2:2141326/1‑101 (MQ=255)
                                                                                            aaGGAGGACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCAAAAACa     >  1:2141329/1‑101 (MQ=255)
                                                                                             aggaggACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCAAGAACat    >  1:2141331/1‑101 (MQ=255)
                                                                                             aggaggACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCAAAAACat    <  2:2140964/101‑1 (MQ=255)
                                                                                              ggaggaCTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATAc                                                  <  2:2140798/54‑1 (MQ=255)
                                                                                              ggaggaCTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTg                                 <  2:2140535/71‑1 (MQ=255)
                                                                                              ggaggaCTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCttgt                       <  2:2140910/81‑1 (MQ=255)
                                                                                              ggaggaCTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCAc                <  2:2140847/88‑1 (MQ=255)
                                                                                                  |                                                                                                 
GGAACCCGCAGGCGCTTCAACACGTTCAATGGCCTGATTGACCCCCGCGACCTGAGCCGCCTGACTTTGCAACCATGCGACCCATTGCGGCGAAGGAGAACTGGCGGTAGCTGCACCCATCAGACGAGTTTCAACGCCGACACCATACATGGTTGCTGCGCCCTGAATCGCTTGTTGTGCACGGTCAAAAACATAT  >  NC_000913/1403444‑1403639

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