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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,0140CA 5.4% 243.7 / 4.7 94A206S (GCA→TCA) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (19/66);  new base A (3/2);  total (26/68)
Fisher's exact test for biased strand distribution p-value = 9.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

CTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAA  >  NC_000913/2476866‑2477163
                                                                                                                                                    |                                                                                                                                                     
cttcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGcc                                                                                                                                                      <  2:1243678/150‑1 (MQ=255)
  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                                                                                             <  2:771771/150‑1 (MQ=255)
                                                         tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                             <  1:242369/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                                                      <  1:1654697/150‑1 (MQ=255)
                                                                                                tAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAtt                                                      <  2:2124730/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            <  1:1221937/150‑1 (MQ=255)
                                                                                                                                          ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc            <  2:158214/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)
                                                                                                                                                    |                                                                                                                                                     
CTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAA  >  NC_000913/2476866‑2477163

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