Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,016 T→A 7.3% 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 7.3% 190.3 / 3.6 73Q205L (CAG→CTG) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (13/52);  new base A (4/2);  total (19/54)
Fisher's exact test for biased strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAG  >  NC_000913/2476868‑2477165
                                                                                                                                                    |                                                                                                                                                      
tcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTACCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctg                                                                                                                                                       <  1:647031/150‑1 (MQ=255)
  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:1607435/1‑150 (MQ=255)
                                                                                                                                           cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca           <  1:519558/150‑1 (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)
                                                                                                                                                    |                                                                                                                                                      
TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTG‑ACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAG  >  NC_000913/2476868‑2477165

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