Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,022 T→G 7.0% D203A (GAC→GCC)  dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,0220TG 7.0% 232.5 / 7.1 91D203A (GAC→GCC) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (23/61);  new base G (4/3);  total (27/64)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGAT  >  NC_000913/2476874‑2477170
                                                                                                                                                    |                                                                                                                                                    
cGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTATCCGCATAATGTTGCGCTCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCtgtctgtc                                                                                                                                                     <  2:1867545/150‑1 (MQ=255)
 gCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTACAGCAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                                    >  2:1542292/1‑150 (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                                                                                        <  1:705837/150‑1 (MQ=255)
                                                             ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                        <  2:185375/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)
                                                                                                                                             ctcactgtcATGGAGAAACGTACAGTGACAAATGTTAATTACCCTGTTGGTTAAAGCATATCTTTGAGCGTATTCCTAGATACATACTGCCACGATTTCATAATACAACAGAAAGTGATGAGTGAGTTGCGCTGACCACGCAACATaaaa        >  2:944158/5‑150 (MQ=255)
                                                                                                                                              tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg       <  2:1177897/150‑1 (MQ=255)
                                                                                                                                              tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg       <  1:2085553/150‑1 (MQ=255)
                                                                                                                                               gtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGt      >  2:1520774/1‑150 (MQ=255)
                                                                                                                                               gtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGt      <  1:1933558/150‑1 (MQ=255)
                                                                                                                                                   gtcATGGAGCAACGTACAGTGTCACAGGTTAATTACCGTGTGGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGAt  <  1:388790/150‑1 (MQ=255)
                                                                                                                                                   gtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGAt  >  2:849898/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGAT  >  NC_000913/2476874‑2477170

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: