Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,014 C→A 8.8% A206S (GCA→TCA)  dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,0140CA 8.8% 191.5 / 6.9 75A206S (GCA→TCA) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (7/60);  new base A (4/2);  total (12/63)
Fisher's exact test for biased strand distribution p-value = 3.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAA  >  NC_000913/2476866‑2477163
                                                                                                                                                    |                                                                                                                                                     
cttcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGcc                                                                                                                                                      <  1:700392/150‑1 (MQ=255)
  tcCCATCGCCACCCCAATATGATTGATCGCGGCAATAACATCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCTTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctg                                                                                                                                                    <  2:1331656/150‑1 (MQ=255)
     cATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCCTGCctgtat                                                                                                                                                 >  1:49601/1‑148 (MQ=255)
      aTCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtctg                                                                                                                                                <  2:588471/150‑1 (MQ=255)
      aTCGCCACCCCAATATGATGCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGAAGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtctg                                                                                                                                                <  2:709008/150‑1 (MQ=255)
         gCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                             <  2:998053/150‑1 (MQ=255)
         gCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                             <  2:213600/150‑1 (MQ=255)
          ccACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGCACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCTTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAt                                                                                                                                            <  2:882738/150‑1 (MQ=255)
           cACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATg                                                                                                                                           <  2:817848/150‑1 (MQ=255)
           cACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGAGCAAAGCCAATTGCTGTAGAGGTCGGCAAATTTACCGCATAATGTTGCGCCCAACTATGCCATTCATGCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATg                                                                                                                                           <  1:219677/150‑1 (MQ=255)
              cccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                        <  2:1456299/150‑1 (MQ=255)
               ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACAATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGc                                                                                                                                       <  1:557270/150‑1 (MQ=255)
                  atatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTACCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAAc                                                                                                                                    <  1:48646/150‑1 (MQ=255)
                   tatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACg                                                                                                                                   <  2:1383012/150‑1 (MQ=255)
                   tatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAACg                                                                                                                                   <  2:1345465/150‑1 (MQ=255)
                      gATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTAc                                                                                                                                <  1:159179/150‑1 (MQ=255)
                       aTTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCACAAGCTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACa                                                                                                                               <  1:1357517/150‑1 (MQ=255)
                            tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGACTGTCATGGAGCAACGTACAGTGac                                                                                                                          <  2:1732388/150‑1 (MQ=255)
                                 gCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTACATGCCAGTCTTTCCTGGAGCAACGGACACGGACAAAgg                                                                                                                     >  1:1356823/1‑150 (MQ=255)
                                  cAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGt                                                                                                                    <  1:1707381/150‑1 (MQ=255)
                                  cAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGt                                                                                                                    <  1:579988/150‑1 (MQ=255)
                                   aaTAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGtt                                                                                                                   <  2:752344/150‑1 (MQ=255)
                                      aaCAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAAt                                                                                                                <  1:1824267/150‑1 (MQ=255)
                                           ctaatacAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATACTGCTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCg                                                                                                           <  2:93242/144‑1 (MQ=255)
                                           cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCTCATAATGTTGCGCCCAACTATGCGATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTAGCg                                                                                                           <  1:1654432/150‑1 (MQ=255)
                                           cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAAACGCATAATGTTGCTCCCAACTATGCCATTCATCCGTACCGGATTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCg                                                                                                           <  1:900234/150‑1 (MQ=255)
                                            tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCgt                                                                                                          <  2:1648660/150‑1 (MQ=255)
                                               aTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTTACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGtt                                                                                                       <  2:1060960/150‑1 (MQ=255)
                                                 cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTgg                                                                                                     <  1:201790/150‑1 (MQ=255)
                                                    aaCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTCCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTa                                                                                                  <  2:1287329/150‑1 (MQ=255)
                                                    aaCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTa                                                                                                  <  1:1550279/150‑1 (MQ=255)
                                                        aTCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGc                                                                                              <  2:1504860/150‑1 (MQ=255)
                                                        aTCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGc                                                                                              <  2:1040746/150‑1 (MQ=255)
                                                         tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCGATGCCTGTCTGTCATGTAGCATCGTACAGTGACACAGGTTAATTAGCGTGTTGGTTAAAGCa                                                                                             <  1:1467461/150‑1 (MQ=255)
                                                         tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                             <  1:1896156/150‑1 (MQ=255)
                                                           aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATg                                                                                           <  2:714418/150‑1 (MQ=255)
                                                           aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATg                                                                                           <  1:197581/150‑1 (MQ=255)
                                                                cAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCGATTCATCCGTACCGGAGTCGTTTCTACCGGGCGGGCTGTACTGGAGCAACCTAAAGTGACCAAGGATAATTACCGTGTTGGTTAAAGCaattattt                                                                                      >  2:281167/1‑143 (MQ=255)
                                                                 aaTTCCAGAAGATGTCTGCAAATTAACCGCATAATGTTGCGCCCAAATATGCCATTCATCAGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTg                                                                                     <  1:98967/150‑1 (MQ=255)
                                                                           gaTGTCTGCTATTTAACCGTATAATGTTGCGCCGTACTATGCCATTCATCCGTACCGTAGTCTTTGCTCCATGTCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTATTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTcc                                                                           <  1:1096545/150‑1 (MQ=255)
                                                                           gaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTcc                                                                           <  1:1542370/150‑1 (MQ=255)
                                                                             tGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCgg                                                                         <  2:1313038/150‑1 (MQ=255)
                                                                                       ttAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACt                                                               <  2:854643/150‑1 (MQ=255)
                                                                                             gCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTTTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAgg                                                         <  1:1399223/150‑1 (MQ=255)
                                                                                                  aTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttc                                                    <  2:587903/150‑1 (MQ=255)
                                                                                                  aTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCGGTCATGGAGCAACGGACAGTGTCACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttc                                                    <  1:1491148/150‑1 (MQ=255)
                                                                                                   tGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                   <  1:1513853/150‑1 (MQ=255)
                                                                                                       gcgcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcat                                               <  1:1343835/150‑1 (MQ=255)
                                                                                                        cgcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatc                                              <  2:1356943/150‑1 (MQ=255)
                                                                                                            cAACTTTGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGTCACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                          <  2:703171/150‑1 (MQ=255)
                                                                                                             accTATGCCATGCATCCGTACCGGAGTCGGTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCa                                         <  2:904803/148‑1 (MQ=255)
                                                                                                                   gCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGt                                   <  2:1812813/150‑1 (MQ=255)
                                                                                                                    ccATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtg                                  >  2:1125099/1‑150 (MQ=255)
                                                                                                                    ccATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtg                                  <  2:851394/150‑1 (MQ=255)
                                                                                                                       ttCATCCGTACCGGAGTCGGTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCTGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                               <  2:726450/150‑1 (MQ=255)
                                                                                                                             cGTACCGGAGTCGTTGCTCCCTGCCGGTCTGTCATGGAGGAACGGACAGTGGCACAGGTTAATTTCCGTGTTTGGTAAAGCATGTGTTTGAGCGTAATTCGGGGTAAATACTGGCAGGCTTTCTTAATCCATCAGAAAGTGAggagtgag                         >  2:1330063/1‑150 (MQ=255)
                                                                                                                               tACCGGAGTCGTTGCTCCATGCCTGGCTGTGATGGGGCAACGTACAGTGACACAGGTTAATTTCCGTGGTGGTTAAAGCATGTCTTTGAGCGTATTACGGGGGACATACTGGCAGGATTTCTTCAAGCATCAGAAAGTGATGAGTGAGtt                       >  2:51470/1‑150 (MQ=255)
                                                                                                                                aCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTg                      <  1:1480470/150‑1 (MQ=255)
                                                                                                                                   ggAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcgc                   <  2:1084150/150‑1 (MQ=255)
                                                                                                                                    gAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCt                  <  1:872365/150‑1 (MQ=255)
                                                                                                                                        gggtGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGGTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACAAACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGGTGCGCTGACg              <  2:1565478/147‑1 (MQ=255)
                                                                                                                                        cggtGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg              <  2:1188944/147‑1 (MQ=255)
                                                                                                                                        cGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg              >  1:119967/1‑150 (MQ=255)
                                                                                                                                           tGCTCGATGCCTGTCTGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg           <  2:1403342/150‑1 (MQ=255)
                                                                                                                                           tGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg           <  2:407683/150‑1 (MQ=255)
                                                                                                                                           tGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg           >  1:59202/1‑150 (MQ=255)
                                                                                                                                           tGCTCAATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg           <  2:1481176/150‑1 (MQ=255)
                                                                                                                                            gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc          >  1:1471442/1‑150 (MQ=255)
                                                                                                                                            gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc          >  1:682014/1‑150 (MQ=255)
                                                                                                                                             cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca         >  1:1624470/1‑150 (MQ=255)
                                                                                                                                             cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca         <  1:1631271/150‑1 (MQ=255)
                                                                                                                                             cTCCATGCCTGTATGTCATGGAGCAACGTACAGTGACACAGGTAAATTACCGTGTTGGTAAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca         <  1:736555/150‑1 (MQ=255)
                                                                                                                                             cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca         >  1:841989/1‑150 (MQ=255)
                                                                                                                                               ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc       <  2:49601/150‑1 (MQ=255)
                                                                                                                                                   gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactcac   <  1:198108/150‑7 (MQ=255)
                                                                                                                                                   gCCTGTCTGTCATGGAGCAACGGACAGTGAGACAGGTTAATTACCGTGTTGGTTATAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca   <  2:45329/150‑1 (MQ=255)
                                                                                                                                                    ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa  >  1:1414943/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAA  >  NC_000913/2476866‑2477163

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