Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,700 T→C 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,7000TC98.6% 490.0 / ‑3.6 142L323P (CTT→CCT) sltlytic murein transglycosylase, soluble
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (69/71);  total (71/71)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 142 total aligned reads displayed.

GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATT  >  NC_000913/4631603‑4631795
                                                                                                 |                                                                                               
ggCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTg                                                                                               >  2:3381914/1‑100 (MQ=255)
 gCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTgg                                                                                              >  1:612151/1‑100 (MQ=255)
  cAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTgg                                                                                              >  1:1090527/1‑99 (MQ=255)
  cAACGATGTCACCGACGAGCAGGCGAAATAGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGc                                                                                             >  2:2400200/1‑100 (MQ=255)
   aaCGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGc                                                                                             <  1:3332461/99‑1 (MQ=255)
     cGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAcc                                                                                          <  2:790520/100‑1 (MQ=255)
     cGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAcc                                                                                          <  1:1065400/100‑1 (MQ=255)
      gATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAACCTACTTCGCTTATTCAACGCCGTGTACGAATGGCGCCTGCCAcc                                                                                          >  1:1149120/1‑99 (MQ=255)
          tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGa                                                                                     >  1:1137086/1‑100 (MQ=255)
          tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGa                                                                                     >  1:1698169/1‑100 (MQ=255)
          tCACCAACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGa                                                                                     >  2:118749/1‑100 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtc                                                                                   <  1:1921796/100‑1 (MQ=255)
              cgacgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcgt                                                                                 <  2:776221/100‑1 (MQ=255)
                acgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcgtcg                                                                               <  2:1577252/100‑1 (MQ=255)
                 cgaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTcgc                                                                              <  2:3028017/100‑1 (MQ=255)
                  gaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTcgcg                                                                             <  2:2187495/100‑1 (MQ=255)
                  gaGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTcgcg                                                                             <  2:1678521/100‑1 (MQ=255)
                   aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCgg                                                                            >  2:949214/1‑100 (MQ=255)
                   aGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCgg                                                                            >  2:444653/1‑100 (MQ=255)
                    gCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGc                                                                           >  2:2977569/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  1:2655489/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  1:1703987/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  1:1584863/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  1:2330608/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        >  2:1428364/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTGCGCTTATTGAGCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        <  1:708253/100‑1 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCACAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                        >  2:61561/1‑100 (MQ=255)
                            aaTGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATAc                                                                   <  2:2163842/100‑1 (MQ=255)
                            aaTGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATAc                                                                   >  1:2517950/1‑100 (MQ=255)
                            aaTGGCGCGATGACGCCATTATGCGCGCGCGATCTACTTCGCTTAGTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATAc                                                                   <  2:3167930/100‑1 (MQ=255)
                               ggcgtcatcgcGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                >  2:1500832/11‑100 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                >  2:2014123/1‑100 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                >  1:1918347/1‑100 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                <  1:2793708/100‑1 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                >  1:3444418/1‑100 (MQ=255)
                               ggCGCGATGACGCCATTATGCGCTCGCAATCTACGTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctg                                                                <  2:1167508/100‑1 (MQ=255)
                                gcgcGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctgg                                                               <  2:1250330/100‑1 (MQ=255)
                                 cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                              >  1:874517/1‑100 (MQ=255)
                                 cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                              >  2:3040391/1‑100 (MQ=255)
                                 cgcgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggc                                                              >  1:1061684/1‑100 (MQ=255)
                                   cgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggctg                                                            >  2:1304667/1‑100 (MQ=255)
                                   cgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggctg                                                            <  1:259970/100‑1 (MQ=255)
                                   cgATGACGCCATTATGCGCTCGCAATCTACGTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggctg                                                            <  1:259741/100‑1 (MQ=255)
                                     atgatgCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCtggctggc                                                          <  2:3410447/100‑1 (MQ=255)
                                     aTGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCtggctggc                                                          <  2:686285/100‑1 (MQ=255)
                                      tGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcg                                                         <  2:1656399/100‑1 (MQ=255)
                                      tGACGCCATTATGCGCTCGCAATCTACTACGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCATGAATACCTCGCTGgcg                                                         >  1:1182665/1‑100 (MQ=255)
                                        aCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                       >  1:2637478/1‑100 (MQ=255)
                                         cGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGgcgc                                                        >  1:557668/1‑98 (MQ=255)
                                         cGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGt                                                      >  2:2282969/1‑100 (MQ=255)
                                            cATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTg                                                   <  1:542105/100‑1 (MQ=255)
                                             aTTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGc                                                  >  2:2100131/1‑100 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                 >  1:1954472/1‑100 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                 <  2:758122/100‑1 (MQ=255)
                                               tATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCt                                                >  2:2637452/1‑100 (MQ=255)
                                                  gcgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                             <  1:2100100/100‑1 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                            >  2:1909301/1‑100 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                            >  2:2086565/1‑100 (MQ=255)
                                                      tCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAg                                         >  2:3225530/1‑100 (MQ=255)
                                                       cGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGc                                        <  2:568836/100‑1 (MQ=255)
                                                       cGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGc                                        <  1:3071448/100‑1 (MQ=255)
                                                        gCAATGTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAGCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCg                                       <  1:1455673/100‑1 (MQ=255)
                                                          aaTCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaa                                     >  1:1558090/1‑100 (MQ=255)
                                                            tCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaag                                   <  1:2315081/100‑1 (MQ=255)
                                                             cTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                                  <  2:344068/100‑1 (MQ=255)
                                                              tACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagag                                 <  2:2572300/100‑1 (MQ=255)
                                                               aCTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagaga                                >  2:2224353/1‑100 (MQ=255)
                                                                  tCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                             <  1:3265684/100‑1 (MQ=255)
                                                                   cGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaaga                            >  1:2727438/1‑100 (MQ=255)
                                                                     cTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                          >  1:421502/1‑100 (MQ=255)
                                                                     cTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                          >  2:888040/1‑100 (MQ=255)
                                                                     cTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATg                          <  1:73850/100‑1 (MQ=255)
                                                                      ttATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGa                         <  1:619046/100‑1 (MQ=255)
                                                                       tATTGAACGCCGTGTACGAATGGCGCCTGGCCCCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTCTGCGAGCGAACGAGAACGCTGac                        >  1:3130400/1‑99 (MQ=255)
                                                                       tATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGaa                        <  1:832075/100‑1 (MQ=255)
                                                                          tGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATgg                     >  1:2828746/1‑100 (MQ=255)
                                                                           gAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                    >  2:44942/1‑100 (MQ=255)
                                                                           gAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGc                    >  1:1098952/1‑100 (MQ=255)
                                                                             aCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGt                  >  2:750511/1‑100 (MQ=255)
                                                                             aCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGt                  >  1:1727804/1‑100 (MQ=255)
                                                                              cGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGtt                 <  1:1768927/100‑1 (MQ=255)
                                                                                ccGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTAc               <  1:1410105/100‑1 (MQ=255)
                                                                                  gtgtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTg             >  1:2791223/1‑100 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTgg            >  1:2358020/1‑99 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc           <  1:2476079/100‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc           <  1:767993/100‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc           <  2:2852473/100‑1 (MQ=255)
                                                                                     tACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCa          >  2:1476252/1‑100 (MQ=255)
                                                                                      aCGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAg         <  1:175826/100‑1 (MQ=255)
                                                                                      aCGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAg         <  1:1022715/100‑1 (MQ=255)
                                                                                       cGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg        >  2:1601325/1‑100 (MQ=255)
                                                                                       cGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg        >  2:2422160/1‑100 (MQ=255)
                                                                                         atgtgcgcCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg      >  1:2075932/5‑100 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg      >  2:1737425/1‑100 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg      <  2:2791267/100‑1 (MQ=255)
                                                                                         aaTGGCGCCTGGCACCGGCGATCGTCGCGACCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCg      >  2:1565222/1‑100 (MQ=255)
                                                                                            ggCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAt   >  1:3436603/1‑100 (MQ=255)
                                                                                            ggCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAt   >  1:1845144/1‑100 (MQ=255)
                                                                                            gacgcCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAt   >  2:2352963/3‑100 (MQ=255)
                                                                                             gcgcCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAtt  >  1:2085991/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATT  >  NC_000913/4631603‑4631795

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