Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 292,607 C→G 5.4% 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 5.4% 299.6 / 6.0 93intergenic (+22/‑38)ssuA/ssuDaliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base G (3/2);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAG  >  NC_002947/292459‑292755
                                                                                                                                                    |                                                                                                                                                     
cccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgcc                                                                                                                                                      >  2:1109836/1‑150 (MQ=255)
cccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgcc                                                                                                                                                      >  2:709354/1‑150 (MQ=255)
cccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgcc                                                                                                                                                      >  2:709355/1‑150 (MQ=255)
 ccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgccg                                                                                                                                                     >  1:1365146/1‑150 (MQ=255)
    gagGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGccc                                                                                                                                                  <  1:545041/150‑1 (MQ=255)
       gtggtgGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCgg                                                                                                                                               <  1:1174717/150‑1 (MQ=255)
        tggtggCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGc                                                                                                                                              >  1:1043133/1‑150 (MQ=255)
            ggCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                          <  2:1878160/150‑1 (MQ=255)
             gCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcactgg                                                                                                                                         <  2:1172014/150‑6 (MQ=255)
             gCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCCGCTGCACCCTGCGGTTAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTagag                                                                                                                                         <  2:1484201/150‑1 (MQ=255)
                 cgcAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCAc                                                                                                                                     <  2:1601198/150‑1 (MQ=255)
                  gcAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCAcc                                                                                                                                    <  2:1385050/150‑1 (MQ=255)
                   cAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGagccagcc                                                                                                                                   <  1:1043553/150‑1 (MQ=255)
                   cAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCAccc                                                                                                                                   <  1:709355/150‑1 (MQ=255)
                   cAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCAccc                                                                                                                                   <  1:709354/150‑1 (MQ=255)
                       agaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                              <  1:1778874/150‑1 (MQ=255)
                       agaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCAttca                                                                                                                               <  1:1006188/150‑3 (MQ=255)
                        gaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATGTGGACACCCCCTGCCAAAGTCGCTCGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCGCCCTTgag                                                                                                                             >  1:357310/1‑150 (MQ=255)
                        gaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCAGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             <  1:1467127/150‑1 (MQ=255)
                          agaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgagga                                                                                                                           <  2:130711/150‑1 (MQ=255)
                           gaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgaggag                                                                                                                          >  2:1145257/1‑150 (MQ=255)
                            aTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             <  1:770359/146‑1 (MQ=255)
                            aTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             >  2:770359/1‑146 (MQ=255)
                               gCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAgacga                                                                                                                      <  2:607054/150‑1 (MQ=255)
                                  gACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTc                                                                                                                   <  1:1399590/150‑1 (MQ=255)
                                  gACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTc                                                                                                                   >  2:905366/1‑150 (MQ=255)
                                       cTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATg                                                                                                              >  2:250088/1‑150 (MQ=255)
                                        ttCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCCAGGGCGTGATCTGGACACCCCCTGCCAAAGTAGCTAGCGCGCCTTGATTCACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGCCGGCTCCAATgg                                                                                                             >  2:698222/1‑149 (MQ=255)
                                        ttCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGa                                                                                                             >  1:111724/1‑150 (MQ=255)
                                            aGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCt                                                                                                         >  2:1511388/1‑150 (MQ=255)
                                            aGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCt                                                                                                         <  1:1145257/150‑1 (MQ=255)
                                            aGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCt                                                                                                         >  2:426056/1‑150 (MQ=255)
                                                 cTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACa                                                                                                    <  1:426056/150‑1 (MQ=255)
                                                      gCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCTGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTc                                                                                               >  2:1036828/1‑150 (MQ=255)
                                                      gCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTc                                                                                               <  2:1248743/150‑1 (MQ=255)
                                                      gCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTc                                                                                               >  1:1248743/1‑150 (MQ=255)
                                                        tGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTg                                                                                             <  1:823677/150‑1 (MQ=255)
                                                         gATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTgg                                                                                            <  2:1378753/150‑1 (MQ=255)
                                                             cccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcc                                                                                        >  2:1135480/1‑150 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcactgg                                                                                                                                         >  2:1590003/1‑96 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcactgg                                                                                                                                         <  1:1590003/101‑6 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGCATGAAGGACTCCAATGAGCCTTAACATTTTCTGGTtcct                                                                                       >  1:98798/1‑150 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcct                                                                                       >  2:484733/1‑150 (MQ=255)
                                                              ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGATGGGTTTGAAATGACGACTCCAATGAGCCTTAACATTTTCTGGTtcct                                                                                       >  1:1867721/1‑150 (MQ=255)
                                                               cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctc                                                                                      <  2:1043133/150‑1 (MQ=255)
                                                                ccggcGCCGAGCATCAAGGACGTGACCTGGACACCCCCTACCAAAGGCGCTAGCGCGCCTTGATTGACCGCAAGTGCCGGGCCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGCGCCTTAACACTTTCTGGTtcctcc                                                                                     <  1:398197/146‑1 (MQ=255)
                                                                aaGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctcc                                                                                     <  1:484733/150‑1 (MQ=255)
                                                                aaGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAATCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcctcc                                                                                     <  1:373210/150‑1 (MQ=255)
                                                                 aGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTccc                                                                                    >  2:43135/1‑150 (MQ=255)
                                                                    cGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCcccac                                                                                 >  2:169498/1‑150 (MQ=255)
                                                                     gCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCcccacc                                                                                <  1:250088/150‑1 (MQ=255)
                                                                            aTCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTg                                                                         >  1:1611393/1‑150 (MQ=255)
                                                                                 ggACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATCGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGAGTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGc                                                                    >  2:977630/1‑150 (MQ=255)
                                                                                   tcGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                          >  2:1355502/2‑79 (MQ=255)
                                                                                   tcGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                          <  1:1355502/78‑1 (MQ=255)
                                                                                          cTGGACACCCCCTGCCAAAGTCGCTCGCGCGCCTTGATTGACTGCCCGTCCCGGGCCGCCGCCCCCGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTg                                                           <  1:698222/150‑1 (MQ=255)
                                                                                           tGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTgg                                                          <  1:899276/150‑1 (MQ=255)
                                                                                           tGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTgg                                                          <  1:8978/150‑1 (MQ=255)
                                                                                            ggACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTggg                                                         >  2:361357/1‑150 (MQ=255)
                                                                                                acCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                              >  1:928588/1‑77 (MQ=255)
                                                                                                acCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                              <  2:928588/77‑1 (MQ=255)
                                                                                                 ccctcTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                            >  1:937731/1‑126 (MQ=255)
                                                                                                 ccctcTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                            <  2:937731/122‑1 (MQ=255)
                                                                                                  ccccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTc                                                   >  1:67978/1‑150 (MQ=255)
                                                                                                   cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTg                                                                                                                                >  2:443092/1‑73 (MQ=255)
                                                                                                   cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTg                                                                                                                                <  1:443092/73‑1 (MQ=255)
                                                                                                   cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCCCATCTCGGGGCACTTcc                                                  >  2:1481769/1‑150 (MQ=255)
                                                                                                   cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTcc                                                  <  1:43135/150‑1 (MQ=255)
                                                                                                      tGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGaa                                               <  2:1365146/150‑1 (MQ=255)
                                                                                                        ccAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAgg                                             <  1:754074/150‑1 (MQ=255)
                                                                                                                gCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTg                                     >  1:1556823/1‑150 (MQ=255)
                                                                                                                 cTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACGTGGGCACTTCCGAAGGCGCACGTGc                                    >  1:1220707/1‑150 (MQ=255)
                                                                                                                      gcgcCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCg                               >  2:1879220/1‑150 (MQ=255)
                                                                                                                        gcCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAgg                                             >  2:1813533/1‑134 (MQ=255)
                                                                                                                        gcCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAgg                                             <  1:1813533/134‑1 (MQ=255)
                                                                                                                          cTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTg                                                           <  1:869900/118‑1 (MQ=255)
                                                                                                                          cTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTg                                                           >  2:869900/1‑118 (MQ=255)
                                                                                                                          cTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCa                           >  1:44655/1‑150 (MQ=255)
                                                                                                                          cTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACGTCCGAAGGCGCACGTGCGGTCGACCa                           >  1:1628763/1‑150 (MQ=255)
                                                                                                                            tgattgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACg                         <  2:729393/150‑1 (MQ=255)
                                                                                                                              attgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGa                              <  2:1199499/143‑1 (MQ=255)
                                                                                                                              attgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGa                              >  1:1199499/1‑143 (MQ=255)
                                                                                                                                   cTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCt                  <  1:1879220/150‑1 (MQ=255)
                                                                                                                                    tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGc                                    >  1:612910/1‑131 (MQ=255)
                                                                                                                                    tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGc                                    <  2:612910/131‑1 (MQ=255)
                                                                                                                                    tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCt                  <  2:197950/149‑1 (MQ=255)
                                                                                                                                    tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCt                  >  1:197950/1‑149 (MQ=255)
                                                                                                                                    tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTg                 <  2:30579/150‑1 (MQ=255)
                                                                                                                                      ccAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCTGTGCCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgca               >  1:1553649/1‑150 (MQ=255)
                                                                                                                                      ccAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACACCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgca               >  1:1787957/1‑150 (MQ=255)
                                                                                                                                           gCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGa          >  1:570827/1‑150 (MQ=255)
                                                                                                                                           gCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGa          >  2:1661382/1‑150 (MQ=255)
                                                                                                                                           gCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGa          <  1:1036828/150‑1 (MQ=255)
                                                                                                                                            ccGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCGTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGAt         >  1:814121/1‑150 (MQ=255)
                                                                                                                                            ccGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGAt         >  1:811171/1‑150 (MQ=255)
                                                                                                                                               gggCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCaa                                                                  <  1:1838178/90‑1 (MQ=255)
                                                                                                                                               gggCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCaa                                                                  >  2:1838178/1‑90 (MQ=255)
                                                                                                                                                  cgccgcCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCa   <  1:977630/150‑1 (MQ=255)
                                                                                                                                                   gccgccCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAg  >  2:1440780/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAG  >  NC_002947/292459‑292755

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: