Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,525,988 C→G 5.3% T85R (ACG→AGG)  murD → UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,525,9880CG 5.3% 1217.1 / 40.9 357T85R (ACG→AGG) murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base C (161/177);  new base G (11/8);  total (172/185)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 359 total aligned reads displayed.

TGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCC  >  NC_002947/1525839‑1526136
                                                                                                                                                     |                                                                                                                                                    
    gTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCGTGGCCAAGGCCAGGCCg                                                                                                                                                  <  1:1712838/150‑1 (MQ=255)
     tCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCGTGGCCAAGGCCAGGCCgg                                                                                                                                                 <  1:3308507/150‑1 (MQ=255)
     tCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCGTGGCCAAGGCCAGGCCgg                                                                                                                                                 <  1:4063106/150‑1 (MQ=255)
     tCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCGTGGCCAAGGCCAGGCCgg                                                                                                                                                 <  1:4282869/150‑1 (MQ=255)
     tCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCgg                                                                                                                                                 >  1:5394971/1‑150 (MQ=255)
      cGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGc                                                                                                                                                <  1:5983259/150‑1 (MQ=255)
      cGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGc                                                                                                                                                <  2:3149568/150‑1 (MQ=255)
        ccGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCt                                                                                                                                              >  2:6128964/1‑150 (MQ=255)
          gATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTg                                                                                                                                            >  1:362565/1‑150 (MQ=255)
            tACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCACCGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgca                                                                                                                                          <  2:3832224/150‑1 (MQ=255)
            tACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgca                                                                                                                                          <  2:5011592/150‑1 (MQ=255)
            tACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGACCACGCCGGCCTTgca                                                                                                                                          <  1:5190346/150‑1 (MQ=255)
               ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                       >  2:5512655/1‑150 (MQ=255)
                 gcgAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAgg                                                                                                                                     <  2:5881238/150‑1 (MQ=255)
                 gcgAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAgg                                                                                                                                     >  2:3385044/1‑150 (MQ=255)
                      cagcCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGcc                                                                                                                                <  1:5138053/150‑1 (MQ=255)
                        gccgccGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCcgc                                                                                                                              >  2:3535938/1‑150 (MQ=255)
                         ccgccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCcgcg                                                                                                                             >  1:111123/1‑150 (MQ=255)
                          cgccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCcgcgc                                                                                                                            <  1:1521285/150‑1 (MQ=255)
                            ccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGGGAACTGGACGTGGACGTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgt                                                                                                                          <  1:38732/150‑1 (MQ=255)
                             cgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtg                                                                                                                         <  1:2046237/150‑1 (MQ=255)
                                 aCTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtggcgt                                                                                                                     <  1:1880836/150‑1 (MQ=255)
                                   tGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGa                                                                                                                   >  2:5751854/1‑150 (MQ=255)
                                    ggAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGaa                                                                                                                  <  2:30472/150‑1 (MQ=255)
                                    ggAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGaa                                                                                                                  <  2:1312758/150‑1 (MQ=255)
                                     gAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAg                                                                                                                 <  1:423958/150‑1 (MQ=255)
                                          ccTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTc                                                                                                            >  1:1441359/1‑150 (MQ=255)
                                           cTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCt                                                                                                           >  1:465233/1‑150 (MQ=255)
                                            tGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTg                                                                                                          >  2:4017348/1‑150 (MQ=255)
                                                ccGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                      >  1:4182462/1‑150 (MQ=255)
                                                 cGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGAc                                                                                                     >  1:2934975/1‑150 (MQ=255)
                                                 cGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGAc                                                                                                     >  2:2844204/1‑150 (MQ=255)
                                                      ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                                <  2:2520097/150‑1 (MQ=255)
                                                      ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                                <  1:1036179/150‑1 (MQ=255)
                                                          ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGTCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTg                                                                                            >  1:3176533/1‑150 (MQ=255)
                                                          ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGAGTGGCGTGAAGCTGTCTGGCGACATAGAGCTg                                                                                            >  1:5058863/1‑150 (MQ=255)
                                                           cGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                           <  2:6560277/150‑1 (MQ=255)
                                                           cGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                           >  2:5989401/1‑150 (MQ=255)
                                                                gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCt                                                                                      <  2:5023732/150‑1 (MQ=255)
                                                                gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCt                                                                                      <  2:3672098/150‑1 (MQ=255)
                                                                gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCt                                                                                      <  2:4461482/150‑1 (MQ=255)
                                                                 tGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                                     >  2:3616565/1‑150 (MQ=255)
                                                                  ggAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCg                                                                                    >  2:2904625/1‑150 (MQ=255)
                                                                      gTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCAt                                                                                >  2:5587807/1‑150 (MQ=255)
                                                                        gcgcTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGc                                                                              <  2:6584932/150‑1 (MQ=255)
                                                                           cTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCaa                                                                           <  2:3176338/150‑1 (MQ=255)
                                                                                gtgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGgcgc                                                                      <  1:2259641/150‑1 (MQ=255)
                                                                                gtgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGgcgc                                                                      <  2:1070105/150‑1 (MQ=255)
                                                                                gtgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGgcgc                                                                      >  1:860493/1‑150 (MQ=255)
                                                                                 tgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGcc                                                                     >  1:4483004/1‑150 (MQ=255)
                                                                                     cTGGACGTGGACGTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATc                                                                 <  1:4873112/150‑1 (MQ=255)
                                                                                        gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                              <  2:1167079/150‑1 (MQ=255)
                                                                                        gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                              <  2:3478777/150‑1 (MQ=255)
                                                                                        gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                              <  1:1379330/150‑1 (MQ=255)
                                                                                        gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                              <  2:2120317/150‑1 (MQ=255)
                                                                                        gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                              <  1:4305571/150‑1 (MQ=255)
                                                                                        gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCAtt                                                              >  2:3332759/1‑150 (MQ=255)
                                                                                         aCGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTg                                                             >  2:629184/1‑150 (MQ=255)
                                                                                           gTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGcc                                                           <  2:4286598/150‑1 (MQ=255)
                                                                                           gTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGcc                                                           <  1:981052/150‑1 (MQ=255)
                                                                                           gTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGcc                                                           >  1:44115/1‑150 (MQ=255)
                                                                                              gACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATc                                                        >  2:4588997/1‑150 (MQ=255)
                                                                                              gACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATc                                                        >  2:4083258/1‑150 (MQ=255)
                                                                                              gACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATc                                                        >  1:1928121/1‑150 (MQ=255)
                                                                                                gttCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAg                                                      >  1:6817861/2‑150 (MQ=255)
                                                                                                cTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAg                                                      >  1:1198450/1‑150 (MQ=255)
                                                                                                 ttCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGt                                                     >  2:520795/1‑150 (MQ=255)
                                                                                                   ccTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTgg                                                   <  2:2323801/150‑1 (MQ=255)
                                                                                                   ccTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTgg                                                   <  2:1799643/150‑1 (MQ=255)
                                                                                                    cTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGt                                                  <  1:2602943/150‑1 (MQ=255)
                                                                                                    cTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGt                                                  >  1:3135843/1‑150 (MQ=255)
                                                                                                     tgtgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGtt                                                 >  2:5397681/1‑150 (MQ=255)
                                                                                                     tgtgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTCGCCACGCCGGCCTTGCAGCAGGCTGCcgcgc                                                                                                                            <  1:4901744/75‑1 (MQ=255)
                                                                                                     tgtgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTCGCCACGCCGGCCTTGCAGCAGGCTGCcgcgc                                                                                                                            >  2:4901744/1‑75 (MQ=255)
                                                                                                       tgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTcc                                               >  1:2791313/1‑150 (MQ=255)
                                                                                                          cgcgCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGTCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAAc                                            <  2:716936/150‑1 (MQ=255)
                                                                                                               cAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCaa                                       >  2:147294/1‑150 (MQ=255)
                                                                                                               cAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCaa                                       <  2:806949/150‑1 (MQ=255)
                                                                                                                aaCGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGcc                                                                     <  1:649624/119‑1 (MQ=255)
                                                                                                                      cTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCAcc                                <  2:3791841/150‑1 (MQ=255)
                                                                                                                      cTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCAcc                                <  2:1607402/150‑1 (MQ=255)
                                                                                                                        gTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCTGTGGTTCCAACGCCAAGAGCACCGt                              >  2:4190473/1‑150 (MQ=255)
                                                                                                                         tACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTc                             <  1:4685052/150‑1 (MQ=255)
                                                                                                                            gTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacg                          >  2:5298131/1‑150 (MQ=255)
                                                                                                                              gAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacgac                        <  1:3342673/150‑1 (MQ=255)
                                                                                                                                 ggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCt                     >  2:4875522/3‑150 (MQ=255)
                                                                                                                                 ccccGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCt                     >  2:652424/1‑150 (MQ=255)
                                                                                                                                  cccGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTg                    <  1:5512551/150‑1 (MQ=255)
                                                                                                                                      gcctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCg                <  2:1852886/150‑1 (MQ=255)
                                                                                                                                         tggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTg             >  1:6318987/1‑150 (MQ=255)
                                                                                                                                         tggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTg             <  1:2198221/150‑1 (MQ=255)
                                                                                                                                         tggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTg             >  2:1909707/1‑150 (MQ=255)
                                                                                                                                         tggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTg             <  2:3342233/150‑1 (MQ=255)
                                                                                                                                           gcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGaa           >  1:5735939/1‑150 (MQ=255)
                                                                                                                                           gcctTGGCCACGCCGGACTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGAGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGaa           >  2:5175959/1‑150 (MQ=255)
                                                                                                                                            cctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGaaa          <  2:6346038/150‑1 (MQ=255)
                                                                                                                                             ctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAAt         >  1:3657225/1‑150 (MQ=255)
                                                                                                                                                 gCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGgcc     <  1:6128964/150‑1 (MQ=255)
                                                                                                                                                 gCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGgcc     <  1:6837767/150‑1 (MQ=255)
                                                                                                                                                 gCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGgcc     >  2:907236/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCC  >  NC_002947/1525839‑1526136

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