Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,449,020 C→A 6.5% A148A (GCC→GCA) ‡ nagZ → beta‑N‑acetylglucosaminidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,449,0200CA 6.5% 605.3 / 22.4 183A148A (GCC→GCA) ‡nagZbeta‑N‑acetylglucosaminidase
Reads supporting (aligned to +/- strand):  ref base C (86/85);  new base A (6/6);  total (92/91)
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.93e-01

Only 100 of 185 total aligned reads displayed.

GAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGC  >  NC_002947/2448871‑2449168
                                                                                                                                                     |                                                                                                                                                    
gAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGcc                                                                                                                                                      <  1:2544022/150‑1 (MQ=255)
  gCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCgg                                                                                                                                                    <  1:1218753/150‑1 (MQ=255)
     gTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCGGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGgcgc                                                                                                                                                 <  1:2678988/150‑1 (MQ=255)
      tGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGGCCCGCTGCTTGCAACGCGCCTGGCCGGcgcg                                                                                                                                                >  1:1251284/1‑150 (MQ=255)
       gCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGGCCCGCTGCTTGCAACGCGCCTGGCCGGCGCGt                                                                                                                                               <  2:1251283/150‑1 (MQ=255)
       gCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGt                                                                                                                                               <  2:2253590/150‑1 (MQ=255)
        cGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGtt                                                                                                                                              <  2:642351/150‑1 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                      >  2:2230334/1‑150 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                      >  1:211223/1‑150 (MQ=255)
                   tGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCa                                                                                                                                   <  2:5414420/150‑1 (MQ=255)
                       gACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGaa                                                                                                                               >  1:1784404/1‑150 (MQ=255)
                            aGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATgcgg                                                                                                                          >  1:4358454/1‑150 (MQ=255)
                             ggTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATgcggc                                                                                                                         <  1:4820816/150‑1 (MQ=255)
                              gTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATgcgtcg                                                                                                                        <  2:2843268/150‑1 (MQ=255)
                               tGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGcggcgg                                                                                                                       <  2:5030355/150‑1 (MQ=255)
                                gCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCggg                                                                                                                      >  1:1609971/1‑150 (MQ=255)
                                gCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGCCGCATTCATgccacgga                                                                                                                      <  2:1212667/150‑7 (MQ=255)
                                gCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGCCGCATTCATgccacgga                                                                                                                      <  2:5113677/150‑7 (MQ=255)
                                  tGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                      <  1:235943/132‑1 (MQ=255)
                                  tGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                      >  2:235943/1‑132 (MQ=255)
                                  tGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGtt                                                                                                                                              <  1:635164/124‑1 (MQ=255)
                                  tGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGtt                                                                                                                                              >  2:634726/1‑124 (MQ=255)
                                    gCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATg                                                                                                                  >  2:5449433/1‑150 (MQ=255)
                                      ggTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGgc                                                                                                                >  2:1794718/1‑150 (MQ=255)
                                           gCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCt                                                                                                           <  2:319511/150‑1 (MQ=255)
                                            ccTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCtg                                                                                                          >  2:1114197/1‑150 (MQ=255)
                                               ggACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTgg                                                                                                       <  1:1234674/150‑1 (MQ=255)
                                                  ccTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCaa                                                                                                    <  1:3921999/150‑1 (MQ=255)
                                                   cTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAg                                                                                                   >  2:2009554/1‑150 (MQ=255)
                                                     cAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCa                                                                                                 <  2:3553542/150‑1 (MQ=255)
                                                     cAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCa                                                                                                 >  1:47256/1‑150 (MQ=255)
                                                         ttCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTc                                                                                             <  2:4849540/150‑1 (MQ=255)
                                                         ttCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTc                                                                                             <  2:3013295/150‑1 (MQ=255)
                                                          tCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTcc                                                                                            <  1:3628311/150‑1 (MQ=255)
                                                          tCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTcc                                                                                            <  2:78658/150‑1 (MQ=255)
                                                          tCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTcc                                                                                            <  2:3858853/150‑1 (MQ=255)
                                                            gCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCg                                                                                          <  2:2522823/150‑1 (MQ=255)
                                                            gCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCg                                                                                          <  2:3892601/150‑1 (MQ=255)
                                                             cTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCgg                                                                                         <  2:2642477/150‑1 (MQ=255)
                                                             cTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCgg                                                                                         >  2:4140262/1‑150 (MQ=255)
                                                              tCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCggg                                                                                        >  2:4980018/1‑150 (MQ=255)
                                                                  gTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCAt                                                                                    <  1:2230334/150‑1 (MQ=255)
                                                                     cTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGc                                                                                 >  1:2026632/1‑150 (MQ=255)
                                                                      tCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCt                                                                                >  1:290384/1‑150 (MQ=255)
                                                                       cGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCgggg                                                                                       <  2:4977671/142‑1 (MQ=255)
                                                                        gACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTgg                                                                              >  1:725754/1‑150 (MQ=255)
                                                                          ccTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGc                                                                            >  1:3325602/1‑150 (MQ=255)
                                                                          ccTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGc                                                                            <  2:3609756/150‑1 (MQ=255)
                                                                          ccTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGc                                                                            <  2:4592026/150‑1 (MQ=255)
                                                                           cTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGcc                                                                           >  1:2805262/1‑150 (MQ=255)
                                                                              gACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGaa                                                                        >  2:3069978/1‑150 (MQ=255)
                                                                              gACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGaa                                                                        <  1:412364/150‑1 (MQ=255)
                                                                                ccaccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGc                                                                      <  2:3439412/150‑1 (MQ=255)
                                                                                 caccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGcc                                                                     >  2:3825721/1‑150 (MQ=255)
                                                                                 caccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGcc                                                                     >  1:1578990/1‑150 (MQ=255)
                                                                                  accaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCg                                                                    <  2:2634930/150‑1 (MQ=255)
                                                                                  accaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCg                                                                    <  1:5086812/150‑1 (MQ=255)
                                                                                   ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGcggcgg                                                                                                                       >  2:4676137/1‑98 (MQ=255)
                                                                                   ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGcggcgg                                                                                                                       <  1:4676137/98‑1 (MQ=255)
                                                                                   ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTg                                                                               >  2:2097198/1‑138 (MQ=255)
                                                                                   ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTg                                                                               <  1:2097198/138‑1 (MQ=255)
                                                                                     aGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGAtt                                                                 <  1:2598347/150‑1 (MQ=255)
                                                                                      gCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTc                                                                <  2:5112629/150‑1 (MQ=255)
                                                                                          aGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCAc                                                            <  2:1090666/150‑1 (MQ=255)
                                                                                             gcCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTg                                                         >  1:5166634/1‑150 (MQ=255)
                                                                                              ccGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTgg                                                        <  1:1455425/150‑1 (MQ=255)
                                                                                              ccGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTgg                                                        >  2:4585630/1‑150 (MQ=255)
                                                                                              ccGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTgg                                                        <  2:2910742/150‑1 (MQ=255)
                                                                                                   gtggGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTCTGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATc                                                   >  1:2565249/1‑150 (MQ=255)
                                                                                                     gggCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATccc                                                 >  1:497099/1‑150 (MQ=255)
                                                                                                      ggCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATcccc                                                >  2:1912574/1‑150 (MQ=255)
                                                                                                       gCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCa                                               <  1:3825721/150‑1 (MQ=255)
                                                                                                         aGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACt                                             >  2:1059996/1‑150 (MQ=255)
                                                                                                          gCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtg                                            <  1:3038215/150‑1 (MQ=255)
                                                                                                           ccGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtga                                           <  2:512129/150‑1 (MQ=255)
                                                                                                            cgcgcgTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtgat                                          <  1:917366/150‑1 (MQ=255)
                                                                                                             ccgcgTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtgatg                                         >  1:5405142/2‑150 (MQ=255)
                                                                                                              cgcgTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtgatga                                        <  2:211223/150‑1 (MQ=255)
                                                                                                                 gTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACg                                     >  1:4560337/1‑150 (MQ=255)
                                                                                                                     gAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGc                                 >  1:3495825/1‑150 (MQ=255)
                                                                                                                      aGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGcc                                <  1:3682347/150‑1 (MQ=255)
                                                                                                                      aGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGcc                                <  1:5303925/150‑1 (MQ=255)
                                                                                                                           gACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaa                           <  2:1354208/150‑1 (MQ=255)
                                                                                                                           gACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaa                           >  2:1718814/1‑150 (MQ=255)
                                                                                                                              ccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaa                        <  1:1761560/150‑1 (MQ=255)
                                                                                                                              ccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaa                        <  1:3355087/150‑1 (MQ=255)
                                                                                                                              ccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaa                        >  2:5396568/1‑150 (MQ=255)
                                                                                                                               cGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaac                       <  2:3617114/150‑1 (MQ=255)
                                                                                                                               cGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaac                       <  2:5414471/150‑1 (MQ=255)
                                                                                                                               cGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaac                       <  2:5414241/150‑1 (MQ=255)
                                                                                                                                  tgctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTgc                    <  2:5202877/150‑1 (MQ=255)
                                                                                                                                       tgaaCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAg               >  2:833877/3‑150 (MQ=255)
                                                                                                                                         aaCGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtgat                                          >  1:5012550/1‑121 (MQ=255)
                                                                                                                                         aaCGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTc             >  2:2575256/1‑150 (MQ=255)
                                                                                                                                            gcgcCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGa          >  2:294268/1‑150 (MQ=255)
                                                                                                                                               ccTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCt       >  1:3052641/1‑150 (MQ=255)
                                                                                                                                                cTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTg      <  2:3495825/150‑1 (MQ=255)
                                                                                                                                                  ggCCGGCGCGATCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGt    >  2:4196348/1‑150 (MQ=255)
                                                                                                                                                   gCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTg   >  1:3897252/1‑150 (MQ=255)
                                                                                                                                                    ccGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGc  <  1:3445287/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGC  >  NC_002947/2448871‑2449168

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: