Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,431,892 T→A 5.1% A36A (GCA→GCT)  PP_2130 ← lytic transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,431,8920TA 5.1% 234.6 / 3.6 78A36A (GCA→GCT) PP_2130lytic transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base A (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GAGTTCGTCGTAGGCCAGGTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTTAGCTTAGCTGTTTGCCGG  >  NC_002947/2431745‑2432036
                                                                                                                                                   |                                                                                                                                                
gAGTTCGTCGTAGGCCAGGTAGGGCGTCCGCGGGGAGACGGTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCctgcgctgc                                                                                                                                                <  2:776619/150‑1 (MQ=255)
  gTTCGTCGTAGGCCAGGTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCtgcgctgcgc                                                                                                                                              >  2:989234/1‑150 (MQ=255)
   ttCGTCGTAGGCCAGGTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCa                                                                                                                                             >  1:1296016/1‑150 (MQ=255)
   ttCGTCGTAGGCCAGGTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCa                                                                                                                                             >  2:1329461/1‑150 (MQ=255)
            ggCCAGGTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTGATAGTACTGGCGTTGCAGGGTCAGGTCGGGGGGCTGCGCTGCGCGGGTGGCAgg                                                                                                                                    >  1:2088810/1‑149 (MQ=255)
              ccAGGTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCgg                                                                                                                                  <  1:890871/150‑1 (MQ=255)
               cAGGTAGGGCGTCAGCGGGTCGTCGCTCAAGGTATGGGCATATCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGt                                                                                                                                 <  2:21722/150‑1 (MQ=255)
                  gTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAag                                                                                                                              >  2:1275038/1‑150 (MQ=255)
                    aGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAagca                                                                                                                            <  2:1913894/150‑1 (MQ=255)
                    aGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGACACCGGTAatca                                                                                                                            <  2:2330278/150‑1 (MQ=255)
                      ggCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAagcagc                                                                                                                          <  2:1784629/150‑1 (MQ=255)
                       gCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCAGCGCAggccacc                                                                                                                                      >  2:565915/1‑132 (MQ=255)
                       gCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCAGCGCAggccacc                                                                                                                                      <  1:565915/137‑6 (MQ=255)
                         gTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAgcagg                                                                                                                       <  1:2184729/150‑1 (MQ=255)
                         gTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGCTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCCgcagg                                                                                                                       <  1:1461284/150‑1 (MQ=255)
                          tCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAgcaggc                                                                                                                      >  1:1163038/1‑150 (MQ=255)
                           cAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAgcaggca                                                                                                                     <  2:1043615/150‑1 (MQ=255)
                           cAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAgcaggca                                                                                                                     <  1:441989/150‑1 (MQ=255)
                           cAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCATTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAgcaggca                                                                                                                     <  1:841482/150‑1 (MQ=255)
                             gCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGcaggcagg                                                                                                                   >  2:1941128/1‑150 (MQ=255)
                              cGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGa                                                                                                                  >  1:1575263/1‑150 (MQ=255)
                                   aGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGGGGAAGCGGTAAGCAGCAGGCAGGATGCGa                                                                                                             <  1:1275038/150‑1 (MQ=255)
                                      cGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGatat                                                                                                          <  1:311070/150‑1 (MQ=255)
                                       gCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATg                                                                                                         >  2:1052439/1‑150 (MQ=255)
                                       gCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATg                                                                                                         >  2:923881/1‑150 (MQ=255)
                                         gcAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTa                                                                                                       >  2:1926204/2‑150 (MQ=255)
                                          cAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTaa                                                                                                      <  2:2266712/150‑1 (MQ=255)
                                          cAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTaa                                                                                                      <  2:939431/150‑1 (MQ=255)
                                           aaGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAAc                                                                                                     <  2:1163038/150‑1 (MQ=255)
                                               ccTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAggc                                                                                                 <  2:1359268/150‑1 (MQ=255)
                                                cTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAggcg                                                                                                >  2:2168926/1‑150 (MQ=255)
                                                 tGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAggcgg                                                                                               >  1:1884663/1‑150 (MQ=255)
                                                         gcgcAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCAt                                                                                       <  1:923881/150‑1 (MQ=255)
                                                           gcAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCg                                                                                                              >  1:1322181/1‑125 (MQ=255)
                                                           gcAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCg                                                                                                              <  2:1322181/125‑1 (MQ=255)
                                                                   gggCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCg                                                                             >  2:388293/1‑150 (MQ=255)
                                                                        tttGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAg                                                                        <  2:1802074/150‑1 (MQ=255)
                                                                         ttGTCGCCTTTGGTCAGCCCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGTGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGt                                                                       <  2:987916/150‑1 (MQ=255)
                                                                           gTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTg                                                                     >  1:216769/1‑150 (MQ=255)
                                                                           gTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTg                                                                     >  2:1001007/1‑150 (MQ=255)
                                                                            tCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTAGTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTcc                                                                              >  2:1350090/1‑140 (MQ=255)
                                                                            tCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTcc                                                                              <  1:1350090/140‑1 (MQ=255)
                                                                            tCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCAGCGCAggcca                                                                                                                                        <  1:1764448/82‑4 (MQ=255)
                                                                            tCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCAGCGCAggcca                                                                                                                                        >  2:1764448/1‑79 (MQ=255)
                                                                                 tttGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCAt                                                               <  2:1296016/150‑1 (MQ=255)
                                                                                     gCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGa                                                           >  1:752325/1‑150 (MQ=255)
                                                                                             tcgctTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGTAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAggcgg                                                                                               <  2:2216252/106‑1 (MQ=255)
                                                                                             tcgctTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGTAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAggcgg                                                                                               >  1:2216252/1‑106 (MQ=255)
                                                                                             tcgctTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGCGTcgg                                                   <  1:2168926/150‑1 (MQ=255)
                                                                                                 ttGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGa                                               <  1:388293/150‑1 (MQ=255)
                                                                                                   ggCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGTAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGc                                             <  2:1507030/150‑1 (MQ=255)
                                                                                                         gTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAgcagca                                                                                                                         <  1:983250/68‑1 (MQ=255)
                                                                                                         gTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAgcagca                                                                                                                         >  2:983250/1‑68 (MQ=255)
                                                                                                            aTAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTAt                                    >  2:577760/1‑150 (MQ=255)
                                                                                                             tAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGCTCTGTAACAGGCGGCTGCGCATGATAGGTCCAGGTATTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTAtt                                   <  1:989234/150‑1 (MQ=255)
                                                                                                             tAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTAtt                                   <  2:216769/150‑1 (MQ=255)
                                                                                                                tACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGag                                <  2:1745739/150‑2 (MQ=255)
                                                                                                                 aCTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGaaa                               >  1:484164/1‑150 (MQ=255)
                                                                                                                 aCTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGaaa                               >  1:125082/1‑150 (MQ=255)
                                                                                                                 aCTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGaaa                               >  2:854360/1‑150 (MQ=255)
                                                                                                                 aCTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGaaa                               <  2:2088810/150‑1 (MQ=255)
                                                                                                                   tGGCGTTGCAGGGTCAGGTCGGTGGCCCGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGc                             >  1:2325696/1‑150 (MQ=255)
                                                                                                                          gCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGcc                      <  2:125082/150‑1 (MQ=255)
                                                                                                                            aGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTcggcgg                                                >  2:2186030/1‑122 (MQ=255)
                                                                                                                            aGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTcggcgg                                                <  1:2186030/122‑1 (MQ=255)
                                                                                                                            aGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCctt                    >  1:1555819/1‑150 (MQ=255)
                                                                                                                                tCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCcttagct                >  1:2304751/1‑150 (MQ=255)
                                                                                                                                tCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCctta                   <  2:1547638/147‑1 (MQ=255)
                                                                                                                                tCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCctta                   >  1:1547638/1‑147 (MQ=255)
                                                                                                                                  aGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCcttagctta              <  1:1926204/150‑1 (MQ=255)
                                                                                                                                    gTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCttagcttagc            >  2:2220024/1‑150 (MQ=255)
                                                                                                                                    gTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTAGTGAAAGCTATAGActtagcttagg            >  2:1302962/1‑149 (MQ=255)
                                                                                                                                     tCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTtagcttagct           <  1:577760/150‑1 (MQ=255)
                                                                                                                                     tCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTtagcttagct           <  1:1052439/150‑1 (MQ=255)
                                                                                                                                        gTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATAc                                                                                  <  2:540124/76‑1 (MQ=255)
                                                                                                                                        gTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATAc                                                                                  >  1:540124/1‑76 (MQ=255)
                                                                                                                                        gTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCcttagcttagggatt        >  2:1668893/1‑145 (MQ=255)
                                                                                                                                         tGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGc                                             <  1:1375173/112‑1 (MQ=255)
                                                                                                                                         tGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGc                                             >  2:1375173/1‑112 (MQ=255)
                                                                                                                                           gCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTTAGCTTAGCTGTTTGc     >  1:936261/1‑150 (MQ=255)
                                                                                                                                             ctgcgctgcgCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTTAGCTTAGCTGTTTGCCg   >  2:2190930/1‑150 (MQ=255)
                                                                                                                                             ctgcgctgcgCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTTAGCTTAGCTGTTTGCCg   >  2:1105982/1‑150 (MQ=255)
                                                                                                                                              tgcgctgcgcAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTTAGCTTAGCTGTTTGCCgg  >  1:472905/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GAGTTCGTCGTAGGCCAGGTAGGGCGTCAGCGGGTAGTCGCTCAAGGCCTGGGCATAGCGCAGGTACGGGCCTTTGTCGCCTTTGGCCAGCGCTCGCTTGGCCTCGTCATAGTACTGGCGTTGCAGGGTCAGGTCGGTGGCCTGCGCTGCGCAGGTGGCAGCGGTAAGCAGCAGGCAGGATGCGATATGTAACAGGCGGCTGCGCATGATACGTCCGGGCAGTTGAACCATGGGAAGTGTCGGCGGAGCCAGCACTATTGAAAGCTATTGCCTTAGCTTAGCTGTTTGCCGG  >  NC_002947/2431745‑2432036

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