Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,645,153 T→C 5.2% E735G (GAG→GGG)  PP_2316 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,645,1530TC 5.2% 228.5 / 2.2 78E735G (GAG→GGG) PP_2316ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base C (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

CAGCACCCAGGTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAGGCATCCACCTGCAGGATGGTCACCG  >  NC_002947/2645008‑2645300
                                                                                                                                                  |                                                                                                                                                    
cAGCACCCAGGTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGTTGGCCGCACCCAGCAGACCGAACTCGATTCGCCGGGCCTTGACCAGCACTACGCGGCCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGGGCTCATc                                                                                                                                                  >  1:1574231/1‑150 (MQ=255)
  gCACCCAGTTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCa                                                                                                                                                <  2:1458539/150‑1 (MQ=255)
    aCCCAGGTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGc                                                                                                                                              >  1:1674260/1‑150 (MQ=255)
      ccAGGTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGt                                                                                                                                            >  1:637703/1‑150 (MQ=255)
      ccAGGTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGt                                                                                                                                            >  1:1571065/1‑150 (MQ=255)
      ccAGGTAATCAGTTCGCAGCCCACGGCCCCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCCGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGt                                                                                                                                            >  1:261035/1‑150 (MQ=255)
          gTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGcc                                                                                                                                        >  1:1172508/1‑150 (MQ=255)
          gTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGACCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGcc                                                                                                                                        <  1:483202/150‑1 (MQ=255)
              tCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGTCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCa                                                                                                                                    >  2:1355948/1‑150 (MQ=255)
               cAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCaa                                                                                                                                   >  1:1599384/1‑150 (MQ=255)
                    cGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCgg                                                                                                                              <  1:598145/150‑1 (MQ=255)
                        gCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTACGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGaa                                                                                                                          <  1:1876859/150‑1 (MQ=255)
                           cACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCTCGCAACAGTCCCCCTTGGAGACTGCGCTCATCCAGCGTCGCCTGCAAGCCGGCGAACAg                                                                                                                       <  1:536131/150‑1 (MQ=255)
                            aCGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGc                                                                                                                      >  1:601744/1‑150 (MQ=255)
                                 cgccAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGc                                                                                                                 <  2:1571065/150‑1 (MQ=255)
                                 caacAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGc                                                                                                                 <  2:268370/147‑1 (MQ=255)
                                  gctaagaGCCCGCTGGCCCTACCCGAGAGGCCGAACGCGATTCGCCGTGCCTTAACCAGCAACGGCCGCGCCGCGCCCACTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGcc                                                                                                                <  2:734328/144‑1 (MQ=255)
                                    cAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCTAACAGCACTGCCAg                                                                                                              <  1:755884/150‑1 (MQ=255)
                                          cccGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCcggc                                                                                                        >  1:136386/1‑150 (MQ=255)
                                               tGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCa                                                                                                  <  2:1172508/150‑1 (MQ=255)
                                                  ccGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCa                                                                                               <  2:961697/150‑1 (MQ=255)
                                                   cGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCAc                                                                                              <  2:756342/150‑1 (MQ=255)
                                                   cGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCAc                                                                                              >  1:1310942/1‑150 (MQ=255)
                                                      aCCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGaa                                                                                           >  1:1648707/1‑150 (MQ=255)
                                                      aCCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGaa                                                                                           >  2:869080/1‑150 (MQ=255)
                                                        ccAGCAGGCCGAACTCTATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATCCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACa                                                                                         >  1:1269802/1‑150 (MQ=255)
                                                        ccAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACa                                                                                         >  2:889213/1‑150 (MQ=255)
                                                            cagGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCaa                                                                                     <  1:889213/150‑1 (MQ=255)
                                                             agGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAac                                                                                    <  2:261035/150‑1 (MQ=255)
                                                             agGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAac                                                                                    >  2:1893629/1‑150 (MQ=255)
                                                              ggCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAaca                                                                                   <  2:914423/150‑1 (MQ=255)
                                                               gCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGc                                                                                                                             <  1:1918290/108‑1 (MQ=255)
                                                               gCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGc                                                                                                                             >  2:1918290/1‑108 (MQ=255)
                                                               gCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGCCCGCACCGACACCCACGCGCACACTCGTGCCCCTTGGCGAGTTCTTTCATCCAGCGGGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAacac                                                                                  <  2:1674260/150‑1 (MQ=255)
                                                                ccGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAAcaca                                                                                 <  2:42687/150‑1 (MQ=255)
                                                                ccGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAAcaca                                                                                 >  1:1093107/1‑150 (MQ=255)
                                                                ccGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAAcaca                                                                                 <  1:38371/150‑1 (MQ=255)
                                                                      tCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTcc                                                                           >  1:827433/1‑150 (MQ=255)
                                                                      tCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTcc                                                                           >  2:1101022/1‑150 (MQ=255)
                                                                        gATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCAc                                                                         >  2:1759140/1‑150 (MQ=255)
                                                                         aTTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACt                                                                                                                   <  1:1223093/108‑1 (MQ=255)
                                                                         aTTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACt                                                                                                                   >  2:1223093/1‑108 (MQ=255)
                                                                          ttCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCg                                                                       <  2:136386/150‑1 (MQ=255)
                                                                              ccGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAg                                                                   <  1:1101022/150‑1 (MQ=255)
                                                                               cGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAgg                                                                  <  2:1679692/150‑1 (MQ=255)
                                                                                    cTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTAc                                                             <  1:1820641/150‑1 (MQ=255)
                                                                                        aCCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTg                                                                                                                                           <  1:59786/69‑1 (MQ=255)
                                                                                        aCCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTg                                                                                                                                           >  2:59786/1‑69 (MQ=255)
                                                                                           aGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGcgcg                                                      >  1:495841/1‑150 (MQ=255)
                                                                                            gCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTTCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCga                                                     >  1:941556/1‑150 (MQ=255)
                                                                                            gCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGt                                                                <  1:875686/139‑1 (MQ=255)
                                                                                            gCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGt                                                                >  2:875686/1‑139 (MQ=255)
                                                                                               aCGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGaa                                                  <  2:290268/150‑1 (MQ=255)
                                                                                               aCGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGaa                                                  >  1:1025058/1‑150 (MQ=255)
                                                                                                 ggCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAAtg                                                <  2:601744/150‑1 (MQ=255)
                                                                                                  gccgcgccgcGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAAtgc                                               >  1:432214/1‑150 (MQ=255)
                                                                                                  cccgcgccgcgCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAAtgt                                               >  1:1827540/2‑149 (MQ=255)
                                                                                                         cgcgccCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGGACCGACATGAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAg                                        >  2:1021538/1‑150 (MQ=255)
                                                                                                          gcgccCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGc                                       >  1:1786474/1‑150 (MQ=255)
                                                                                                            cgccCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCt                                      <  1:1355948/150‑1 (MQ=255)
                                                                                                               ccAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCCCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCcggcg                                                                                                       >  1:161629/1‑83 (MQ=255)
                                                                                                               ccAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCCCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCcggcg                                                                                                       <  2:161779/83‑1 (MQ=255)
                                                                                                               ccAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCCCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCcggcg                                                                                                       <  2:161629/83‑1 (MQ=255)
                                                                                                               ccAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCCCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCcggcg                                                                                                       >  1:161779/1‑83 (MQ=255)
                                                                                                               caaGTGCGCGCAACAGTGCTCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGt                                   >  1:630371/3‑150 (MQ=255)
                                                                                                                cAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTc                                  <  2:1599384/150‑1 (MQ=255)
                                                                                                                      gcgcAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCaa                            <  2:637703/150‑1 (MQ=255)
                                                                                                                       cgcAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAg                           <  2:1648707/150‑1 (MQ=255)
                                                                                                                       cgcAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAg                           <  2:1025058/150‑1 (MQ=255)
                                                                                                                        gcAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAgg                          >  2:1860884/1‑150 (MQ=255)
                                                                                                                            cAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAGGCATc                      <  2:1093107/150‑1 (MQ=255)
                                                                                                                                   ccTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTAGTCGAGCAAGGCATCCACCTGc               >  2:1699810/1‑150 (MQ=255)
                                                                                                                                    cTTGGCGCATGCGCGGAGCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGACAGCAGGGGTGCGTAGCGGATGGCGAACACATACTCCGCCGACAGGGTTACCTGCGCGGGCCAGCTGCGCAGTTGGTCGAGCAAGGCATCCACCTGCa              <  1:1893629/150‑1 (MQ=255)
                                                                                                                                    cTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAGGCATCCACCTGCa              <  2:1310942/150‑1 (MQ=255)
                                                                                                                                     ttGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAATAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAGGCATCCACCTGCAg             <  2:1695912/150‑1 (MQ=255)
                                                                                                                                                gcTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTa                                                              >  1:998625/1‑90 (MQ=255)
                                                                                                                                                gcTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTa                                                              <  2:998625/90‑1 (MQ=255)
                                                                                                                                                gcTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAGGCATCCACCTGCAGGATGGTCAccg  >  2:382449/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGCACCCAGGTAATCAGTTCGCAGCCCACGGCCGCCAATAGCCCGCTGGCCGCACCCAGCAGGCCGAACTCGATTCGCCGGGCCTTGACCAGCAACGGCCGCGCCGCGCCCAGTGCGCGCAACAGTGCCCCTTGGCGAATGCGCTCATCCAGCGTGGCCTGCAAGCCGGCGAACAGCACTGCCAGCCCGGCGGCCAGCACGAACAGCAACACATACTCCACCGCCAGGGTTACCTGCGCGAGAATGCTGCGCAGCTGGTCGAGCAAGGCATCCACCTGCAGGATGGTCACCG  >  NC_002947/2645008‑2645300

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