Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 690,948 G→A V279V (GTC→GTT)  ybeX ← CorC‑HlyC family protein YbeX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913690,9480GA97.9% 156.2 / ‑5.6 47V279V (GTC→GTT) ybeXCorC‑HlyC family protein YbeX
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/23);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCTTCTTC  >  NC_000913/690709‑691193
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
gCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCGGAAGCCCCAACAGCGAAATAAGCGCCGCATGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTg                                                                                                                                                                                                                                            <  1:403512/251‑1 (MQ=255)
         gggAGGGGACGGCGGTTAAAGGGCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACg                                                                                                                                                                                                                                   <  2:402294/250‑1 (MQ=255)
          ggAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACg                                                                                                                                                                                                                                   <  1:358933/250‑1 (MQ=255)
                         tAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCAt                                                                                                                                                                                                                    <  2:464704/250‑1 (MQ=255)
                           aaGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGc                                                                                                                                                                                                                 >  2:113800/1‑251 (MQ=255)
                                  aCGCCTGAAGGACCATCAGAGCACTAATCGCCGAAGGCCAGACAGCGACAGGGGAGAAGGCAAGCGTTCCGCAGGCACCGAATAATAACGCAAGCAGCAGAAGAATGCACTGGCGTTCAATTAATGAGGCAAACGCCATATAGTTATATATCCAGTTTCGGTTTTATTCATACAGCTTCGACTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATGGACACtttt                                                                                                                                                                                                          >  2:28436/1‑250 (MQ=255)
                                        gAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTgg                                                                                                                                                                                                    >  1:501816/1‑251 (MQ=255)
                                                     cGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCgtcgatgt                                                                                                                                                                                       <  1:324108/251‑1 (MQ=255)
                                                      gAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCAAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCgtcgatgtc                                                                                                                                                                                      >  1:61162/1‑251 (MQ=255)
                                                           aatCGCCGCAGGCCAGACGTCGTAAGGAGAGACGGCCAGCGTTACGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATCACTGAGGCAAAAGCCCTGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGATTAGACTGGGGTGAGTCATCCGGGATTTTAACATGCACCTGAATAATTCGCCGACTGTCGGCCATCGCCACTTTGACCTGGTGACCgtcgaggtcgag                                                                                                                                                                                   >  2:133709/1‑248 (MQ=255)
                                                            atCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGt                                                                                                                                                                                >  1:212544/1‑251 (MQ=255)
                                                                    aGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTAAATTAATGAGGCAAAAGCCATGTCGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGCTGTCGATAGTTTCGCCAc                                                                                                                                                                        >  2:21768/1‑251 (MQ=255)
                                                                                      agagaAGGCCAGCGTTCCGCAGGCAACGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTGCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCaa                                                                                                                                                      <  2:344648/251‑1 (MQ=255)
                                                                                         gaAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATg                                                                                                                                                   <  1:264926/251‑1 (MQ=255)
                                                                                           aGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGcc                                                                                                                                                 <  2:413768/251‑1 (MQ=255)
                                                                                             gCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTg                                                                                                                                               <  1:5915/251‑1 (MQ=255)
                                                                                               cAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCa                                                                                                                                             <  1:266910/251‑1 (MQ=255)
                                                                                                 gcggtCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCAGCCAGCGCCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATc                                                                                                                                           <  2:187689/247‑1 (MQ=255)
                                                                                                     tCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCAcc                                                                                                                                        <  2:567054/250‑1 (MQ=255)
                                                                                                      ccGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCTCCa                                                                                                                                       <  1:146177/250‑1 (MQ=255)
                                                                                                              aCCGAATAATAACGCCAGCAGCAGGCGAATGCGATGGCGTTCAATTAATGAGGCAAAAGCCATGGAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGATGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGGTGTCGATAGTTTTGCCACGCTCCGGAAGATGCCCAATTGCCTTCATAACCAGACAACCgc                                                                                                                              >  2:245897/1‑250 (MQ=255)
                                                                                                                cGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCCCGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGAt                                                                                                                             >  2:49493/1‑250 (MQ=255)
                                                                                                                cGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGAt                                                                                                                             >  1:3127/1‑250 (MQ=255)
                                                                                                                           gCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGCcttc                                                                                                                 >  2:117136/1‑251 (MQ=255)
                                                                                                                                   cagGAGAATGAGCTGTCGTTCAATTAATGAAGCAAAAGCCATGTAGTTATCTATCCATGTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGGCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCgt                                                                                                         <  2:452242/251‑2 (MQ=255)
                                                                                                                                          aTGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCCCCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                                                                  >  2:115081/1‑251 (MQ=255)
                                                                                                                                                   cGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGaa                                                                                         <  1:28252/251‑1 (MQ=255)
                                                                                                                                                    gTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAAc                                                                                        >  2:506240/1‑251 (MQ=255)
                                                                                                                                                        aaTTAATGAGGTAAAAGCCATGTAGTTAACTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTTGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGGCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCtt                                                                                    <  1:416246/251‑1 (MQ=255)
                                                                                                                                                                     aaaGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                                                                       >  2:530635/1‑251 (MQ=255)
                                                                                                                                                                      aaGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTTAAGTCTTc                                                                       >  1:526119/1‑250 (MQ=255)
                                                                                                                                                                      aaGCCATGTAGTTATCTATCCAGTTTCGGTTTTATGCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCa                                                                      <  2:571523/251‑1 (MQ=255)
                                                                                                                                                                           aTGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATgg                                                                  >  1:302090/1‑250 (MQ=255)
                                                                                                                                                                                                cGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGGCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCa                                            <  2:302090/251‑1 (MQ=255)
                                                                                                                                                                                                                       cTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCCATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGCCgg                      >  2:299950/1‑250 (MQ=255)
                                                                                                                                                                                                                        tGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCCGTGCGCGTACGGTCCAGGTCTGACGACTCCGCTGACGGa                     >  2:308500/1‑250 (MQ=255)
                                                                                                                                                                                                                         ggggTGATTCATCCGGGAGTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTGGAACTGGTAACCGTCGATGTCGATAGTTTCGCCTCGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGaa                    <  2:417960/250‑1 (MQ=255)
                                                                                                                                                                                                                         ggggTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGGCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAg                   <  2:98979/251‑1 (MQ=255)
                                                                                                                                                                                                                          gggTGAGTCATCCGGGCTTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAACCAGTGCGCGTACGGTCCCGGTATGAAGACGCAgatgccgggagt                  >  2:251289/1‑240 (MQ=255)
                                                                                                                                                                                                                           ggTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTc                 <  2:108152/251‑1 (MQ=255)
                                                                                                                                                                                                                            gTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCg                >  1:459001/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTCATCAGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGCCCACCGATAGTGTCGACTACTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTACATGGAAGCCAGTGCGCGTACGGTCACGGTATGACGACTCAGCGGACGGAAGCCGatat            >  2:425486/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aTCCGGGATTTTAACATGAACCTGAATAAAACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATatcat         <  2:501816/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGAAGGAAGTCGATATCAtctt      >  2:292186/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCttcttc  >  1:329656/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCttcttc  >  1:70597/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTTTTAACATGAACCTGAACAATCCGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCttcttc  <  2:546535/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
GCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCTTCTTC  >  NC_000913/690709‑691193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑