Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,016 T→A 9.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 9.3% 227.5 / 6.6 91Q205L (CAG→CTG) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (17/62);  new base A (8/2);  total (27/64)
Fisher's exact test for biased strand distribution p-value = 4.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAG  >  NC_000913/2476868‑2477165
                                                                                                                                                    |                                                                                                                                                     
tcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctg                                                                                                                                                      <  2:1121756/150‑1 (MQ=255)
  ccATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCATTTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtc                                                                                                                                                    <  2:1354408/150‑1 (MQ=255)
    aTCGCCACCGCATTATGATTCATCGCGGCAATATCAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCCGCAAATTAACCGCATAATGTTGCGCCGATCTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtctg                                                                                                                                                  <  1:1925612/150‑1 (MQ=255)
      cGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTATCCGCATAATGTTGCGCTCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCtgtctgtc                                                                                                                                                <  2:1867545/150‑1 (MQ=255)
       gCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTACAGCAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                               >  2:1542292/1‑150 (MQ=255)
       accACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGCTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCcggtctggaa                                                                                                                                               >  1:26199/2‑147 (MQ=255)
          aCCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGCACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATgg                                                                                                                                            <  2:1118773/150‑1 (MQ=255)
          aCCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATgg                                                                                                                                            <  1:1057338/150‑1 (MQ=255)
           ccccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                          <  2:717670/150‑1 (MQ=255)
            cacAATATGAATAATCGCGTCAATAACAGCTAAATCAGAACGATCAAAGCCATTTACATAATATGTCGGCAAATTAACCGCATAATGTTTCGCACAACTATGCCATTCATCCGTACCCTTGTCTTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                          <  2:302918/148‑1 (MQ=255)
             ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCATTGAGTCATGGAGc                                                                                                                                         >  2:390932/1‑150 (MQ=255)
             ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGc                                                                                                                                         <  1:1550434/150‑1 (MQ=255)
                atatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAAc                                                                                                                                      <  2:1612994/150‑1 (MQ=255)
                    gATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTAc                                                                                                                                  <  2:589614/150‑1 (MQ=255)
                    gATTCATCGCCGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTAc                                                                                                                                  <  1:70162/150‑1 (MQ=255)
                     aTTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTACATGGCTGTGTGTAATGGAGCAAcattca                                                                                                                                 >  1:1480874/1‑145 (MQ=255)
                       tCATCGCGGCAATAACAGCTAAATCAGACCGATCAAAGCCAATTCCAGCAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                               <  1:1183495/150‑1 (MQ=255)
                       tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAATCCACTTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                               <  1:1191301/150‑1 (MQ=255)
                       tCATCGAGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                               <  2:534926/150‑1 (MQ=255)
                         aTCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTGCCAGGCGGTGTGTGATGGAGCCACGGACACTGa                                                                                                                             >  2:182018/1‑150 (MQ=255)
                          tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGac                                                                                                                            >  2:360543/1‑150 (MQ=255)
                          tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCAGCATGCCTGTCTGTGATGGAGCAACGTACAggggc                                                                                                                            >  1:1286332/1‑146 (MQ=255)
                            gcgGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCAACAAGGCAGGGTGTCATGGGGCCAACTACAGTGacaa                                                                                                                          >  1:1177897/1‑149 (MQ=255)
                              ggCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTCCAAGCCTGTCTGTCATTGAGCAACGTACAGTGACACAg                                                                                                                        >  2:246159/1‑150 (MQ=255)
                                cAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGt                                                                                                                      <  2:2032894/150‑1 (MQ=255)
                                        gCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTAcc                                                                                                              <  1:400220/150‑1 (MQ=255)
                                               cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTgg                                                                                                       <  2:765022/150‑1 (MQ=255)
                                                aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTATCCATGCCTGTATGTCAAGGAGCCACGGTCAGTGACACCGGGTTATTACCGTGTTTGt                                                                                                      >  1:28610/1‑150 (MQ=255)
                                                aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGt                                                                                                      <  2:1908569/150‑1 (MQ=255)
                                                aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCCAGTGGGTGATGGAGCCACGGTCAGGGCCCCAGGGTAATTAACGTGTTGGt                                                                                                      >  1:1256778/1‑150 (MQ=255)
                                                    cGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAAGCCAGTGGGTCATCGAGCAACCTACCAGGAGACAGGGTAATTACCGGGTTGGGTaaa                                                                                                  >  2:398030/1‑150 (MQ=255)
                                                      aTCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGc                                                                                                <  1:774168/150‑1 (MQ=255)
                                                       tCAAAGCCCATTCCAGAAGATGTCGGCAAATTAACCGCATACTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                               <  2:1194536/150‑1 (MQ=255)
                                                       tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                               <  1:242369/150‑1 (MQ=255)
                                                       tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                               <  2:771771/150‑1 (MQ=255)
                                                         aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCGGTCTGTCCTGGTGCCACGTACAGTGACACCGGTTTATTACCGTATTGGTTATAGCATg                                                                                             >  2:331289/1‑150 (MQ=255)
                                                                aTTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGa                                                                                      <  2:1476798/150‑1 (MQ=255)
                                                                   ccAGAAGATTTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                   <  2:1004722/150‑1 (MQ=255)
                                                                   ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                   <  2:185375/150‑1 (MQ=255)
                                                                   ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                   <  1:705837/150‑1 (MQ=255)
                                                                   acAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                   <  2:1499997/149‑1 (MQ=255)
                                                                      gaagaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTAt                                                                                <  2:1215380/150‑1 (MQ=255)
                                                                         gaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCCTGGCTGTCTGTCATGGGGCAACGTACAATGACCCAGGTTTATTTACGTGTTGGTTTATTCATGTGTTTTAGCGTAttgc                                                                             >  1:24631/1‑148 (MQ=255)
                                                                               ggCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTa                                                                       <  1:840662/150‑1 (MQ=255)
                                                                                gCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCACCATGCCTGTCTGTCATGGAGCAACGTACACTGACACCGGTTAATTAACGTGTTGGTTAAAGCATTTGTTTGAGCGTATTCCGGGCTaa                                                                      >  2:1489369/1‑149 (MQ=255)
                                                                                 cAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                                     <  2:1246568/150‑1 (MQ=255)
                                                                                     ttAACCGCATAATGTTGCGCCCATCTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCCTGTCTGTCATGGAGCAACGTACAGTGACACATGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACt                                                                 <  2:1267953/150‑1 (MQ=255)
                                                                                     ttAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACACACt                                                                 <  1:1348888/150‑1 (MQ=255)
                                                                                      tAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCGATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTg                                                                <  1:1473562/150‑1 (MQ=255)
                                                                                       aaCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTgg                                                               <  2:1443508/150‑1 (MQ=255)
                                                                                        aCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGc                                                              <  2:824499/150‑1 (MQ=255)
                                                                                           gCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCAGTCTGTCATGGAGCAAAGGACACTGACACAGGTTTATTACCGTGGTTGTTAAATCATGTCTTTGAGCCTATTCCCGGCTTAATACTGGAAgt                                                           >  1:363720/1‑149 (MQ=255)
                                                                                            cATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGa                                                          <  1:1005297/150‑1 (MQ=255)
                                                                                             aTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAGGCCTGTCTGTCATGGAGCAACGTACAGTGACACATGTTAATTACCGTGTTGGTTAAAGCATGTATTTGAGAGTATTACGGGGAACATACTGTCAGGAt                                                         >  2:344834/1‑150 (MQ=255)
                                                                                             aTAATGTTGCGACCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGGCAAGGAGCACCGAACAGTGACACAAGGTAATTACCGTGTTGGTTAAAGCATGTCTTTGGCCGAATTGCCGGCTACATACATGCAGGAt                                                         >  2:109044/1‑150 (MQ=255)
                                                                                              tAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAtt                                                        <  2:2124730/150‑1 (MQ=255)
                                                                                              tAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAtt                                                        <  1:1654697/150‑1 (MQ=255)
                                                                                               aaTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAttt                                                       <  1:798679/150‑1 (MQ=255)
                                                                                               aaTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGGTTAAATACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAttt                                                       <  2:185845/150‑1 (MQ=255)
                                                                                                 tGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                     <  1:2067779/150‑1 (MQ=255)
                                                                                                   ttGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAAGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttcttc                                                   <  2:1901897/150‑1 (MQ=255)
                                                                                                    tGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTca                                                  <  1:656871/150‑1 (MQ=255)
                                                                                                       gcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcc                                               <  1:953252/150‑1 (MQ=255)
                                                                                                        cccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcca                                              <  2:599461/150‑1 (MQ=255)
                                                                                                          cAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTTACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                            <  1:1461870/150‑1 (MQ=255)
                                                                                                          cAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                            <  1:782428/150‑1 (MQ=255)
                                                                                                           aaCTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCa                                           <  1:360543/150‑1 (MQ=255)
                                                                                                                  ccATTCATCCGTACCTGAGTAGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGGTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTTCATACTGGCAGGATTTCTTCATCCATCAGTAAGtg                                    <  2:653706/150‑1 (MQ=255)
                                                                                                                   cATTCATCCGTACCGGAGTCGTGGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACTCAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtga                                   <  2:643991/150‑1 (MQ=255)
                                                                                                                        aTCCGTCCCGGCGTCGGTGCTCCAAGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGGTAATTACCGTGTTGGTTAAAGCATCTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagt                              <  2:173562/150‑1 (MQ=255)
                                                                                                                         tCCGTACCGGAGTCGTTTCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                             <  2:1534717/150‑1 (MQ=255)
                                                                                                                         tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                             <  2:363720/150‑1 (MQ=255)
                                                                                                                          ccGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTTCAGTGACACAGGTTAATTTCCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtga                            <  1:621453/150‑1 (MQ=255)
                                                                                                                          ccGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtga                            <  1:811487/150‑1 (MQ=255)
                                                                                                                              aCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTg                        <  1:729593/150‑1 (MQ=255)
                                                                                                                              aCCGGAGTCGTAGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTg                        <  1:390932/150‑1 (MQ=255)
                                                                                                                                cGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGAAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcg                      <  2:1579202/150‑1 (MQ=255)
                                                                                                                                  gAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAGTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCt                    <  2:1488382/150‑1 (MQ=255)
                                                                                                                                   aGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTg                   <  1:1264069/150‑1 (MQ=255)
                                                                                                                                    gTCGTTGCTCCATGCCTGTCTTGCATGGAGCAACGTACAGAGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGa                  <  1:873356/150‑1 (MQ=255)
                                                                                                                                       gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg               >  2:590660/1‑150 (MQ=255)
                                                                                                                                        ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc              <  2:158214/150‑1 (MQ=255)
                                                                                                                                        ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc              <  1:1221937/150‑1 (MQ=255)
                                                                                                                                          gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc            >  1:1606470/1‑150 (MQ=255)
                                                                                                                                          gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc            >  1:1034145/1‑150 (MQ=255)
                                                                                                                                           cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca           >  2:1287508/1‑150 (MQ=255)
                                                                                                                                           cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca           >  1:1571954/1‑150 (MQ=255)
                                                                                                                                           cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGATGGTTAAAGCATGTCTATGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca           >  1:1466640/1‑150 (MQ=255)
                                                                                                                                             ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc         <  1:761599/150‑1 (MQ=255)
                                                                                                                                              cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaact        <  1:1757985/150‑4 (MQ=255)
                                                                                                                                               aTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactc       <  1:901383/150‑5 (MQ=255)
                                                                                                                                                tGCCTGTCTGTCCTGGAGCAACGTACAGTGACACAGGTTAATTACCGCGGTGGTTAAAGCATGTCTTTGAGCGGATTCCGGGCTACATACTGGCCGGATTTCTTCATCCAGCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc      <  1:1211691/150‑1 (MQ=255)
                                                                                                                                                tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc      >  1:671356/1‑150 (MQ=255)
                                                                                                                                                 gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca     >  2:1922748/1‑150 (MQ=255)
                                                                                                                                                  ccTGTCTGTCATGGAGCAACGTACAGTTACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCATTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa    >  2:1686105/1‑150 (MQ=255)
                                                                                                                                                  ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa    >  2:614135/1‑150 (MQ=255)
                                                                                                                                                   ctgtctgtCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaaa   <  2:855531/150‑1 (MQ=255)
                                                                                                                                                    tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg  <  2:1177897/150‑1 (MQ=255)
                                                                                                                                                    tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg  <  1:2085553/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAG  >  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: