Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,764 G→A S181S (TCC→TCT)  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,7640GA87.8% 104.3 / 5.0 49S181S (TCC→TCT) arcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (26/17);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

GGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGT  >  NC_000913/4639523‑4640008
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
ggAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCTCGTCTACAGTACGGTCGTGTGGTTTCAGCTCACTGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAAt                                                                                                                                                                                                                                                 <  1:497951/250‑1 (MQ=255)
       agggtGAAAAATAAATACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCGTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATTTTTACGAATACGGCGGATCGTCACGTCTACAGTACGTTCGTGTGGTTTCAGCTCTCGGCCGGTCAGTTGCTGCAGCAGTTCAGCACGAGACTGAATTTTGCCTg                                                                                                                                                                                                                                         <  1:187269/247‑1 (MQ=255)
            gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGttt                                                                                                                                                                                                                                    >  1:482474/1‑251 (MQ=255)
            gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGttt                                                                                                                                                                                                                                    <  1:114320/251‑1 (MQ=255)
                     aaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAg                                                                                                                                                                                                                           <  1:37466/251‑1 (MQ=255)
                         ggCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCTTGGGTTTTCACAGAAGTGaa                                                                                                                                                                                                                       <  1:63999/251‑1 (MQ=255)
                            gcTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGc                                                                                                                                                                                                                     >  1:264188/1‑250 (MQ=255)
                              tAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACCGAAGTGAAGCAT‑g                                                                                                                                                                                                                 >  1:217239/1‑251 (MQ=255)
                                  aaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑Ggcgc                                                                                                                                                                                                             >  1:271737/1‑251 (MQ=255)
                                              ttttGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTCCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTg                                                                                                                                                                                                  <  1:190141/250‑1 (MQ=255)
                                              ttttGAAGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTACCGTCTACAGTAAGGTCGTGTGGGTTTAGCTCACGGACGGTCCTTTTCTTCAGCCGATCAGCACCAGGCTGAACTTTGCCTGGGGTTTTACCAAAGTGAAGCAT‑ggaggggaaatccctgc                                                                                                                                                                                                 >  2:473346/1‑236 (MQ=255)
                                                  gACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCgg                                                                                                                                                                                             >  1:119919/1‑251 (MQ=255)
                                                        ggtAAAGCCGATTAATCTTCCAGACCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCt                                                                                                                                                                                        <  2:406261/250‑1 (MQ=255)
                                                           aaaGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGt                                                                                                                                                                                     >  1:326595/1‑250 (MQ=255)
                                                             ggCCGATTAATCTGCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACt                                                                                                                                                                                  <  2:119919/250‑1 (MQ=255)
                                                                        tttCCAGATCCCCCCAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATc                                                                                                                                                                        <  2:271737/249‑1 (MQ=255)
                                                                                  accacAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCa                                                                                                                                                            <  1:484733/247‑1 (MQ=255)
                                                                                     acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCCCAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCaacg                                                                                                                                                         >  1:495230/2‑251 (MQ=255)
                                                                                      cAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCTCAGAAGTGAAGCAT‑GGCGCGGAACTCGCTTCGCGGCAGCTTGTACTTCTCGCCATCA‑GGGCCGATCaaca                                                                                                                                                         >  1:69549/1‑249 (MQ=255)
                                                                                         aaGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGGGAAGCAT‑GGCGCGGAACTCTCTGCGCGGCAGCTTGTTCTGCTCGCCATCA‑GGGCCGCTCaaagaacg                                                                                                                                                     >  2:493189/1‑251 (MQ=255)
                                                                                            cGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGACTACGGCGGATCGTCACGTCTACAGTACGGGCGGGTGGGTTCAGCTCACGGGCGGTCATTTTCTTCAGCAGTTCCGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATTGGC‑CGGAACTCGCTGCGCGGCAACTTGTACTTCCCGCCATCT‑GGGCCGATCaccgaaaggat                                                                                                                                                  >  2:297535/1‑246 (MQ=255)
                                                                                               aaaCCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGtt                                                                                                                                               <  2:396217/250‑1 (MQ=255)
                                                                                                      cACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAACGTCACGTCTCAAGTACGGTCGTGTGGTTTCAACTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGGTTTCACAGAAGGGAAGCCT‑GGAGCGGAACTCGCTGCGCGGCAGCTTGTTCTGCTCGCCCTAC‑GGGCAGATCACAGAACGGCTGTTGgtgttc                                                                                                                                        >  2:423915/1‑249 (MQ=255)
                                                                                                                              ccggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑t                                                                                                                 >  1:293955/1‑249 (MQ=255)
                                                                                                                               cggcgtatccggcgtaGATTCGCAATGTTTACGAATACGGCGGCTCGTCACGTCTACAGTACGGGCGTGTGGTTTCAGCTCACGGCCGGTCCTTTTCTTCAGAAGTTCAGCACAGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGGACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCAGTTGATGTCCACGTCCCCACCATACAACTTG‑TAGc                                                                                                              >  2:309927/1‑251 (MQ=255)
                                                                                                                                                        tGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCCGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTAACAGAAGTGAAGCAT‑GGCGCGGCACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGACGTCCAGTTCCCAACCATTGTACTTG‑TAGCTTTCAACGCTACGAAGttcttccct                                                                                     >  1:281554/1‑248 (MQ=255)
                                                                                                                                                            tACGAATACGGCGGATCGTCCCGTCTACAGAACGGTCGGGTGGTTTCAGATCACGGCCGGTCATTTTCTTCAGCAGTTCCGAAAGAGACTGAGTTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCA‑CATGGGCCGATCAACGAACGGCGGATGATGTCCCGTTCCCCACCATTGAAACTG‑GAGCTGTCAACGCTACGACGTTCTTCGcgggac                                                                                 >  2:341375/1‑246 (MQ=255)
                                                                                                                                                            tACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGAc                                                                                  <  2:159240/250‑1 (MQ=255)
                                                                                                                                                              cGAATACGGCGGATCGTCACGTATACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGCATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGCACGGCTGTTGATGTCCAGTTCCCAACCATTAAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGt                                                                               >  1:460936/1‑251 (MQ=255)
                                                                                                                                                                 aTACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTAcc                                                                            >  1:460011/1‑251 (MQ=255)
                                                                                                                                                                         gATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTc                                                                     >  1:290659/1‑250 (MQ=255)
                                                                                                                                                                             gTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGt                                                                >  1:305426/1‑251 (MQ=255)
                                                                                                                                                                              tCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGCTTCATGGt                                                                >  1:187402/1‑250 (MQ=255)
                                                                                                                                                                                      aCAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGAc                                                         >  1:182848/1‑249 (MQ=255)
                                                                                                                                                                                      aCAGTACGGACGTGAGGATTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGATTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTAGCCATAA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTACACACAATTGAACTTGATA‑CTTTAAAAGCTACGGCGTTCTTAGCTGACAGTACCCAGATTCCTGGGACgagtc                                                         >  2:331660/1‑245 (MQ=255)
                                                                                                                                                                                       aaGTACGGGCGTGGGGGTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAgc                                                      <  2:182848/250‑2 (MQ=255)
                                                                                                                                                                                             ggTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCATCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCTATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGTCAGCAGGTTg                                                <  1:371925/251‑1 (MQ=255)
                                                                                                                                                                                                cGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGACCGGCTGTTGATTTCCAGTTCCCAACCGTTGAACTTG‑TTGCTTTCAACGCTACGACGTTCTTCGCTGACCGTACCCATATTCATGGTACGGGACAGCCGGTTGCgg                                             >  2:115121/1‑250 (MQ=255)
                                                                                                                                                                                                 gtgtggTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTcacaa                                                                                                                                                                                                                              >  2:382756/5‑75 (MQ=255)
                                                                                                                                                                                                 gtgtggTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAg                                                                                                                                                                                                                              <  1:382756/72‑1 (MQ=255)
                                                                                                                                                                                                   tgggTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATGTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGc                                           <  2:115596/248‑1 (MQ=255)
                                                                                                                                                                                                       tttCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGATCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTGCCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGa                                       >  1:186863/1‑250 (MQ=255)
                                                                                                                                                                                                         tCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAAt                                     <  1:441974/250‑1 (MQ=255)
                                                                                                                                                                                                                                 tcAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGt            >  1:486214/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cagcagTTCAGCACGAGACTGAATTTTGCCTGGGTTGTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGt            <  2:293955/250‑1 (MQ=255)
                                                                                                                                                                                                                                   agcagTTCAGCACGAGCCTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTACAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTAAGTTCACGCGGGTTGAACGGttt          >  1:243880/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cagTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGAACTGCTCGCCATCA‑GGGACGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGAACCCAGATTCATGGTACGGGGCAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTgg        >  1:250004/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTTCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATg   <  2:6052/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCCTTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCCCGAATCGTCAGTTCACGCGGGTTGATCGGTTTGGGGATGt  >  2:230616/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
GGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGT  >  NC_000913/4639523‑4640008

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