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,9480GA93.0% 121.1 / ‑0.9 43V279V (GTC→GTTybeXCorC‑HlyC family protein YbeX
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (19/21);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCAATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCAT  >  NC_000913/690701‑691186
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
aGCCAATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGCATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGCGTCATCCGGGATTTTGAcc                                                                                                                                                                                                                                             >  2:182904/1‑250 (MQ=255)
         cGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAAAGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGGGAGTCATCCGGGATTTTAACATGAACCTg                                                                                                                                                                                                                                     >  1:16861/1‑250 (MQ=255)
              gACTGGAGGGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGTTGAGTCATCCGGGATTTTAACATGAACCTGaataat                                                                                                                                                                                                                               <  2:185924/251‑1 (MQ=255)
                                              cTGAAGCCCCATCAGCGACATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACt                                                                                                                                                                                               >  1:113318/1‑251 (MQ=255)
                                                aataccccAGCAGCGAAATAATCGCCGGAGACCAGACGTCGTAAGGAGAGACGTCCAGCGGTCCGCAGGCCACGAATAATAACGCCAGCAGTGAGCGAATGCGCTGGCGTTCAAGTAATGAGGCAAAAGCCATGTAGTTATTTATCCAGTGGCGGTTTTATTCATCCAGCTTCGGCTGGGTTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTTAACTgt                                                                                                                                                                                             <  2:331367/247‑2 (MQ=255)
                                                     cccAGCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAAcc                                                                                                                                                                                        <  1:78848/251‑1 (MQ=255)
                                                        aTCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCgtc                                                                                                                                                                                     >  1:95240/1‑251 (MQ=255)
                                                               aaaTAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCgtcgatgtcg                                                                                                                                                                              <  2:347216/251‑1 (MQ=255)
                                                                       gccgccagCCAGCCGTAGTACGGAGAGAAGGCCCGCGTTCCGCAGGCACCGACTAATAACGCCAGCCGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAACCCATGTAGTTATCTATCCCGTTTCGGGTTTATTCATCCAGCTTAGGCTGGGTGGAGTCATCCGGGATTTTAACATGAACCTGACTAATACGCCGACTGTCGGCCATCGCCCCTTTTAACTGGTAACCGTCGACGTCGATAGTTTc                                                                                                                                                                      >  2:221302/8‑251 (MQ=255)
                                                                         cgcAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGc                                                                                                                                                                    <  1:240548/251‑1 (MQ=255)
                                                                                 agacgcggtAAGGAGAGAAAGCCAGCGTTCCGCAGGAACCGAATAAAAACGCCAGCAGCAGGCGAAGTCGCTGGCGTTCAATTTATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCAACCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCGGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGcc                                                                                                                                                            <  2:16861/244‑1 (MQ=255)
                                                                                          aaGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATg                                                                                                                                                    <  1:325019/250‑1 (MQ=255)
                                                                                          aaGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGc                                                                                                                                                   <  1:25605/251‑1 (MQ=255)
                                                                                           aGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGcc                                                                                                                                                  <  1:10345/251‑1 (MQ=255)
                                                                                                    ggCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCt                                                                                                                                         <  1:90164/251‑1 (MQ=255)
                                                                                                     gccagcgtcccgacggaactgAATAATACCCCCAGCAGCAGGCAAACGCCGTGGCGTTGTATTAAGGAGGCAAAAGCCAGTGAGTTAGCTATCAATTTGCGATTTTATCCATCCAGCTTCGGCGGGGTTGGGTCATCCGGGATTTTAACATGAACCTGTATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGGAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAATATTCCCAAATGccgg                                                                                                                                        <  2:247482/231‑3 (MQ=255)
                                                                                                               cGCAGGCACCGAATAATAACGCCAGCAGCAGGCGACTGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCAcc                                                                                                                                 >  1:330750/1‑248 (MQ=255)
                                                                                                                         gAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAACTGCCTGCATCACCAGCCCACCgct                                                                                                                      >  1:153548/1‑247 (MQ=255)
                                                                                                                                                gAATGGGCGGGCGTTCAATTTATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTAAATAATACGCCGGCTGTCGGCCATCGCCACTTTGCACTAGTAACCGTCGATGTCGATAGTTTCGCCGCGCGCCTGAAGATGCCCAAATGCCTGCATCACCAGAACACCGATAGTGTCGACTTCTTCATCGCTAAAg                                                                                             <  2:113318/251‑1 (MQ=255)
                                                                                                                                                gAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAg                                                                                             <  1:437315/251‑1 (MQ=255)
                                                                                                                                                    gcgcTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTggg                                                                                         >  2:150502/1‑251 (MQ=255)
                                                                                                                                                             ttCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTACATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGaccg                                                                                >  2:186856/1‑248 (MQ=255)
                                                                                                                                                                          gCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTTGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTc                                                                   >  1:200504/1‑251 (MQ=255)
                                                                                                                                                                            aaaaGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCtt                                                                 <  1:182904/251‑1 (MQ=255)
                                                                                                                                                                                        aGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGcc                                                      >  1:439146/1‑251 (MQ=255)
                                                                                                                                                                                           tatctatcCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAg                                                    <  1:419367/250‑1 (MQ=255)
                                                                                                                                                                                                         cGGTTTTATTCATCCAGCTGCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGACCGCTTCGTTGAAGTCTTCAATGGAAGCCATTGCGCGTACGGTCCa                                     <  1:391934/251‑1 (MQ=255)
                                                                                                                                                                                                          ggTTTTATTCATCCAGCTTCGGCTTGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAg                                    <  1:442995/251‑1 (MQ=255)
                                                                                                                                                                                                          ggTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTcc                                      >  1:95808/1‑249 (MQ=255)
                                                                                                                                                                                                          ggTTTTATTCACCCATCTTCGTCTGGGGTGAGTCATCCGCGATTTGAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACGGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCATACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTGCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAg                                    <  1:264455/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTCCAGATTCGGCTGGGGTGAGACATCCGGGATTTTAACATGACCCTGACTAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAAACGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCACACCACCGATAGTGTCGACTTCTTCATCGCTAAGGTGGGTGCCGAACGCTTCGTTGAAGACTTCACTGGATGCCAGTGCGGGTAAGGTCCAGGCATgacgcc                          >  2:196254/1‑249 (MQ=255)
                                                                                                                                                                                                                      ccAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAACGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTCCGGTCCAGGTATGACGACTc                        >  2:355631/1‑251 (MQ=255)
                                                                                                                                                                                                                          cTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGCCTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTc                        >  1:31843/1‑247 (MQ=255)
                                                                                                                                                                                                                           ttCGGCTGGGGAGAGTCATCCGGGAATTTAACATGAAACTGCATAATACGCCGACTGTCGGCCCTCGCCACTTTGAACTGGTCACAGTCGATGTCGATAGTTTCGCCACACGCCGGAAGATGCCCAACTTCCTGCATCACCAGCCCACCGATAGTGTCGACTTCTTTATCGCTCAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGACAGTGCGGGTACGGTCCATGTCTGAAGACTCAGCTg                   >  2:431027/1‑251 (MQ=255)
                                                                                                                                                                                                                           gtCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTCGAACTGGTAACCGCCGATGTCGATAGTTTCGCCACTCGCCGGAAGATGCCCAAATGCCGTCATCACGAGACCACCGATAGTGTCGACTTCTGCATCGCATAAGGGGGTGCCGAACGCATCGTGGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGTCGACTCAGCTg                   <  1:431027/250‑1 (MQ=255)
                                                                                                                                                                                                                              ggCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGACGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCAAGTGCGCGTACGGTCCAGGAATGACGACTCAGCTGACg                >  1:89217/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggggTTAGTCATCCGGGATTTTAACATGAACCTGATTAATACTCCGACTGTCGGCCCTCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAATATGCCCAAATGCCTGCATCACCAGGCCACCGATAGTGTCGACTTCTTCATCTCTAAAGTGGTTGCCGAACGCTTCGTTGAAGTCTTCAATGTAAGCCAGTGCGCGTCCGGTCCAGGTATGACGACTCAGCTGACGGAAg            <  1:418590/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggggTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGa              >  1:95824/1‑249 (MQ=255)
                                                                                                                                                                                                                                  ggggTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGa              >  1:357167/1‑249 (MQ=255)
                                                                                                                                                                                                                                   gggTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGt           >  2:69428/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCAGGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGata      <  1:332333/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGGGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGat       >  2:126618/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aTCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATatcat  <  2:230561/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
AGCCAATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCAT  >  NC_000913/690701‑691186

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