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,9480GA98.0% 168.5 / ‑4.0 52V279V (GTC→GTT) ybeXCorC‑HlyC family protein YbeX
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/19);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

AATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCTTCTTC  >  NC_000913/690705‑691193
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
aaTAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGa                                                                                                                                                                                                                                                   >  1:131641/1‑250 (MQ=255)
         aGACTGGAGTGGACGGCGGTTAAAGGTCTACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATTGGGCAAAAGCCATGT‑AGTTAACTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGCTTTTAACATGAACCTgcatt                                                                                                                                                                                                                                          >  2:225736/1‑246 (MQ=255)
                                  gTCAACGCATGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCCGCCGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGCATAATACAGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGCGTCATCCGGGGTTTTAACATGCACCTGAATTATCAC‑CCGACTGTCGGCCATCGcccc                                                                                                                                                                                                               >  2:190816/1‑249 (MQ=255)
                                            gAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACt                                                                                                                                                                                                       <  1:317504/249‑1 (MQ=255)
                                                    aTCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGAGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCgtc                                                                                                                                                                                             >  1:111050/1‑251 (MQ=255)
                                                       aGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCgtcgat                                                                                                                                                                                          <  2:131641/251‑1 (MQ=255)
                                                                    ccgcAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCTCCGAATAATAACGCCAGCAGCAGTCTAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCAGCGCCACTTTGAACTGGTAACCGTCGATGTCTATAGTTTCg                                                                                                                                                                             <  2:362639/251‑1 (MQ=255)
                                                                                 gtcgtAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCCTGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGa                                                                                                                                                                 >  1:356683/1‑250 (MQ=255)
                                                                                                        ttCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTCACCGTCGATGTCGATAGTTTCGCCACGCACCGGAAGCTGCCCAAATGCCTGCA‑TCAcc                                                                                                                                        >  2:224351/1‑251 (MQ=255)
                                                                                                          ccGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAg                                                                                                                                      <  2:12170/251‑1 (MQ=255)
                                                                                                              aGGCACCGATTAATAACGCCAGCCGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGAGAGTCATCCGGGATTTTCACATGCACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTCCGCCAACCGCCGGAAGATGCCCAACTGCCTGCA‑TCCCCAgcccc                                                                                                                                  >  2:351756/1‑247 (MQ=255)
                                                                                                                               cgcagcagaagGCGAATGCGCTGGCGTCCAATTAATGAGGCAAAAGCCA‑GTGAGTTATCTATCCAGTTTCGGTTTTATCCATCCAGCTTCGGCTGGGGTGAGTAATCCGGGATTTTAACATGAACCGGAATAAT‑ACGCCGACTGTCGTCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTGCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCAActtc                                                                                                                 <  2:347484/249‑1 (MQ=255)
                                                                                                                                 cagcagcagGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGActtctt                                                                                                               <  1:40394/251‑1 (MQ=255)
                                                                                                                                    cagcagGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCAt                                                                                                            >  1:265746/1‑251 (MQ=255)
                                                                                                                                           cGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTaa                                                                                                      >  1:7217/1‑250 (MQ=255)
                                                                                                                                            gAATGCTCGGGCGTTCAATTAATGAGGAAAAAGCCAGGT‑AGTTAGCTATCCAGTTTTGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACAGGAACCGGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGAGGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCATCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAg                                                                                                    <  2:363786/251‑1 (MQ=255)
                                                                                                                                                   cTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATCGTGTCGACTTCTTCATCGCTAAAGTGGGTGc                                                                                             >  2:151561/1‑251 (MQ=255)
                                                                                                                                                   cTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTAATCCAGCTTCGGCTAGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAATTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCCCCGATAGAGTCGACTTCTTCATCGCTAACGTGGGTGc                                                                                             >  1:351640/1‑251 (MQ=255)
                                                                                                                                                        gTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGCTTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCTCCACGCGCCGGAAAATGCCCAAATGCCT‑CAGTCACCAGACCACCGATAGTGTCAAACTCTTAATCCCTAAAGTGGGTGCCGacc                                                                                        >  2:261669/1‑249 (MQ=255)
                                                                                                                                                          tCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                                                                       <  1:190289/250‑1 (MQ=255)
                                                                                                                                                                aaTGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGttt                                                                                                                                                                               >  1:193650/1‑157 (MQ=255)
                                                                                                                                                                aaTGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGttt                                                                                                                                                                               <  2:193650/157‑1 (MQ=255)
                                                                                                                                                                 aTGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTg                                                                               >  1:242420/1‑251 (MQ=255)
                                                                                                                                                                     ggCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGt                                                                           >  1:212816/1‑251 (MQ=255)
                                                                                                                                                                      gCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTc                                                                          >  1:221179/1‑251 (MQ=255)
                                                                                                                                                                        aaaaGCCATGT‑AGTTATCTCTCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCCCTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGACACGCGCCGGAAGATGCCCAAATGACTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTAATCGCAAAAGTGGGTGCCGCACCCCTCGCTGAAGTCtt                                                                        >  2:378296/1‑251 (MQ=255)
                                                                                                                                                                        aaaaGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCtt                                                                        <  1:370770/251‑1 (MQ=255)
                                                                                                                                                                              cATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCAACCAGCTTCGGCTGGGGTGAGGCATCCGGGATTTTAACATGAACCTGAATAAG‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATgg                                                                  <  2:102099/251‑1 (MQ=255)
                                                                                                                                                                               aTGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGACCGCTTCGTTGAAGTCTTAAATgg                                                                  >  1:260035/1‑250 (MQ=255)
                                                                                                                                                                                  t‑AGTTATCTCTCACGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGCTTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCACCACTTTTCACTGGTAACCACCGATGTCGATAGGTTCGCCACGCGCCGGAAGATGACAAAATGCCTGCA‑TCACCCGACCACCGATAGTGTCTACTTCTTCATCCCTAACATTGGTGCCGAACGCCTCGTTGAAGGCTTCactgacacc                                                              >  2:125218/1‑243 (MQ=255)
                                                                                                                                                                                       tatctatcCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGCTTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCC‑TCACCAGACCACCGCTAGTGTCGCCTTCTTCCTCGCTAAAGTGGGTCCCGACCGCTTCGTTGAAGTCTTCACTGGAACCCAGt                                                          >  2:114299/1‑251 (MQ=255)
                                                                                                                                                                                       tatctatcCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGCCCACCGCTAGGGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCACTGGAAGcacgt                                                          >  2:200228/1‑247 (MQ=255)
                                                                                                                                                                                        ttctatcCAGTTTCTGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCACCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGt                                                          <  2:356683/249‑1 (MQ=255)
                                                                                                                                                                                        atctatcCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTg                                                         >  1:108553/1‑251 (MQ=255)
                                                                                                                                                                                         tctatcCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGACCCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCCCCAGACCACCGCTAGTGTCGACTTCTTCATCGCTCAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTgc                                                        >  2:342554/1‑251 (MQ=255)
                                                                                                                                                                                            atcCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCTATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTgc                                                        <  1:340279/248‑1 (MQ=255)
                                                                                                                                                                                            atcCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCAGCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTgc                                                        <  2:139310/248‑1 (MQ=255)
                                                                                                                                                                                                aGTTTCGGTTTTATTCATCCAGCTTCGGCTGGTGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACg                                                 <  1:332582/251‑1 (MQ=255)
                                                                                                                                                                                                           tATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGGATg                                      >  1:81644/1‑251 (MQ=255)
                                                                                                                                                                                                              tCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTTGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATgacg                                   >  2:86524/1‑251 (MQ=255)
                                                                                                                                                                                                              tCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATgacg                                   >  1:291628/1‑251 (MQ=255)
                                                                                                                                                                                                              tCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTAAATGGAAGCCAGTGCGCGTACGGTCCAGGTATgacg                                   >  1:237557/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACt                                <  1:151561/250‑1 (MQ=255)
                                                                                                                                                                                                                         cGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGAc                        <  2:212816/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGCCTCAGCTGACg                       >  2:342761/1‑251 (MQ=255)
                                                                                                                                                                                                                            cTGGGGTGAGTCATCCGGGCTTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGCTAGTGTCGACTTCTTCATCGCTAAAGTGGATGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTACAGGTATGACGACTCAGCTGACgt                      >  2:349746/1‑249 (MQ=255)
                                                                                                                                                                                                                             tGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCCTCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGAAGGa                     >  1:250196/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aGTCATCCGGGATTTTAACATGAACCTGACTAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATCGAGTCGACTTCTTCATCGCTACAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGAAGACTCAGCTGCAGGAAGTCGaa              >  2:218217/1‑249 (MQ=255)
                                                                                                                                                                                                                                      tCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAg                                                                                                                           >  1:64658/1‑139 (MQ=255)
                                                                                                                                                                                                                                      tCATCCGGGATTTTAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAg                                                                                                                           <  2:64658/139‑1 (MQ=255)
                                                                                                                                                                                                                                        aTCCGGGATTTTAACATGAACCTGAATAAC‑TCGCCGACTGTCGGACCGCGCCACTTTGAACTGGTAACCGTCGATGTCGAGAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCAGCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCCTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATatcat         <  2:73585/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ttttAACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCttcttc  <  2:242420/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
AATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGT‑AGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAAT‑ACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCA‑TCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCTTCTTC  >  NC_000913/690705‑691193

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