Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,511,474 T→G 7.5% M435R (ATG→AGG)  PP_4839 → iron‑regulated membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,511,4740TG 7.5% 503.4 / 2.8 216M435R (ATG→AGG) PP_4839iron‑regulated membrane protein
Reads supporting (aligned to +/- strand):  ref base T (62/121);  new base G (3/20);  total (65/151)
Fisher's exact test for biased strand distribution p-value = 5.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

TTCGGTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCC  >  NC_002947/5511325‑5511617
                                                                                                                                                     |                                                                                                                                               
 tCGGTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtg                                                                                                                                                <  2:598749/150‑1 (MQ=255)
   ggTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggt                                                                                                                                              >  1:5651753/1‑150 (MQ=255)
    gTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggtg                                                                                                                                             <  1:5478433/150‑1 (MQ=255)
    gTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggtg                                                                                                                                             >  1:6346930/1‑150 (MQ=255)
    gTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTGGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggtg                                                                                                                                             <  2:714754/150‑1 (MQ=255)
    gTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCATGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggtg                                                                                                                                             <  2:3850279/150‑1 (MQ=255)
     tGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCGTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGAAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggtgg                                                                                                                                            <  2:4442816/150‑1 (MQ=255)
     tGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggtgg                                                                                                                                            <  1:4453335/150‑1 (MQ=255)
      gAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATggtggt                                                                                                                                           <  2:3946851/150‑1 (MQ=255)
       aGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGgtggtg                                                                                                                                          <  1:684301/150‑1 (MQ=255)
       aGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGgtggtg                                                                                                                                          <  1:6091811/150‑1 (MQ=255)
       aGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGgtggtg                                                                                                                                          <  1:1717205/150‑1 (MQ=255)
         cGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAt                                                                                                                                        <  2:5618984/150‑1 (MQ=255)
           ggCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAtgt                                                                                                                                      <  1:5844981/150‑1 (MQ=255)
           ggCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAtgt                                                                                                                                      <  1:1904905/150‑1 (MQ=255)
            gCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAtgtg                                                                                                                                     >  1:722689/1‑150 (MQ=255)
             cTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTgg                                                                                                                                    <  2:906872/150‑1 (MQ=255)
              tGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTggg                                                                                                                                   >  1:3592293/1‑150 (MQ=255)
              tGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTggg                                                                                                                                   >  2:172045/1‑150 (MQ=255)
               ggTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTgggg                                                                                                                                  <  1:4956412/150‑1 (MQ=255)
                    tgtgGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTgg                                                                                                                             <  2:3489412/150‑1 (MQ=255)
                         tggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGTGGATAGtt                                                                                                                       <  2:3510360/150‑1 (MQ=255)
                          ggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCTGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGtt                                                                                                                       <  1:3045417/150‑1 (MQ=255)
                          ggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGTTGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGtt                                                                                                                       <  1:2308547/150‑1 (MQ=255)
                           gAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGttt                                                                                                                      <  1:961666/150‑1 (MQ=255)
                           gAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGttt                                                                                                                      >  1:2853699/1‑150 (MQ=255)
                              gcgcCGTCGTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggt                                                                                                                   <  1:6777359/150‑1 (MQ=255)
                              gcgcCGTCCTGAATGCGGCCTGGGCGTACCGCCTTTGCGTCTTGACCTGCCACGCTGGCAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggt                                                                                                                   <  1:55277/150‑1 (MQ=255)
                              gcgcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCTGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggt                                                                                                                   <  2:2652540/150‑1 (MQ=255)
                              gcgcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggt                                                                                                                   <  2:2770794/150‑1 (MQ=255)
                              gcgcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCGGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggt                                                                                                                   <  2:138831/150‑1 (MQ=255)
                               cgcCGTCCTGAAGGCGGCCTGGGCGTCCCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggtg                                                                                                                  <  2:2080209/150‑1 (MQ=255)
                               cgcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCTGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGCTGGTGGTGATGTGGGGGGTGGATAGTTtggtg                                                                                                                  <  2:6229929/150‑1 (MQ=255)
                               cgcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggtg                                                                                                                  <  1:3254004/150‑1 (MQ=255)
                                gcCGTCCTGATGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggtgg                                                                                                                 <  1:5042618/150‑1 (MQ=255)
                                gcCGTCCTGAAGGCGGCCTGGGCGTCCCGCCTTTGCGTCCTGCCCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggtgg                                                                                                                 <  1:6161891/150‑1 (MQ=255)
                                   gTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTTGGGGGTGGATAGTTTGGTGGTGa                                                                                                              <  1:5073326/150‑1 (MQ=255)
                                    tCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAg                                                                                                             <  2:4745197/150‑1 (MQ=255)
                                     ccTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAgg                                                                                                            <  2:395322/150‑1 (MQ=255)
                                      cTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGgc                                                                                                           <  2:6862967/150‑1 (MQ=255)
                                            gcggcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCgccgcc                                                                                                     >  1:3636114/1‑150 (MQ=255)
                                            gcggcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGGGGTGAGGGGGGGGGTGGATCGGGTGGTGGTGAGGCgccgcc                                                                                                     >  1:6393428/1‑150 (MQ=255)
                                            gcggcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGGGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGGGGGGGGGTCGATGGGGGGGATGTGGGGGGGGGATCGGTTGGTGGTGAGGCgccgcc                                                                                                     >  2:2923159/1‑150 (MQ=255)
                                               gcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCCGGGGGTGATGTGCCCGCTGGTGGGGGTGTCGCTGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCggg                                                                                                  <  2:2463112/150‑1 (MQ=255)
                                                ccTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAtggtg                                                                                                                                             >  1:1006273/1‑106 (MQ=255)
                                                    ggCGTACCGCCTTTGCGTCATGACCTGCCACGCGGGAAGGCGGCAGTGGCCGTGATGCGGGTGCTGGGGGTGATGTTCCCGCGGGTGGGGGTGCCGAGGGTGGTGCTGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTg                                                                                             <  2:4696297/150‑1 (MQ=255)
                                                    ggCGGACCGCCTGTGCGTCAGGCCCCGCCCCGCTGGAAGGCGGCAGGGGCCGTGATGCTGGCGCTGGGGGTGATGCTCCCGCGGGCGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTg                                                                                             <  1:4168265/150‑1 (MQ=255)
                                                     cactgccGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTg                                                                                                                              <  2:4260120/111‑1 (MQ=255)
                                                      cGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGc                                                                                           <  2:3636114/150‑1 (MQ=255)
                                                      cGTACCGCCCTCGCGTCAGGCCCCGCCCCGCTGGCCGGCGGCAGTGGCCGTGCCGCTGGGGCGGGGGGCGCGGGGCCCGCGGGGGGGGGTGTCGAGGGTGGTGACGCGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGc                                                                                           <  1:3073708/150‑1 (MQ=255)
                                                        tACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCGGGTGCGGGGGGTGAGGGTCCCGCTGGTGGGGGTGTCGATGGTGGCGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCa                                                                                         <  2:1478007/150‑1 (MQ=255)
                                                              cTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCtg                                                                                   >  2:185742/1‑149 (MQ=255)
                                                                ttGCGTCATGACCTGCCACGCTGGACGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGCGATGTTCCCGCTGGTGGGGGTGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGa                                                                                 <  2:2488010/150‑1 (MQ=255)
                                                                    gTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCt                                                                             >  1:2593464/1‑150 (MQ=255)
                                                                    gTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCt                                                                             >  1:765569/1‑150 (MQ=255)
                                                                     tcctGACCTGCCACGCTGGCAGGCGGCAGTGGCCGTGATGCTGGCGCGGGGGGTGCTGTGCCCGCTGGGGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCtt                                                                            <  2:1081794/147‑1 (MQ=255)
                                                                     gcATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGGTGGGGGTGATGTGCCCGCTGGTGGGGGTGTCGATGGTGCTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCGtt                                                                            <  2:4827922/149‑1 (MQ=255)
                                                                       aggACCTGCCACGCTGGCAGGCGGCAGTGGCCGTGCTGCTGGTGCTGGGGGTGATGGTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTtgt                                                                          <  2:6072387/148‑1 (MQ=255)
                                                                          cccggCCCCGCGGGACGGCGGCCGGGGCCGGGCTGCGGGGGCCGGGGGGGCGGGCCCCGCTGGGGGGGGCGGCGAGGGTGGGGCTGTGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGtt                                                                       <  2:1819155/146‑1 (MQ=255)
                                                                          cccggCCACGCGGGAAGGCGGCAGTGGCCGTGATGCTGGTGCGGGGGGGGATGGTCCCGCTGGGGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGtt                                                                       <  2:5768233/146‑1 (MQ=255)
                                                                            cTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGc                                                                     >  2:513462/1‑150 (MQ=255)
                                                                                cACGCTGGACGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCtgt                                                                 <  2:5193178/150‑1 (MQ=255)
                                                                                    cTGGCAGGCGGCCGTGGCCGTGATGCTGGTGCTGGGGGTGAGGTTCCCGCTGGTGGGGGTGTCGACGGGGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTg                                                             <  2:5045969/150‑1 (MQ=255)
                                                                                     tGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTCGGGGTGAGGTTCCCGCAGGTGGGGGTCTCGAGCCTGGTGATGTGGGGCGTGGATAGTTTGGGGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTCTGCAgg                                                            >  2:6269038/1‑150 (MQ=255)
                                                                                      ggcaGGCGGCAGGGGCCGCGCCGCGGGCGCCGGGGGGGCTGGTCCCGCCGGCGGGGGGGTCGATGGGGGCGACGCGGCGGCGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGc                                                           <  1:1360523/147‑1 (MQ=255)
                                                                                           gcggcAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCtt                                                      >  1:3708559/1‑150 (MQ=255)
                                                                                                cggggCCGCGATCCGGGCGCGGGGGGTGCCGGCCCCGCTGGTGGGGGGGTCGAGGGCGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCgg                                                 <  2:3887665/147‑1 (MQ=255)
                                                                                                cgcggCCGTGCCGCTGGAGCCGGGGGCGCTGCTCCCGCTGGCGGGCGTGTCGACGGCGGTGCCGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCgg                                                 <  1:629160/147‑1 (MQ=255)
                                                                                                cgTGGCCGTGACGCGGGTGCGGGGGGTGACGCTCCCGCGGGGGGGGGCCTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCgg                                                 <  2:3297468/149‑1 (MQ=255)
                                                                                                aGTGGCCGTGCTGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCgg                                                 <  2:3163758/150‑1 (MQ=255)
                                                                                                   ggCCGTGAGGCTGGTGCTGGGGGTGACGTGCCCGCTGGTGGGGGGGTCGAGGGTGGGGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTa                                              <  2:257228/150‑1 (MQ=255)
                                                                                                      cGTGATGCTGGGGCGGGGGGTGATGGTCCCGCGGGGGGGGGGGTCGATGGTGGTGATGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                           <  2:715843/150‑1 (MQ=255)
                                                                                                      cGTGATGCTGCCGCCGGGGCCGATCCCCCCGCGGGCGGGCCTGACGCCGGTGGCGAGGCGGGGGGTGGATAGTTTGGTGGTGCGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                           <  1:2240793/150‑1 (MQ=255)
                                                                                                       gtgcggctggtgctggGGGCGATGGTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAcc                                          <  2:5498970/145‑1 (MQ=255)
                                                                                                       gTGATGCTGGCGCGGGGGGTGAGGCCCCCGCGGGGGGGGGAGACGATGGTGGCGCTGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAcc                                          <  1:1147668/150‑1 (MQ=255)
                                                                                                          aTGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTggtggt                                                                                                                >  1:3555644/1‑77 (MQ=255)
                                                                                                                   ctggGGGTGCCGCCCCCGCGGGCGGGGGCGTCGCCGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaac                              <  1:2427067/150‑4 (MQ=37)
                                                                                                                   ctggGGGTGATGGTCCCGCGGGGGGGGGGGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAg                              <  2:6440314/150‑1 (MQ=39)
                                                                                                                   ccggggggGCTGTGCCCGCGGGGGGGGGCGGCGAGGGGGGCGATGCGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaac                              <  1:4023042/148‑4 (MQ=25)
                                                                                                                          cgccggccccGCTGGGGGGGGCGCCGAGGGGGGGGCGGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACcgc                       <  2:4284960/143‑1 (MQ=37)
                                                                                                                             tGTTCCCGCTGGTGGGGGGGTCGAGGGTGGCGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGtt                    <  2:1884207/150‑1 (MQ=255)
                                                                                                                             cgTTCCCGCTGGGGGGGGCGTCGCTGGTCGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGtt                    <  2:4152773/149‑1 (MQ=255)
                                                                                                                              gTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTg                   >  2:153099/1‑150 (MQ=255)
                                                                                                                               ttCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCACCAACGAGGGACcgccgtgt                  >  2:659796/1‑145 (MQ=255)
                                                                                                                               ttCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGGAAAACCCCACCCAACAGGGACCCCCTTGt                  >  1:4588053/1‑150 (MQ=255)
                                                                                                                               ttCCCGCTGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCGCCTACGAGGGACcgcggggt                  >  2:1020598/1‑146 (MQ=255)
                                                                                                                               gccccGCGGGCGGGGCGGCCGAGGGCGGCGATGGGGCGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGt                  <  1:895735/148‑1 (MQ=37)
                                                                                                                                tCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTa                 <  2:4530942/150‑1 (MQ=255)
                                                                                                                                gcccGCGGGGGGGGGCGGCGCGGGTGGTGCGGTGGGGGGGGGATAGTTTGGGGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGGTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTa                 <  2:861362/149‑1 (MQ=37)
                                                                                                                                gcccGCCGGGGGGGGGGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTa                 <  2:5906242/149‑1 (MQ=255)
                                                                                                                                 cccGCGGGGGGGGGGGTCGACGGTGGCGATGGGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAg                <  2:3708559/150‑1 (MQ=38)
                                                                                                                                  ccgccgGCGGGGGTGTCGATGGCGCCGATGCGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACACGCTCCTACGAGGGACCGCGTTGTAgg               <  2:1604527/150‑1 (MQ=38)
                                                                                                                                      tgggggGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCg           <  2:6748689/150‑1 (MQ=255)
                                                                                                                                      tggcggGGGTGTCGATGGTGGTGCTGCGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCg           <  2:4527557/150‑1 (MQ=255)
                                                                                                                                      gggggggggTGGCGAGGGGGGAGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCg           <  2:132243/146‑1 (MQ=38)
                                                                                                                                       ggtggGGGCGCCGCTGGGGGTGAGGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg          <  2:1419481/150‑1 (MQ=38)
                                                                                                                                       ggcgggggCGCCGATGGGGGTGACGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg          <  2:2355456/147‑1 (MQ=38)
                                                                                                                                        gggggggTGCCGCGGGTGGTGCTGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCggg         <  2:1277204/148‑1 (MQ=255)
                                                                                                                                            gggTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCCGGTAAACCCCCTCCAACGACGGACCCCGTTGTAGGAGCGGGTTTa     >  1:1011906/1‑150 (MQ=255)
                                                                                                                                               tgtCGATGGAGGCGATGCGGGGGGCGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTAccc  <  2:2566760/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TTCGGTGAGCGGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCC  >  NC_002947/5511325‑5511617

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