Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,631,700 T→C 100% L323P (CTT→CCT)  slt → lytic murein transglycosylase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,7000TC97.4% 628.7 / 1.7 184L323P (CTT→CCT) sltlytic murein transglycosylase, soluble
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base C (75/104);  total (76/108)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 184 total aligned reads displayed.

GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATT  >  NC_000913/4631603‑4631795
                                                                                                 |                                                                                               
ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTg                                                                                               >  1:1158320/1‑100 (MQ=255)
ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTg                                                                                               <  2:143224/100‑1 (MQ=255)
ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCAATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTg                                                                                               <  1:126307/100‑1 (MQ=255)
 gCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTgg                                                                                              <  1:1298286/100‑1 (MQ=255)
 gCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTgg                                                                                              >  2:2150924/1‑100 (MQ=255)
  cAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGc                                                                                             >  2:1100966/1‑100 (MQ=255)
  cAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGc                                                                                             <  1:201538/100‑1 (MQ=255)
    aCGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAc                                                                                           >  1:2323548/1‑100 (MQ=255)
     cGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAcc                                                                                          <  2:1074890/100‑1 (MQ=255)
        tGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGc                                                                                       <  1:1143518/100‑1 (MQ=255)
        tGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGc                                                                                       <  2:3262145/100‑1 (MQ=255)
        tGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGc                                                                                       <  2:1176296/100‑1 (MQ=255)
         gTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCg                                                                                      <  1:1045905/100‑1 (MQ=255)
          tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGa                                                                                     >  2:502834/1‑100 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtc                                                                                   >  1:1097676/1‑100 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtc                                                                                   <  2:412886/100‑1 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGCACGAATGGCGCCTGGCACCGGCGAtc                                                                                   <  2:2057364/100‑1 (MQ=255)
             ccGACGAGCAGGCGAAATGGCGCGATGGCGCCCTTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGCATGGCGCCCGGCACCGGCGAtcg                                                                                  >  1:3290869/1‑100 (MQ=255)
             ccGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcg                                                                                  >  2:1396869/1‑100 (MQ=255)
             ccGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcg                                                                                  >  2:2138333/1‑100 (MQ=255)
             ccGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcg                                                                                  <  1:275526/100‑1 (MQ=255)
              cgacgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcgt                                                                                 >  1:2746349/1‑100 (MQ=255)
               gacgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcgtc                                                                                >  2:499644/1‑100 (MQ=255)
                   aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCgg                                                                            <  1:3383297/100‑1 (MQ=255)
                     cAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGcc                                                                          <  1:3203331/100‑1 (MQ=255)
                      aGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCt                                                                         <  2:1752176/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        >  2:3161171/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        >  1:2009646/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  2:2323484/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        >  2:1641022/1‑100 (MQ=255)
                       ggCGAAATGGAGCTATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  1:3102688/100‑1 (MQ=255)
                        gcataaTGGCGTCATCGCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACACCGTGTACGAATGGCGCCTGGCACCGACGATCGTCGCGGCCTGa                                                                       >  2:1123159/5‑100 (MQ=255)
                        gCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCTCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  1:1259310/99‑1 (MQ=255)
                        gCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGGCTGGCACCGGCGATCGTCGCGGCCt                                                                         >  1:3367772/1‑98 (MQ=255)
                          gAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAAt                                                                     <  1:2360553/100‑1 (MQ=255)
                          gAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAAt                                                                     <  1:212491/100‑1 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATAcccg                                                                >  1:918439/1‑98 (MQ=255)
                                 cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                              <  2:2442413/100‑1 (MQ=255)
                                 cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                              <  2:2671389/100‑1 (MQ=255)
                                 cgcgATGACGCCATTATGCGCACGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                              >  1:427655/1‑100 (MQ=255)
                                    gATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggctgg                                                           <  2:1724561/100‑1 (MQ=255)
                                     aTGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCtggctggc                                                          <  1:2450178/100‑1 (MQ=255)
                                      tGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcg                                                         <  2:2209703/100‑1 (MQ=255)
                                          gCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcg                                                         <  1:894388/96‑1 (MQ=255)
                                            ccgatctcgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                                   <  1:2763734/93‑1 (MQ=255)
                                             aTTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGc                                                  >  2:915808/1‑100 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                 <  2:1016969/100‑1 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                 <  2:2899695/100‑1 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                 <  2:236906/100‑1 (MQ=255)
                                               tATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCt                                                <  2:1801209/100‑1 (MQ=255)
                                                aTGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTa                                               <  2:250951/100‑1 (MQ=255)
                                                 tGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTAt                                              <  2:200332/100‑1 (MQ=255)
                                                  gcgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                             <  2:328386/100‑1 (MQ=255)
                                                  gcgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                             >  1:1746445/1‑100 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTgagatcgg                                            >  2:2763757/1‑93 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                            <  2:2977753/100‑1 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGGGTACGAAAGGCGCCTGGCCCCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGGCTGCCGATgg                                            >  1:1506131/1‑100 (MQ=255)
                                                   cgcTCGCAAGCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                            <  1:71525/100‑1 (MQ=255)
                                                    gcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGa                                           >  1:2404813/1‑100 (MQ=255)
                                                     cTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGaa                                          >  1:1835448/1‑100 (MQ=255)
                                                     cTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGaa                                          >  1:263988/1‑100 (MQ=255)
                                                       cGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGc                                        <  1:2198067/100‑1 (MQ=255)
                                                        gCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCg                                       <  2:836444/100‑1 (MQ=255)
                                                         cAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCga                                      <  2:1812728/100‑1 (MQ=255)
                                                         cAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCg                                       >  2:2776321/1‑99 (MQ=255)
                                                            tCTACGTCGCTTATTGAACGCCGTGGACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaag                                   <  1:2711877/100‑1 (MQ=255)
                                                             cTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                                  <  1:1641015/100‑1 (MQ=255)
                                                              tACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                                  <  1:3309486/99‑1 (MQ=255)
                                                                cTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgagagagaa                               >  2:1559581/1‑100 (MQ=255)
                                                                cTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaa                               >  2:306454/1‑100 (MQ=255)
                                                                cTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaa                               <  1:608090/100‑1 (MQ=255)
                                                                  tCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                             >  1:46397/1‑100 (MQ=255)
                                                                  tCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                             <  1:299926/100‑1 (MQ=255)
                                                                  tCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                             <  1:960132/100‑1 (MQ=255)
                                                                   cGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAATCgaaagagaaaga                            <  2:2276446/100‑1 (MQ=255)
                                                                   cGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaaga                            >  1:2500083/1‑100 (MQ=255)
                                                                   cGCTTATTCAACGCCGTGCACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaaga                            <  2:2081177/100‑1 (MQ=255)
                                                                    gCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGAt                           >  2:1496628/1‑100 (MQ=255)
                                                                    gCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGAt                           >  2:144959/1‑100 (MQ=255)
                                                                     ctaaTTGATAGCCGTGTACGAACGGAGACGGGGATCGGCGGTCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaag                                   <  1:3247389/88‑1 (MQ=255)
                                                                         ttGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATg                      >  2:38175/1‑100 (MQ=255)
                                                                         ttGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATg                      <  1:747140/100‑1 (MQ=255)
                                                                          tggaCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATgg                     <  1:2839973/97‑1 (MQ=255)
                                                                          tGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATgg                     <  2:2808999/100‑1 (MQ=255)
                                                                            aaCGCTGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCg                   >  1:2888594/1‑100 (MQ=255)
                                                                             aCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGt                  <  1:786941/100‑1 (MQ=255)
                                                                              cGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGtt                 <  1:2782269/100‑1 (MQ=255)
                                                                                 cGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACt              <  2:2427438/100‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc           <  2:1223745/100‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc           <  2:396366/100‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc           <  1:1286048/100‑1 (MQ=255)
                                                                                     tACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATAAATGGCGTTACTGGCa          >  2:702626/1‑100 (MQ=255)
                                                                                      aCGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAg         <  1:2420415/100‑1 (MQ=255)
                                                                                       cGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg        <  1:3062102/100‑1 (MQ=255)
                                                                                       cGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAAGGGCGTTACTGGCAgg        >  2:3333781/1‑100 (MQ=255)
                                                                                        gAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGc       >  2:2829630/1‑100 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg      >  2:2822912/1‑100 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg      >  2:628173/1‑100 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg      <  2:874638/100‑1 (MQ=255)
                                                                                             gcgcCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAtt  <  1:3381581/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATT  >  NC_000913/4631603‑4631795

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