Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,719 C→T P196P (CCG→CCA)  arcA ← DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,7190CT94.4% 146.8 / NA 58P196P (CCG→CCA) arcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (24/27);  minor base G (3/0);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGT  >  NC_000913/4639472‑4639956
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
cGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGt                                                                                                                                                                                                                                                <  1:209645/251‑1 (MQ=255)
  aaaaTGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCAACGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTCCAGTC‑GTGAGGttt                                                                                                                                                                                                                                              >  1:655240/1‑251 (MQ=255)
    aaTGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGCTCGTCACGTCTACCGTACGGTC‑GTGTGGTTTc                                                                                                                                                                                                                                             >  1:680361/1‑250 (MQ=255)
      tGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGGTTCGAAATGTTTACGAATACGGCGGATCGTCACCTCTACAGTACGGGT‑GTGGGGTTTCatg                                                                                                                                                                                                                                          >  2:355967/1‑249 (MQ=255)
      tGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GGGTGGTTTTAGc                                                                                                                                                                                                                                          >  2:231789/1‑251 (MQ=255)
         aaGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCa                                                                                                                                                                                                                                       <  1:356523/251‑1 (MQ=255)
         aaGCCAGTAAAGAAGGTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCa                                                                                                                                                                                                                                       <  2:613519/251‑1 (MQ=255)
                               ggACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTCAAAAGCGCCGTTTTTTTTGACGGTGGAAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCAAGTCAACAGTACGGTC‑GTGTGGATTCAGCTAACGGCCGG‑TCATTttctttggcc                                                                                                                                                                                                                >  1:430142/1‑246 (MQ=255)
                                  cgatcaattACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGtt                                                                                                                                                                                                             <  2:411397/244‑1 (MQ=255)
                                  cgatcaattACGTAGCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGtt                                                                                                                                                                                                             <  1:524876/244‑1 (MQ=255)
                                     tcaattACGTATCTGGAAATAAGGTTGAAACATAAAAACGGCGCTAAAACGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGAGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGAGTAGATTCGAAATGTTTACGAATACGGCGGCCAGTCACGTCTTCAGTCCGGTCTGTG‑GGTTTCAGCTCACGGCCGA‑TATTTTTCTTcagaggttca                                                                                                                                                                                                           >  2:314456/5‑243 (MQ=255)
                                         ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACg                                                                                                                                                                                                      >  1:418390/1‑251 (MQ=255)
                                            cGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGTTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACgaga                                                                                                                                                                                                   <  1:175628/251‑4 (MQ=255)
                                                    gAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAAttt                                                                                                                                                                                           <  1:309043/251‑1 (MQ=255)
                                                           ggTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑gg                                                                                                                                                                                   <  1:644765/251‑1 (MQ=255)
                                                               gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑GCATTTTCTTCAGCAGTTCAGCCCGCGGCTGAATTTTTACT‑GGGttt                                                                                                                                                                               >  2:184132/1‑251 (MQ=255)
                                                               aaaaaaTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGttt                                                                                                                                                                               <  2:144326/250‑1 (MQ=255)
                                                                    aTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGGGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTcac                                                                                                                                                                           >  1:554439/1‑250 (MQ=255)
                                                                       aaaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGGCTACAGGACGGGC‑GTGGGGTTTCATCCCACGGCCGGATC‑TTTTCTTCCGCCGTTCAGCCCGGGACTGAATTTTGGCT‑TGGTTTTcaccga                                                                                                                                                                        >  2:134791/1‑247 (MQ=255)
                                                                                  aaaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATgg                                                                                                                                                            <  1:507914/251‑1 (MQ=255)
                                                                                   aaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTc                                                                                                                                                                                                                                             <  1:240354/171‑1 (MQ=255)
                                                                                     aaGCGCCGTTTTTTTTGCCGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTAGGGTTTT‑ACAGAAGTGAAGCATGgcgc                                                                                                                                                         >  1:537919/1‑251 (MQ=255)
                                                                                        cgccggtgttttttACGGTGTTACAGCCGTTTAATTTGCCAGGTCACCACAGCATCGTTATCCCTCACCTTGGATGGTGGCGATGATTTTCGGCGTATCCGGCGTAGATTCAAAATTTTTCCTAATCCGGCGGATCGTCCCGTCTTCAGTACTGTC‑CTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GCGTTTTCACTGAAGTGAAGCATGGCGCGGa                                                                                                                                                      <  2:176875/243‑1 (MQ=255)
                                                                                                        ggtggtAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGGTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGGGAAGCATGGCGCGGAACTAGCTGggcggc                                                                                                                                        >  2:524625/1‑249 (MQ=255)
                                                                                                                aGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTc                                                                                                                                                                                                                                        <  1:364380/147‑1 (MQ=255)
                                                                                                                   cGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTATATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTCTCACAGAAGTGAAGCATg                                                                                                                                                             <  2:97576/217‑1 (MQ=255)
                                                                                                                   cGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATg                                                                                                                                                             >  1:97576/1‑217 (MQ=255)
                                                                                                                    gATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTc                                                                                                                          >  1:331189/1‑251 (MQ=255)
                                                                                                                       tAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGGCCCGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTTAGCAGTTCAGGACGAGACTGAATTTTTCCT‑GGGTTTTTACCGGAGGTAAGCATGTCGCGGGACACACTGCGCGGcggtttttacgcctcgcc                                                                                                                       >  2:418005/1‑234 (MQ=255)
                                                                                                                                aGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGcc                                                                                                             >  1:641416/1‑251 (MQ=255)
                                                                                                                                        acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCaacg                                                                                                     <  1:484887/250‑1 (MQ=255)
                                                                                                                                                gATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGACACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCCAGCGACTGAATTTTGCCT‑GGGTTTTCCCAGAAATGAAGCATGGCGCGGGACTCGCTGCGCGGCCGCTTGTACTGCTCGCCCT‑CAGGGCCGGTCAACGCACGGCTg                                                                                             >  1:375605/1‑251 (MQ=255)
                                                                                                                                                      cTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATg                                                                                       >  1:315858/1‑251 (MQ=255)
                                                                                                                                                              gggaTGGTGGCGATGCTTTCCGGCGTATCCGGCGTAGATTCGAAACGTTTACGAATACGGCGGATCGTCACGTCTTCAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCATCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATTGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTc                                                                               <  1:560877/248‑1 (MQ=255)
                                                                                                                                                               ggatgtgggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCTCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTcc                                                                              <  2:107832/244‑1 (MQ=255)
                                                                                                                                                                    gtggCGATGATTGCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCTG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GTGTTTTCACAGAAGTGAAGCATGGCGCGGAACGCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAAcc                                                                         <  1:609585/251‑1 (MQ=255)
                                                                                                                                                                      ggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAATTCAGCACGAGACTGAATTTTGCCT‑GGGGTTTCACAGACATGAAGCATGGCCAGGAACTCGCTGCGCGGCCGCGTGTACGGATCGCCATTCAGT‑CCGATCAACGTACGGCCGTTGATGTACAGTTCCCCACCAt                                                                       >  2:315418/1‑251 (MQ=255)
                                                                                                                                                                       gCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCAtt                                                                      >  1:558460/1‑251 (MQ=255)
                                                                                                                                                                         gatgatTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTGCGGTC‑GTGTGGCTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGGTCAGCACGGGAATGAATTTTGCCT‑GGGTTTTCACATAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTCCTGATCGCCAT‑CCGGGCCGATCAACGAACGGCTGATTAAGTCCAATTCCCCACCATTg                                                                     >  2:622466/1‑250 (MQ=255)
                                                                                                                                                                                tccggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGt                                                              <  1:524625/250‑1 (MQ=255)
                                                                                                                                                                                   ggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCt                                                         >  1:603927/1‑251 (MQ=255)
                                                                                                                                                                                    gcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTTAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCGCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCtt                                                        <  1:120998/251‑1 (MQ=255)
                                                                                                                                                                                     cggatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGTGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCttt                                                       <  2:641416/248‑1 (MQ=255)
                                                                                                                                                                                           cggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGc                                                 >  1:416708/1‑251 (MQ=255)
                                                                                                                                                                                              cgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTT‑TATGCATTCAACGCTac                                              >  1:540172/1‑251 (MQ=255)
                                                                                                                                                                                                      cGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGGA‑GCTTTAAACGCTACGACGttct                                      >  1:399669/1‑251 (MQ=255)
                                                                                                                                                                                                       gAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGttctt                                     <  1:184132/251‑1 (MQ=255)
                                                                                                                                                                                                        aaaTGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGttcttc                                    <  1:241149/251‑1 (MQ=255)
                                                                                                                                                                                                                aCGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACa                             <  2:280951/250‑1 (MQ=255)
                                                                                                                                                                                                                       cgtcggttcttatcGTCTACAGGTCGGGC‑GTGTGGTTTTAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCTGCTTTTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGTCGTTCTTCGCTGACAGTACCCAg                     <  2:358760/238‑1 (MQ=255)
                                                                                                                                                                                                                       cgggggATCGTCACGTCTACAGTACGGGC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAg                     <  2:244225/247‑1 (MQ=255)
                                                                                                                                                                                                                       cggcggATCGTAACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGATCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGGTGTCCAGATGCCAACCATTGAGCTTGTC‑GCTTTCAACGCTACGACTTTTTTCGCTGCCAGTAccc                       >  2:195702/1‑249 (MQ=255)
                                                                                                                                                                                                                           ggATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTc                 >  1:634926/1‑251 (MQ=255)
                                                                                                                                                                                                                               cGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATg              <  1:612807/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGACCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGATGACAGTACCCCGATTCATGGTAAggg       >  2:569550/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGCAAGTACCCAGATTCATGGTACgggg      >  2:640817/1‑250 (MQ=255)
                                                                                                                                                                                                                                         aCAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAg   <  1:158408/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAgc  >  1:510294/1‑248 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
CGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGT  >  NC_000913/4639472‑4639956

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: ‑