Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,690,454 G→A 6.2% T140M (ACG→ATG) ‡ ydcS ← polyamine ABC transporter, periplasmic polyamine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,690,4540GA 6.2% 274.7 / 4.2 94T140M (ACG→ATG) ‡ydcSpolyamine ABC transporter, periplasmic polyamine‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (42/46);  new base A (2/4);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ATCGGCCCGTCATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCGGTGACC  >  NC_002947/1690307‑1690603
                                                                                                                                                   |                                                                                                                                                     
aTCGGCCCGTCATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTg                                                                                                                                                     >  2:1620983/1‑150 (MQ=255)
     cccGTTATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCCTGCTGCTTTTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTa                                                                                                                                                >  2:2300225/1‑150 (MQ=255)
     cccGTTATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCGGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTa                                                                                                                                                >  2:2297091/1‑150 (MQ=255)
     cccGTCATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTa                                                                                                                                                >  1:616404/1‑150 (MQ=255)
        gTCATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTacac                                                                                                                                             >  2:49735/1‑150 (MQ=255)
        gTCATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTacac                                                                                                                                             >  2:49653/1‑150 (MQ=255)
           aTAGGCTTGCACGCGGCCCTTGGTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCGTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCACGCCGTaccagtg                                                                                                                                          <  1:946961/150‑6 (MQ=255)
             aGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCt                                                                                                                                        <  1:2023834/150‑1 (MQ=255)
             aGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCACGCCGTaccagtggg                                                                                                                                        <  2:1027447/150‑8 (MQ=255)
             aGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCACGCCGTaccagtggg                                                                                                                                        <  2:41425/150‑8 (MQ=255)
                cTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGt                                                                                                                                     <  2:1447728/150‑1 (MQ=255)
                 ttGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTc                                                                                                                                    <  2:1383732/150‑1 (MQ=255)
                    cACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGATAGACCTCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGAc                                                                                                                                 <  1:2133231/150‑1 (MQ=255)
                        cgGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACcaca                                                                                                                             <  2:616404/150‑1 (MQ=255)
                        cgGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACcaca                                                                                                                             <  1:1727794/150‑1 (MQ=255)
                           cccTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATAc                                                                                                                          <  1:1288281/150‑1 (MQ=255)
                             cTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTg                                                                                                                                                     >  2:2422530/1‑121 (MQ=255)
                             cTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTg                                                                                                                                                     <  1:2422530/121‑1 (MQ=255)
                             cTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCa                                                                                                                        <  1:1259442/150‑1 (MQ=255)
                             cTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCa                                                                                                                        >  2:841314/1‑150 (MQ=255)
                                 tttGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTa                                                                                                                    <  2:361707/150‑1 (MQ=255)
                                       ttGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGATATCGGTGTTGTACAGCAGCAAGTATGCCCCCCACTCGTAAGGCGCGCAGCACACCGGCCTTTCTACCACATACCAGCTACCGttt                                                                                                              >  2:1205145/1‑149 (MQ=255)
                                       ttGCCGTCGGAACGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTc                                                                                                              >  2:1166693/1‑150 (MQ=255)
                                        tGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCt                                                                                                             <  1:841314/150‑1 (MQ=255)
                                           cGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCa                                                                                                          >  1:2229305/1‑150 (MQ=255)
                                             tCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAgg                                                                                                        <  1:111530/150‑1 (MQ=255)
                                                  cAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCgg                                                                                                   <  1:383599/150‑1 (MQ=255)
                                                   aGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGAc                                                                                                                                 >  1:1570622/1‑119 (MQ=255)
                                                   aGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGAc                                                                                                                                 <  2:1570622/119‑1 (MQ=255)
                                                      ttCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCg                                                                                               >  1:250222/1‑150 (MQ=255)
                                                        cTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTATTTGAAGATATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGAt                                                                                             <  1:1205145/150‑1 (MQ=255)
                                                        cTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCACGTGTAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGAt                                                                                             <  1:1841585/150‑1 (MQ=255)
                                                         tGCGGCTCGAACACCACGCTCCAGCCGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCCCTGTGCCGGCGCGGGTCGATg                                                                                            >  1:721916/1‑150 (MQ=255)
                                                          gCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCACGCCGTacca                                                                                                                                             >  2:739724/1‑98 (MQ=255)
                                                          gCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCACGCCGTacca                                                                                                                                             <  1:739724/100‑3 (MQ=255)
                                                          gCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAACACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTAGATGc                                                                                           >  1:2435521/1‑150 (MQ=255)
                                                                cGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTcc                                                                                     <  1:151196/150‑1 (MQ=255)
                                                                 gAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCa                                                                                    <  1:49653/150‑1 (MQ=255)
                                                                 gAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCa                                                                                    <  1:49735/150‑1 (MQ=255)
                                                                      ccacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTgg                                                                               <  1:1166693/150‑1 (MQ=255)
                                                                      ccacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTgg                                                                               <  2:1141075/150‑1 (MQ=255)
                                                                       cacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATAc                                                                                                                          >  1:936952/1‑106 (MQ=255)
                                                                       cacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATAc                                                                                                                          <  2:936952/106‑1 (MQ=255)
                                                                       cacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATg                                                                                            <  1:1990610/136‑1 (MQ=255)
                                                                       cacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATg                                                                                            >  2:1990610/1‑136 (MQ=255)
                                                                       cacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTggg                                                                              >  1:963953/1‑150 (MQ=255)
                                                                       cacGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTACAGGCGCGGGTCGATGCTCTTCCAGTTggg                                                                              >  2:1974362/1‑150 (MQ=255)
                                                                            tCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCa                                                                         >  1:268066/1‑150 (MQ=255)
                                                                            tCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGAGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCCGTTGGGGATCa                                                                         >  1:760144/1‑150 (MQ=255)
                                                                               aGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAggg                                                                      >  1:1267496/1‑150 (MQ=255)
                                                                               aGCNGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCTACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAggg                                                                      >  2:1508702/1‑150 (MQ=255)
                                                                                 cTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCa                                                                    <  2:1267496/150‑1 (MQ=255)
                                                                                 cTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTGGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTGCAGGTGGGGAGCAGGGCa                                                                    >  1:1074003/1‑150 (MQ=255)
                                                                                        gcgcCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGa                                                             <  1:1914041/150‑1 (MQ=255)
                                                                                         cgcCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGAt                                                            <  1:2211068/150‑1 (MQ=255)
                                                                                          gcCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGCAGTCCAGGCGCGGGTCGATGCTCGTCCAGTTGGGGATCAGGGCAGTGTTGATc                                                           <  1:300047/150‑1 (MQ=255)
                                                                                                 tGAAGATATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGGCCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTa                                                    <  1:136756/150‑1 (MQ=255)
                                                                                                          cGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCa                                                                                    >  1:1251850/1‑109 (MQ=255)
                                                                                                          cGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCa                                                                                    <  2:1251850/109‑1 (MQ=255)
                                                                                                           ggTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTCCACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCgg                                                                                                   >  2:1226732/1‑93 (MQ=255)
                                                                                                           ggTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCgg                                                                                                   <  1:1226732/93‑1 (MQ=255)
                                                                                                           ggTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGcc                                          >  2:672518/1‑150 (MQ=255)
                                                                                                            gtgtTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCt                                         <  1:2430113/150‑1 (MQ=255)
                                                                                                              gttgtACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGc                                       >  1:720419/1‑150 (MQ=255)
                                                                                                                tgtACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCaa                                     <  1:2423446/150‑1 (MQ=255)
                                                                                                                     agcagcACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTc                                                                                                                                    <  2:1568066/50‑1 (MQ=255)
                                                                                                                     agcagcACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTc                                                                                                                                    <  2:1568125/50‑1 (MQ=255)
                                                                                                                     agcagcACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTc                                                                                                                                    >  1:1568066/1‑50 (MQ=255)
                                                                                                                     agcagcACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTc                                                                                                                                    >  1:1568125/1‑50 (MQ=255)
                                                                                                                      gcagcaCGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGcc                               <  1:932168/150‑1 (MQ=255)
                                                                                                                          caCGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAg                           >  2:1198389/1‑150 (MQ=255)
                                                                                                                             gTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGa                        >  2:527730/1‑150 (MQ=255)
                                                                                                                              tttGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTgg                                                                               >  1:2135402/1‑94 (MQ=255)
                                                                                                                              tttGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTgg                                                                               <  2:2135402/94‑1 (MQ=255)
                                                                                                                              tttGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTAGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGAt                       <  2:2435521/150‑1 (MQ=255)
                                                                                                                                tGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGc                     <  2:720419/150‑1 (MQ=255)
                                                                                                                                tGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGc                     <  2:250222/150‑1 (MQ=255)
                                                                                                                                    ccccACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCg                 >  1:559198/1‑150 (MQ=255)
                                                                                                                                    ccccACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCg                 >  1:1468633/1‑150 (MQ=255)
                                                                                                                                    ccccACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCg                 >  1:559450/1‑150 (MQ=255)
                                                                                                                                    ccccACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGAAGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCg                 >  1:2303462/1‑150 (MQ=255)
                                                                                                                                     cccACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTgg                                                                               >  1:1768248/1‑87 (MQ=255)
                                                                                                                                     cccACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTgg                                                                               <  2:1768248/87‑1 (MQ=255)
                                                                                                                                     cccACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTACAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGc                >  1:2385203/1‑150 (MQ=255)
                                                                                                                                          gtgTACGGCACGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCAGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGCGTTGATCGGCTGTACGCGCTTGCCTGCACTAAGCCGAAGCGCTGCATCGCCCGAg           >  2:1695802/3‑150 (MQ=255)
                                                                                                                                             tACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCaa                                     >  2:1422873/1‑121 (MQ=255)
                                                                                                                                             tACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCaa                                     <  1:1422873/121‑1 (MQ=255)
                                                                                                                                             tACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCg        <  2:1468633/150‑1 (MQ=255)
                                                                                                                                             tACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCg        <  1:998988/150‑1 (MQ=255)
                                                                                                                                               cGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTTCCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCGGt      <  1:1643864/150‑1 (MQ=255)
                                                                                                                                                ggCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCTGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCGGTg     <  2:1371059/150‑1 (MQ=255)
                                                                                                                                                ggCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCGGTg     >  1:2074241/1‑150 (MQ=255)
                                                                                                                                                 gCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCt                                                       <  2:2229218/99‑1 (MQ=255)
                                                                                                                                                 gCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCt                                                       >  1:2229218/1‑99 (MQ=255)
                                                                                                                                                   gTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCGGTGAcc  >  2:877659/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATCGGCCCGTCATAGGCTTGCACGCGGCCCTTGTTTGGCTTGCCGTCGGGCAGGTTCTGCGGCTCGAACACCACGCTCCAGCTGGTGGGCGCCTGCTTGAAGGTATCGGTGTTGTACAGCAGCACGTTTGGCCCCCACTGGTACGGCGTGCCGTACACCTGCTTGTCGACCACATACCAGCTACCGTTCTGCAGGCGCGGGTCGATGCTCTTCCAGTTGGGGATCAGGGCAGTGTTGATCGGCTGTACGCGCTTGCCTGCAATCAGCCGCAGCGATGCATCGCCCGAGGCGGTGACC  >  NC_002947/1690307‑1690603

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: