Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 292,607 C→G 6.9% intergenic (+22/‑38) ssuA → / → ssuD aliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947292,6070CG 6.9% 223.6 / 7.2 72intergenic (+22/‑38)ssuA/ssuDaliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base G (2/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGA  >  NC_002947/292459‑292747
                                                                                                                                                    |                                                                                                                                            
cccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgcc                                                                                                                                             >  1:377853/1‑150 (MQ=255)
 ccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCTCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgccg                                                                                                                                            >  1:615290/1‑150 (MQ=255)
  cAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgccgc                                                                                                                                           >  2:786428/1‑150 (MQ=255)
  cAGAGGTGGTGACCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgccgc                                                                                                                                           >  1:77911/1‑150 (MQ=255)
    gagGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGccc                                                                                                                                         <  2:258908/150‑1 (MQ=255)
       gtggtgGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCgg                                                                                                                                      <  1:239814/150‑1 (MQ=255)
       gtggtgGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCgg                                                                                                                                      <  2:377853/150‑1 (MQ=255)
            ggCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcactg                                                                                                                                 <  1:472066/150‑5 (MQ=255)
             gCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTagag                                                                                                                                >  1:872170/1‑150 (MQ=255)
               cgcgCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGcc                                                                                                                              >  1:1015165/1‑150 (MQ=255)
                gcgcAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCa                                                                                                                             >  1:527825/1‑150 (MQ=255)
                    aaCAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCGCCGGATAGAGACACCCt                                                                                                                         >  2:491635/1‑150 (MQ=255)
                      cAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTg                                                                                                                       <  1:424635/150‑1 (MQ=255)
                       agaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                      <  1:710892/150‑1 (MQ=255)
                        gaagaTCGCCGCCACCTTACCGGCGCTCAAGCTGTTACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                     <  2:615290/150‑1 (MQ=255)
                          agaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                   >  1:250472/1‑150 (MQ=255)
                          agaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                   <  2:448273/150‑1 (MQ=255)
                          agaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                   <  2:1015165/150‑1 (MQ=255)
                                    cacCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCa                                                                                                         >  2:498629/1‑150 (MQ=255)
                                     acCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                      <  2:961451/136‑1 (MQ=255)
                                     acCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                      >  1:961451/1‑136 (MQ=255)
                                       cTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                 >  2:1154721/1‑123 (MQ=255)
                                       cTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                 <  1:1154721/123‑1 (MQ=255)
                                             ggCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCtt                                                                                                >  2:488497/1‑150 (MQ=255)
                                             ggCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCtt                                                                                                >  2:720530/1‑150 (MQ=255)
                                                gcTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAAc                                                                                             <  2:743886/150‑1 (MQ=255)
                                                gcTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAAc                                                                                             <  2:872170/150‑1 (MQ=255)
                                                 cTCAAGCTGATACCCAAGCCGCTGAGCATCAAGAACATGATAAGAAAATGTTACGTCAAATCCACTAGCGCCCCGTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACa                                                                                            >  1:626267/1‑150 (MQ=255)
                                                       cTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCt                                                                                      <  1:498629/150‑1 (MQ=255)
                                                        tGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                     >  2:738745/1‑150 (MQ=255)
                                                        tGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                     <  1:639785/150‑1 (MQ=255)
                                                         gATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCt                                                                                      >  1:719489/1‑148 (MQ=255)
                                                         gATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCt                                                                                      <  2:719489/148‑1 (MQ=255)
                                                          aTACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGt                                                                                   <  1:720530/150‑1 (MQ=255)
                                                             cccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcc                                                                                <  2:527825/150‑1 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGGGCCGGGGCGCCGCCCCGGCTAGAGCCACCGTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcct                                                                               >  1:683365/1‑150 (MQ=255)
                                                               cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctc                                                                              >  2:7809/1‑150 (MQ=255)
                                                                aaGCCGCTGAGCATCAAGGACGTGATCTGCACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctcc                                                                             <  1:491635/150‑1 (MQ=255)
                                                                          gCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACgg                                                                   <  1:488497/150‑1 (MQ=255)
                                                                           cATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGt                                                                  >  2:377262/1‑150 (MQ=255)
                                                                           cATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGt                                                                  >  2:722199/1‑150 (MQ=255)
                                                                             tCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGc                                                                                                   >  1:477561/1‑115 (MQ=255)
                                                                             tCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGc                                                                                                   <  2:477561/115‑1 (MQ=255)
                                                                              cAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGAc                                                               <  2:77911/150‑1 (MQ=255)
                                                                                aGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCgg                                                                                                                                      <  2:297675/77‑1 (MQ=255)
                                                                                aGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCgg                                                                                                                                      >  1:297675/1‑77 (MQ=255)
                                                                                 ggACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATg                                                                                                      >  1:167341/1‑108 (MQ=255)
                                                                                 ggACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATg                                                                                                      <  2:167341/108‑1 (MQ=255)
                                                                                     gTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                     >  1:376638/1‑121 (MQ=255)
                                                                                     gTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                     <  2:376638/121‑1 (MQ=255)
                                                                                     gTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGt                                                        >  1:600858/1‑150 (MQ=255)
                                                                                         tCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgacga                                                                                                              <  2:79736/92‑1 (MQ=255)
                                                                                         tCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgacga                                                                                                              >  1:79736/1‑92 (MQ=255)
                                                                                             gACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGc                                                <  1:377262/150‑1 (MQ=255)
                                                                                             gACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGc                                                >  1:123371/1‑150 (MQ=255)
                                                                                                   cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAAt                                                                                                       >  1:721631/1‑89 (MQ=255)
                                                                                                   cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAAt                                                                                                       <  2:721631/89‑1 (MQ=255)
                                                                                                   cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTcc                                          >  2:640609/1‑150 (MQ=255)
                                                                                                    ccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCg                                         >  1:740470/1‑150 (MQ=255)
                                                                                                       gCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAg                                      <  1:786428/150‑1 (MQ=255)
                                                                                                        ccAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGCAgg                                     >  2:183580/1‑150 (MQ=255)
                                                                                                         cAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACAtttt                                                                                        <  1:1007333/98‑1 (MQ=255)
                                                                                                         cAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACAtttt                                                                                        >  2:1007333/1‑98 (MQ=255)
                                                                                                                gCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTg                             >  2:183526/1‑150 (MQ=255)
                                                                                                                            tgattgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                 <  1:183526/150‑1 (MQ=255)
                                                                                                                            tgattgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCATTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                 <  2:402565/150‑1 (MQ=255)
                                                                                                                                 gaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTAc            <  2:626267/150‑1 (MQ=255)
                                                                                                                                    tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTACCATTTTCTGGTTCCTCCCCACCCCCCGTGACGGCAAGTACCTGGCCACCTCCGAAGGCCCACGCGCGGTCCACCACGGTTACCtc         >  2:543366/1‑149 (MQ=255)
                                                                                                                                    tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTg         <  1:722199/150‑1 (MQ=255)
                                                                                                                                       tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTAc            >  1:925252/5‑144 (MQ=255)
                                                                                                                                       tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTAc            <  2:925252/140‑1 (MQ=255)
                                                                                                                                       tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgc        <  1:48653/144‑1 (MQ=255)
                                                                                                                                       tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgc        >  2:48653/5‑148 (MQ=255)
                                                                                                                                       cAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAGCGTTTTCTGGTTCCTCCCCACCCACGGGGACGGCAAGTCACTGGGCACTTCCGAAGGCGCACGTGCGGTCGACGACGGTTACCTgcag      >  2:488903/1‑150 (MQ=255)
                                                                                                                                       cAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgcag      <  1:738745/150‑1 (MQ=255)
                                                                                                                                           gCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGa  <  2:250472/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGA  >  NC_002947/292459‑292747

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