Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 292,607 C→G 6.6% 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.6% 189.4 / 4.7 61intergenic (+22/‑38)ssuA/ssuDaliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (24/33);  new base G (2/2);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAGG  >  NC_002947/292459‑292756
                                                                                                                                                     |                                                                                                                                                     
cccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCGAGCTCGCCTGGATGGACTGCCAGTGCCGCGGcgcc                                                                                                                                                      >  2:1366053/1‑150 (MQ=255)
 ccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgccg                                                                                                                                                     >  2:1040333/1‑150 (MQ=255)
   agagGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCgccgcc                                                                                                                                                   >  1:1262705/1‑150 (MQ=255)
          gtggCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGTACGTGATCTGGACACCCCCTGCCAAAGTCGCTATCGCGCCTTGATTGGCGTCCAGTGCCTGGGCGCCGCCCCGGCTa                                                                                                                                            >  2:299723/1‑150 (MQ=255)
                     aCAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCtt                                                                                                                                 <  1:299723/150‑1 (MQ=255)
                        gaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                              <  1:990259/150‑1 (MQ=255)
                        gaagaTCGCCGACACCTTCCAGGCGCCACCGATGATCCATAGCTCGTGGAGCATCAAGGACGTGATCTGGACACCGCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                              <  1:1027015/150‑1 (MQ=255)
                            aTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGgaga                                                                                                                          <  2:1041834/150‑1 (MQ=255)
                            aTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGgaga                                                                                                                          <  1:204006/150‑1 (MQ=255)
                             tCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCTCCGGCTAGAGCCACCCTTGAGGAgac                                                                                                                         <  1:750293/150‑1 (MQ=255)
                               gCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgacga                                                                                                                       <  2:1262705/150‑1 (MQ=255)
                                        ttCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                            >  2:1353309/1‑136 (MQ=255)
                                        ttCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                            <  1:1353309/136‑1 (MQ=255)
                                        ttCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGa                                                                                                              >  1:931947/1‑150 (MQ=255)
                                                  tCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGTCTCCAATGAGCCTTAACAt                                                                                                    <  2:1163120/150‑1 (MQ=255)
                                                    aaGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACAttt                                                                                                  >  1:1396665/1‑150 (MQ=255)
                                                    aaGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTCGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACAttt                                                                                                  >  1:88984/1‑150 (MQ=255)
                                                      gCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcact                                                                                                                                           >  2:1079833/1‑104 (MQ=255)
                                                      gCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcact                                                                                                                                           <  1:1079833/107‑4 (MQ=255)
                                                        tGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCt                                                                                               <  1:45135/149‑1 (MQ=255)
                                                        tGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCt                                                                                               >  2:45135/1‑149 (MQ=255)
                                                          aTACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGCCTAGCGATGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGGCGACTCCACTGAGCCTTAACATCTTCTGGt                                                                                            >  1:1026606/1‑150 (MQ=255)
                                                            aCCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtc                                                                                          <  1:590663/150‑1 (MQ=255)
                                                              ccgaGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTACCATTTTCTGGt                                                                                            <  2:20301/143‑1 (MQ=255)
                                                              ccgaGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGt                                                                                            >  1:20301/4‑146 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcct                                                                                        >  1:693574/1‑150 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcct                                                                                        >  1:248717/1‑150 (MQ=255)
                                                               cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                            >  1:734611/1‑113 (MQ=255)
                                                               cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                            <  2:734611/113‑1 (MQ=255)
                                                               cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctc                                                                                       >  1:399571/1‑150 (MQ=255)
                                                               cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctc                                                                                       <  2:997692/150‑1 (MQ=255)
                                                                aaGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCtt                                                                                                                                 <  2:1246900/107‑1 (MQ=255)
                                                                aaGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCtt                                                                                                                                 >  1:1246900/1‑107 (MQ=255)
                                                                aaGCCGCTGAGCATCAAGGACGTGATCTGGACACCACCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctcc                                                                                      <  1:1366053/150‑1 (MQ=255)
                                                                          aGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                             <  2:719773/150‑1 (MQ=255)
                                                                            cATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                              >  1:790200/1‑131 (MQ=255)
                                                                            cATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                              <  2:790200/131‑1 (MQ=255)
                                                                              tCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGa                                                                         <  1:235073/150‑1 (MQ=255)
                                                                              tCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGa                                                                         >  1:128760/1‑150 (MQ=255)
                                                                                aaGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                               >  1:490018/1‑94 (MQ=255)
                                                                                aaGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                               <  2:490018/94‑1 (MQ=255)
                                                                                aaGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACg                                                                       <  1:1040333/150‑1 (MQ=255)
                                                                                      gTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGt                                                                 >  2:1152696/1‑150 (MQ=255)
                                                                                          tCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTATCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCt                                                             <  1:1152696/150‑1 (MQ=255)
                                                                                            tGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTgg                                                           <  2:931947/150‑1 (MQ=255)
                                                                                             ggACACCCCCTGCCAAAGTCGCTAGCGCGCCTCGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTggg                                                          <  2:88984/150‑1 (MQ=255)
                                                                                             gatctcccccTGCCAAAGTCGCTAGCGAGCCTTGATTGCCTGCCAGTGCCGGGGCCCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTggg                                                          <  2:784900/145‑1 (MQ=255)
                                                                                                   ccccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCGGCCCCACCCACGGTGACGGCAAGTACCTGGGCagatc                                                    >  1:784900/1‑146 (MQ=255)
                                                                                                   ccccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTc                                                    >  1:369293/1‑150 (MQ=255)
                                                                                                           aaaGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGc                                                         >  2:103797/1‑137 (MQ=255)
                                                                                                           aaaGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGc                                                         <  1:103797/137‑1 (MQ=255)
                                                                                                              gTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCAc                                         <  2:693574/150‑1 (MQ=255)
                                                                                                              gTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCCCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCAc                                         <  2:875921/150‑1 (MQ=255)
                                                                                                                  cTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTggg                                                          <  2:1211185/129‑1 (MQ=255)
                                                                                                                  cTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTggg                                                          >  1:1211185/1‑129 (MQ=255)
                                                                                                                   tAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                             >  2:1043073/1‑109 (MQ=255)
                                                                                                                   tAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                             <  1:1043073/109‑1 (MQ=255)
                                                                                                                                ttgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                          <  1:2136/147‑1 (MQ=255)
                                                                                                                                ttgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                          >  2:2136/1‑147 (MQ=255)
                                                                                                                                 tgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTa                      >  2:1277748/1‑150 (MQ=255)
                                                                                                                                   aCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTAcc                    <  1:528819/150‑1 (MQ=255)
                                                                                                                                        tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGt                                                                           >  1:151726/5‑90 (MQ=255)
                                                                                                                                        tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGt                                                                           <  2:151726/86‑1 (MQ=255)
                                                                                                                                        cAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgcag               <  2:1396665/150‑1 (MQ=255)
                                                                                                                                                     ccgccCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAgg  >  2:390817/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAGG  >  NC_002947/292459‑292756

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: