Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,016 T→A 9.4% Q205L (CAG→CTG)  dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,0160TA 9.4% 174.4 / 8.2 72Q205L (CAG→CTG) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (5/58);  new base A (5/2);  total (12/60)
Fisher's exact test for biased strand distribution p-value = 3.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAA  >  NC_000913/2476868‑2477164
                                                                                                                                                    |                                                                                                                                                    
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                                                                                            <  1:197581/150‑1 (MQ=255)
                                                         aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATg                                                                                            <  2:714418/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)
                                                                                                                                                   ctgtctgtCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaaa  <  1:1553756/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAA  >  NC_000913/2476868‑2477164

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: