Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,018 T→A 5.6% R204S (AGA→AGT) ‡ dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,0180TA 5.6% 184.3 / 2.4 72R204S (AGA→AGT) ‡dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (14/50);  new base A (2/2);  total (18/54)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

CCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGT  >  NC_000913/2476870‑2477166
                                                                                                                                                    |                                                                                                                                                     
ccATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCAAAATCAGAACGATCAAAGCCAATTCCAGACGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtc                                                                                                                                                      <  2:314660/150‑1 (MQ=255)
  aTCGCCACCCCAATATGAATCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtctg                                                                                                                                                    <  2:494441/150‑1 (MQ=255)
     gCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAAGCCAgtcagtca                                                                                                                                                 >  1:41482/1‑150 (MQ=255)
          cccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                            <  1:441386/150‑1 (MQ=255)
           ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGc                                                                                                                                           <  1:1314858/150‑1 (MQ=255)
           ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGc                                                                                                                                           <  2:1807729/150‑1 (MQ=255)
              atatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAAc                                                                                                                                        <  1:1539764/150‑1 (MQ=255)
               tatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATAAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAACg                                                                                                                                       <  1:1734988/150‑1 (MQ=255)
                   aTTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGTCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACa                                                                                                                                   <  1:874952/150‑1 (MQ=255)
                     tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                                 <  1:1332348/150‑1 (MQ=255)
                       aTCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGAATAATGTTGCGCCCAACTATGCCATACATCCGTAACGGAGTCGTTGCTCAAAGCCAGTCTTTCATGGAGCAACCTTAAGTG‑a                                                                                                                              >  1:1098918/1‑150 (MQ=255)
                        tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ac                                                                                                                             <  1:438949/150‑1 (MQ=255)
                          gcgGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATTCCTGGCTGTGATGGAGCAACGTACAGGG‑accc                                                                                                                           >  1:1020906/1‑148 (MQ=255)
                          gcgGCAATAACAGCTAAATCAGAACGAGCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑acac                                                                                                                           <  1:1531089/150‑1 (MQ=255)
                             gCAATAACAGCTAAATCAGAACGATCAAAGCCAAATCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAgg                                                                                                                        <  1:1164987/150‑1 (MQ=255)
                               aaTAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGtt                                                                                                                      <  1:177729/150‑1 (MQ=255)
                                 tAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTaa                                                                                                                    <  1:1611719/150‑1 (MQ=255)
                                 tAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCGGTGTGTCATTGAGCAACGTACAGTG‑ACACAGGTTaa                                                                                                                    >  2:1051515/1‑150 (MQ=255)
                                    cAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTa                                                                                                                 <  1:828934/150‑1 (MQ=255)
                                      gCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTGCATGCCTGTCATGCATGGAGCAACTTACAGTG‑ACACAGGTTAATTaac                                                                                                               >  2:121040/1‑148 (MQ=255)
                                      gCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGGAACGTAAAGTG‑ACACAGGTTAATTaac                                                                                                               >  2:125712/1‑148 (MQ=255)
                                       cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCg                                                                                                              <  2:1055332/150‑1 (MQ=255)
                                         aaaTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCgtg                                                                                                            <  2:757284/150‑1 (MQ=255)
                                         aaaTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCCTGTCTGTCATGGAGCATCGTACAGTG‑ACACAGGTTAATTTCCgtt                                                                                                            >  2:83029/1‑149 (MQ=255)
                                            tCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTg                                                                                                         <  2:1278171/150‑1 (MQ=255)
                                             cAGAACGAGCAAAGGCAATTCCTGAAGAGGTCGGCTAATTAACCGCTTAATGTTGATCCCGACTATGCCATTGATCTGTACAGGAGTCGTGGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTgg                                                                                                        <  1:171748/150‑1 (MQ=255)
                                             cAGAACGAGCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTAGGCCATTCATCCGTACCGGAGTCGTCGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTgg                                                                                                        <  2:219540/150‑1 (MQ=255)
                                              aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGt                                                                                                       <  1:1430174/150‑1 (MQ=255)
                                                aaCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTa                                                                                                     <  1:585989/150‑1 (MQ=255)
                                                  cGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTGGCCAGGCCGGGATTTGAAAGAGCAACCTAAAGTGTAACCAGTTTAATAACCGGTTTGGTTaa                                                                                                    >  2:769958/1‑150 (MQ=255)
                                                   gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTCCCTGCCTGTCTGGGATGGGGCAACGTACAGTT‑ACACACGTTATTTTACGTGTTGGTTaaaa                                                                                                  >  2:1131985/1‑149 (MQ=255)
                                                   gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAg                                                                                                  <  1:1563624/150‑1 (MQ=255)
                                                    aTCACAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGc                                                                                                 <  2:1726606/150‑1 (MQ=255)
                                                    aTCAAAGCCAACTCCCGACGCTGTCGGCAACTTAACCGCATAATGTGGCGCCCAACTATGCCATTCATCCGTACCGGAGTCTTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGc                                                                                                 <  2:86337/150‑1 (MQ=255)
                                                      cAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCAt                                                                                               <  2:1361824/150‑1 (MQ=255)
                                                       aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATg                                                                                              <  1:1586326/150‑1 (MQ=255)
                                                       aaaGCAAATTCCAGAAGATGTCGGCAAATTATCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATg                                                                                              <  1:208581/150‑1 (MQ=255)
                                                          gCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAAGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCt                                                                                           >  1:1316233/1‑150 (MQ=255)
                                                           ccAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCtt                                                                                          <  2:1830657/150‑1 (MQ=255)
                                                               ttCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                      <  1:818045/150‑1 (MQ=255)
                                                                 ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATAAGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                    <  1:1812117/150‑1 (MQ=255)
                                                                     aagaTGTAGGCAAATTATCCGCATAATGTTGCGCGCAACAATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTAtt                                                                                <  2:807563/150‑1 (MQ=255)
                                                                         tGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTGATGGCGCCAAGGACAGTG‑ACACAGGTTTATTTCCGGGTTGGTTAAAGCATGTCGTTGAGCGTATTACgg                                                                            >  2:1318929/1‑150 (MQ=255)
                                                                         tGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATGCCgg                                                                            <  2:1361463/150‑1 (MQ=255)
                                                                          gTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                           <  2:518854/150‑1 (MQ=255)
                                                                           tCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGc                                                                          <  1:1054948/150‑1 (MQ=255)
                                                                               cAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                                      <  2:1350306/150‑1 (MQ=255)
                                                                                               tGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                      <  2:1267570/150‑1 (MQ=255)
                                                                                               tGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGGCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                      <  1:125712/150‑1 (MQ=255)
                                                                                               tGTTCCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTAAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                      <  1:121040/150‑1 (MQ=255)
                                                                                                gTTGCGCCCAACTATGCCAGTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATACCGGGCTACATACTGGCAGGATttctt                                                     <  2:1808256/150‑1 (MQ=255)
                                                                                                      cccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTGCGGGCTTCAAACTGGCAGGATttcttcctgca                                               >  1:50075/1‑145 (MQ=255)
                                                                                                               gCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGt                                      <  1:692281/150‑1 (MQ=255)
                                                                                                                ccATTCATCCGTCCCGGAGTCGTTGCTCCATGCCTGGCTGTCATTGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGGTGGTTAAAGCATGTCGTTGAGCGGATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtg                                     <  2:154644/150‑1 (MQ=255)
                                                                                                                   ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                                  <  1:698926/150‑1 (MQ=255)
                                                                                                                   ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                                  <  2:1572671/150‑1 (MQ=255)
                                                                                                                       tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                              <  1:83029/150‑1 (MQ=255)
                                                                                                                        ccGTACCTGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtga                             <  1:1871572/150‑1 (MQ=255)
                                                                                                                            aCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTg                         <  1:1538023/150‑1 (MQ=255)
                                                                                                                               ggAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcgc                      <  1:1763442/150‑1 (MQ=255)
                                                                                                                                gAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCt                     <  1:1834750/150‑1 (MQ=255)
                                                                                                                                 tgTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTg                    <  1:694110/149‑1 (MQ=255)
                                                                                                                                     gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                >  1:643300/1‑150 (MQ=255)
                                                                                                                                       tGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg              >  2:496833/1‑150 (MQ=255)
                                                                                                                                        gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGGTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc             <  1:532445/150‑1 (MQ=255)
                                                                                                                                        gCTCCATGACAGACAGGCATGGAGCAACGTACAGTG‑ACACAGGTAAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc             >  2:1120561/1‑150 (MQ=255)
                                                                                                                                         cTCCATGCCTGTCTGTCATGTAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCTTCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca            <  1:1500721/150‑1 (MQ=255)
                                                                                                                                         cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca            <  1:519558/150‑1 (MQ=255)
                                                                                                                                         cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca            >  1:1607435/1‑150 (MQ=255)
                                                                                                                                         cTCCATGACAGACAGGCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca            >  2:1821093/1‑150 (MQ=255)
                                                                                                                                          tCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcat           >  2:1505233/1‑150 (MQ=255)
                                                                                                                                              tGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc       <  1:1051515/150‑1 (MQ=255)
                                                                                                                                              tGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc       <  2:1515454/150‑1 (MQ=255)
                                                                                                                                               gCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca      <  2:496225/150‑1 (MQ=255)
                                                                                                                                                ccTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa     >  2:958312/1‑150 (MQ=255)
                                                                                                                                                 ctgtctgtCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaaa    >  1:1740790/1‑150 (MQ=255)
                                                                                                                                                  tgtctgtcATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg   >  1:912384/1‑150 (MQ=255)
                                                                                                                                                  tgtctgtcATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGGCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg   <  1:1605669/150‑1 (MQ=255)
                                                                                                                                                   gtctgtcATGGAGCAACGTACAGTG‑ACACAGGTTAATTATCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCAATGAGTAAGTGATTAGTGAGTTGCGCTGTCGGCGCATCATCAAAGt  <  2:670206/150‑1 (MQ=255)
                                                                                                                                                   gtctgtcATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGt  >  2:184675/1‑150 (MQ=255)
                                                                                                                                                   gtctgtcATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGt  <  2:780538/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGT  >  NC_000913/2476870‑2477166

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: ‑