Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,332,871 A→C 6.7% S134R (AGT→AGG)  basS ← sensory histidine kinase in two‑component regulatory system with BasR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,8710AC 6.7% 310.9 / 11.9 116S134R (AGT→AGG) basSsensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (57/51);  new base C (4/4);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

Only 100 of 116 total aligned reads displayed.

CAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTG  >  NC_000913/4332723‑4333017
                                                                                                                                                    |                                                                                                                                                  
cAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                                                                                                   <  2:1011957/150‑1 (MQ=255)
 aaGCCGGGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGGCAGGCGACt                                                                                                                                                  <  1:1266841/150‑1 (MQ=255)
  aGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTATCGCCAGCAGTTCCAGATGCAAACGCACCCCCGACAGTGGCGTTCGCAGTTCGTGCGCGACGGCAGCGGTAAACAACCTTACCGTATCCAGCGTGCTGGTCCGGCGAcag                                                                                                                                                 >  2:1651405/1‑148 (MQ=255)
    ccGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAc                                                                                                                                               <  2:1277091/150‑1 (MQ=255)
       tGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAg                                                                                                                                            <  2:1699318/150‑1 (MQ=255)
       ggCCACTCACGGTGCTACATCACTGTGATGCGTTTGCGCCAGCAGTGCCAGAGGCAAACGCACCCCCGCCAGGGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTGATCCAGCGTGCTGGTCAGGCGACTGACCAg                                                                                                                                            <  1:858751/149‑1 (MQ=255)
        gCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGAGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGa                                                                                                                                           >  1:1110212/1‑150 (MQ=255)
        gCCACTAACGGTGATACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgc                                                                                                                                           <  1:264649/150‑2 (MQ=255)
          cACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGCCTGACCAgcac                                                                                                                                         <  1:869556/150‑4 (MQ=255)
          cACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                                                                                                         >  1:1386757/1‑150 (MQ=255)
          cACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                                                                                                         <  2:636977/150‑1 (MQ=255)
           aCTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCg                                                                                                                                        <  2:1395099/150‑1 (MQ=255)
             tAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGtt                                                                                                                                      <  1:1323108/150‑1 (MQ=255)
              aaCGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                                                                                                                     >  2:853907/1‑150 (MQ=255)
                  gTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACg                                                                                                                                 <  1:34777/150‑1 (MQ=255)
                  gTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATCCAAACGCACCCCCGCCCGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGGAAACACCCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACg                                                                                                                                 <  1:733681/150‑1 (MQ=255)
                     cTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                                                                                                              >  2:1011700/1‑150 (MQ=255)
                      tACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                                                                                                             >  2:1276060/1‑150 (MQ=255)
                        cATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaaa                                                                                                                           >  1:62103/1‑150 (MQ=255)
                         aTCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                                                                                                          <  1:552032/150‑1 (MQ=255)
                           cAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacca                                                                                                                        >  1:153724/1‑150 (MQ=255)
                               gtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAAAAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGGACGTTTAACGCCGATAAAACCGc                                                                                                                    >  1:808720/1‑150 (MQ=255)
                                 gATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                                                                                                  >  2:289596/1‑150 (MQ=255)
                                 gATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTAGCAGTTCGTGCGAGACGTCCGCGGTAAACACCCTTTCGTTATCAAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCCCCGCtt                                                                                                                  >  1:1248929/1‑150 (MQ=255)
                                  aTGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTc                                                                                                                 >  1:81737/1‑150 (MQ=255)
                                      gTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGCCTGACCAGATCGTTTAACGCCGAAACCACcgattcgatt                                                                                                             >  2:761198/1‑150 (MQ=255)
                                       ttttCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt                                                                                                            >  2:182735/1‑150 (MQ=255)
                                        tttCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                                                                                                           <  2:1386757/150‑1 (MQ=255)
                                           cGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAg                                                                                                        >  2:1370071/1‑150 (MQ=255)
                                            cccAGCAGTTCCAGCTGCAAACGCACCCCCGCCAGTGGCGTTCGCAGGTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                                                                                                       <  1:1651405/149‑1 (MQ=255)
                                             ccAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                                                                                                      <  2:385966/150‑1 (MQ=255)
                                                  agTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcg                                                                                                 >  1:670215/1‑150 (MQ=255)
                                                   gTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                                                                                                >  1:151730/1‑150 (MQ=255)
                                                     tCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg                                                                                              <  2:211069/150‑1 (MQ=255)
                                                     tCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg                                                                                              >  2:42469/1‑150 (MQ=255)
                                                       cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAACCCACCGCTTCGATTTCGAGGGTGGCGCtgtg                                                                                            >  1:847621/1‑150 (MQ=255)
                                                       cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                                                                                            <  1:742191/150‑1 (MQ=255)
                                                           tGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGTGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAGTTCGAGGGTGGCGCTGTGAATg                                                                                        >  2:1259251/1‑150 (MQ=255)
                                                           tGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg                                                                                        >  1:1712519/1‑150 (MQ=255)
                                                             cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCCCCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                                                                      >  1:464585/1‑150 (MQ=255)
                                                             cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                                                                      >  2:1026057/1‑150 (MQ=255)
                                                               aaCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACACCCTTTCGTTTTCCAGCGTGCTGGTCGGGCGACTGCCCAGAGCGTTTTACGCCGAACCCACCGCTTGGAGTGCGAGGGTGGCGGTGGGATTGGCaa                                                                                    >  1:1420723/1‑150 (MQ=255)
                                                                 cGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATg                                                                                  >  1:407734/1‑150 (MQ=255)
                                                                  gCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATgg                                                                                 <  1:1662418/150‑1 (MQ=255)
                                                                   cACCCCCGCCAGTGGCGTTCGCAGTTCGTGCTCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATggg                                                                                <  1:267924/150‑1 (MQ=255)
                                                                   cACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATggg                                                                                <  1:1698592/150‑1 (MQ=255)
                                                                   cACCCCCGACAGTGGTGTTCGCAGTTCGTGAGTGTCGTCAGCGGTAAAATACCTTTCGTTATCCATCGTGCTGGTCAGGCGACTGACCAGATCGTTTATCGCCGAAAACACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATggg                                                                                <  2:859174/150‑1 (MQ=255)
                                                                     cccccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGCTCGTTTAACGCCGAACCCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCg                                                                              >  1:533709/1‑150 (MQ=255)
                                                                     cccccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGAAAACAACCTTACGTGATCCAGCGTGCTGGTCAGGCGACTGACCAGAACGTTTACCGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGCCg                                                                              >  1:121019/1‑150 (MQ=255)
                                                                      ccccGGCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCGGTCGGTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCGTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGt                                                                             >  1:40986/1‑150 (MQ=255)
                                                                       cccGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGtt                                                                            <  2:81737/150‑1 (MQ=255)
                                                                          gCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTAGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAg                                                                         >  1:949662/1‑150 (MQ=255)
                                                                             aGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTg                                                                      >  2:317987/1‑150 (MQ=255)
                                                                               tGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTc                                                                    <  1:1370071/150‑1 (MQ=255)
                                                                                      cGCAGTTCGTTCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTg                                                             <  1:853907/150‑1 (MQ=255)
                                                                                      cGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTg                                                             >  2:450451/1‑150 (MQ=255)
                                                                                      cGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAGCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGCGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTg                                                             >  1:1565704/1‑150 (MQ=255)
                                                                                        cAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcg                                                           <  1:578346/150‑1 (MQ=255)
                                                                                               tGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTcc                                                    >  2:1206657/1‑150 (MQ=255)
                                                                                                gcgcGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCa                                                   <  2:153724/150‑1 (MQ=255)
                                                                                                  gcgACGTCAGCTGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGc                                                 <  1:445466/150‑1 (MQ=255)
                                                                                                  gcgACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGc                                                 <  2:62103/150‑1 (MQ=255)
                                                                                                  gcgACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGc                                                 >  2:1019303/1‑150 (MQ=255)
                                                                                                   cgACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGct                                                >  1:1462986/1‑150 (MQ=255)
                                                                                                    gACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGctc                                               >  1:1344391/1‑150 (MQ=255)
                                                                                                      cGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCtt                                             <  2:377885/150‑1 (MQ=255)
                                                                                                      cGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCtt                                             <  1:584034/150‑1 (MQ=255)
                                                                                                       gTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCttt                                            >  2:1279846/1‑150 (MQ=255)
                                                                                                           gCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGc                                        <  2:307276/150‑1 (MQ=255)
                                                                                                            cGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCa                                       <  1:1276060/150‑1 (MQ=255)
                                                                                                            cGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGGGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCa                                       <  1:965175/150‑1 (MQ=255)
                                                                                                               tAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCt                                    >  2:1681439/1‑150 (MQ=255)
                                                                                                                aaaCAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTc                                   >  1:1708782/1‑150 (MQ=255)
                                                                                                                 aaaaacCTTTCGTTATCCAGCGTCCAGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTACGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTAAAGATGTCGGCGGAGCGCGCTTCCAGCTCTTATTGCAGCTcc                                  >  1:23606/4‑150 (MQ=255)
                                                                                                                   caacCTTTCGTTATCCAGCGTACTGGTCAGGCGACTAACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGc                                <  2:1354400/150‑1 (MQ=255)
                                                                                                                      ccTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATATCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAg                             >  2:1483032/1‑150 (MQ=255)
                                                                                                                       cTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAGCAGATCGTTTAACGCCGATTCCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgc                            <  1:52131/150‑1 (MQ=255)
                                                                                                                        tttCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcg                           <  2:802311/150‑1 (MQ=255)
                                                                                                                           cGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcgggc                        >  2:1442402/1‑150 (MQ=255)
                                                                                                                           cGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcgggc                        >  2:1448894/1‑150 (MQ=255)
                                                                                                                           cGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcgggc                        >  2:1132353/1‑150 (MQ=255)
                                                                                                                             ttATCCAGCGTTCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGACACGACCGCTTCGATTTCGAGGGTGGCGCTGTGAATTGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGGAGATCGGCCAgcgggcgg                      <  2:1248929/150‑1 (MQ=255)
                                                                                                                             ttATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGTTGGAGCTGTGAATGGAAATGGGCGTTAAGTAGTCGGCGGTGCGAGCTTCCAGCTCTTTTTGCAGCTCCGCAAgcgggcgg                      >  2:292859/1‑150 (MQ=255)
                                                                                                                              tATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGcgggcggg                     >  2:884667/1‑150 (MQ=255)
                                                                                                                              tATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGcgggcggg                     <  1:729256/150‑1 (MQ=255)
                                                                                                                               aTCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGCTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGt                    <  1:103961/150‑1 (MQ=255)
                                                                                                                                 ccAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGa                  <  1:1026057/150‑1 (MQ=255)
                                                                                                                                     cgatcTGGTCAGTCGCCTGCCCAGATCGTTTAACGCCGAAACCACCGCTTCGATGTCGAGGGTGGCGCTGGGAATGGCAATGGGCGGTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATgcg              >  1:314364/5‑150 (MQ=255)
                                                                                                                                     cGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATgcg              >  2:1670282/1‑150 (MQ=255)
                                                                                                                                     cGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATgcg              >  1:1308325/1‑150 (MQ=255)
                                                                                                                                      gTGCTGGTCAGGCGACTGACCAGATCGTTTTACGCCGAAACCACCGCTTCGTTTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATgcgg             <  2:1110212/150‑1 (MQ=255)
                                                                                                                                       atcTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATgcggc            >  1:466999/3‑150 (MQ=255)
                                                                                                                                         cTGGTCAGTCGCCTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGt          >  2:1622851/1‑150 (MQ=255)
                                                                                                                                         cTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGt          <  2:151730/150‑1 (MQ=255)
                                                                                                                                           ggTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTAc        >  2:236578/1‑150 (MQ=255)
                                                                                                                                           ggTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTAc        >  2:1412962/1‑150 (MQ=255)
                                                                                                                                             tCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCg      <  2:1341054/150‑1 (MQ=255)
                                                                                                                                              cAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGc     <  1:884667/150‑1 (MQ=255)
                                                                                                                                                 gCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTg  <  2:464585/150‑1 (MQ=255)
                                                                                                                                                 gCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTg  <  2:468343/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTG  >  NC_000913/4332723‑4333017

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