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→TCTarcADNA‑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

GGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGT  >  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)
                              tAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACCGAAGTGAAGCATg                                                                                                                                                                                                                 >  1:217239/1‑251 (MQ=255)
                                  aaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGgcgc                                                                                                                                                                                                             >  1:271737/1‑251 (MQ=255)
                                              ttttGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTCCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTg                                                                                                                                                                                                  <  1:190141/250‑1 (MQ=255)
                                              ttttGAAGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTACCGTCTACAGTAAGGTCGTGTGGGTTTAGCTCACGGACGGTCCTTTTCTTCAGCCGATCAGCACCAGGCTGAACTTTGCCTGGGGTTTTACCAAAGTGAAGCATggaggggaaatccctgc                                                                                                                                                                                                 >  2:473346/1‑236 (MQ=255)
                                                  gACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCgg                                                                                                                                                                                             >  1:119919/1‑251 (MQ=255)
                                                        ggtAAAGCCGATTAATCTTCCAGACCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCt                                                                                                                                                                                        <  2:406261/250‑1 (MQ=255)
                                                           aaaGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGt                                                                                                                                                                                     >  1:326595/1‑250 (MQ=255)
                                                             ggCCGATTAATCTGCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACt                                                                                                                                                                                  <  2:119919/250‑1 (MQ=255)
                                                                        tttCCAGATCCCCCCAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATc                                                                                                                                                                        <  2:271737/249‑1 (MQ=255)
                                                                                  accacAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCa                                                                                                                                                            <  1:484733/247‑1 (MQ=255)
                                                                                     acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCCCAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCaacg                                                                                                                                                         >  1:495230/2‑251 (MQ=255)
                                                                                      cAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCTCAGAAGTGAAGCATGGCGCGGAACTCGCTTCGCGGCAGCTTGTACTTCTCGCCATCAGGGCCGATCaaca                                                                                                                                                         >  1:69549/1‑249 (MQ=255)
                                                                                         aaGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGGGAAGCATGGCGCGGAACTCTCTGCGCGGCAGCTTGTTCTGCTCGCCATCAGGGCCGCTCaaagaacg                                                                                                                                                     >  2:493189/1‑251 (MQ=255)
                                                                                            cGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGACTACGGCGGATCGTCACGTCTACAGTACGGGCGGGTGGGTTCAGCTCACGGGCGGTCATTTTCTTCAGCAGTTCCGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATTGGCCGGAACTCGCTGCGCGGCAACTTGTACTTCCCGCCATCTGGGCCGATCaccgaaaggat                                                                                                                                                  >  2:297535/1‑246 (MQ=255)
                                                                                               aaaCCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGtt                                                                                                                                               <  2:396217/250‑1 (MQ=255)
                                                                                                      cACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAACGTCACGTCTCAAGTACGGTCGTGTGGTTTCAACTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGGTTTCACAGAAGGGAAGCCTGGAGCGGAACTCGCTGCGCGGCAGCTTGTTCTGCTCGCCCTACGGGCAGATCACAGAACGGCTGTTGgtgttc                                                                                                                                        >  2:423915/1‑249 (MQ=255)
                                                                                                                              ccggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGt                                                                                                                 >  1:293955/1‑249 (MQ=255)
                                                                                                                               cggcgtatccggcgtaGATTCGCAATGTTTACGAATACGGCGGCTCGTCACGTCTACAGTACGGGCGTGTGGTTTCAGCTCACGGCCGGTCCTTTTCTTCAGAAGTTCAGCACAGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGGACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCAGTTGATGTCCACGTCCCCACCATACAACTTGTAGc                                                                                                              >  2:309927/1‑251 (MQ=255)
                                                                                                                                                        tGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCCGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTAACAGAAGTGAAGCATGGCGCGGCACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGACGTCCAGTTCCCAACCATTGTACTTGTAGCTTTCAACGCTACGAAGttcttccct                                                                                     >  1:281554/1‑248 (MQ=255)
                                                                                                                                                            tACGAATACGGCGGATCGTCCCGTCTACAGAACGGTCGGGTGGTTTCAGATCACGGCCGGTCATTTTCTTCAGCAGTTCCGAAAGAGACTGAGTTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCACATGGGCCGATCAACGAACGGCGGATGATGTCCCGTTCCCCACCATTGAAACTGGAGCTGTCAACGCTACGACGTTCTTCGcgggac                                                                                 >  2:341375/1‑246 (MQ=255)
                                                                                                                                                            tACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGAc                                                                                  <  2:159240/250‑1 (MQ=255)
                                                                                                                                                              cGAATACGGCGGATCGTCACGTATACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGCATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGCACGGCTGTTGATGTCCAGTTCCCAACCATTAAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGt                                                                               >  1:460936/1‑251 (MQ=255)
                                                                                                                                                                 aTACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTAcc                                                                            >  1:460011/1‑251 (MQ=255)
                                                                                                                                                                         gATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTc                                                                     >  1:290659/1‑250 (MQ=255)
                                                                                                                                                                             gTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGt                                                                >  1:305426/1‑251 (MQ=255)
                                                                                                                                                                              tCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGCTTCATGGt                                                                >  1:187402/1‑250 (MQ=255)
                                                                                                                                                                                      aCAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGAc                                                         >  1:182848/1‑249 (MQ=255)
                                                                                                                                                                                      aCAGTACGGACGTGAGGATTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGATTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTAGCCATAAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTACACACAATTGAACTTGATACTTTAAAAGCTACGGCGTTCTTAGCTGACAGTACCCAGATTCCTGGGACgagtc                                                         >  2:331660/1‑245 (MQ=255)
                                                                                                                                                                                       aaGTACGGGCGTGGGGGTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAgc                                                      <  2:182848/250‑2 (MQ=255)
                                                                                                                                                                                             ggTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCATCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCTATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGTCAGCAGGTTg                                                <  1:371925/251‑1 (MQ=255)
                                                                                                                                                                                                cGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGACCGGCTGTTGATTTCCAGTTCCCAACCGTTGAACTTGTTGCTTTCAACGCTACGACGTTCTTCGCTGACCGTACCCATATTCATGGTACGGGACAGCCGGTTGCgg                                             >  2:115121/1‑250 (MQ=255)
                                                                                                                                                                                                 gtgtggTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTcacaa                                                                                                                                                                                                                              >  2:382756/5‑75 (MQ=255)
                                                                                                                                                                                                 gtgtggTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAg                                                                                                                                                                                                                              <  1:382756/72‑1 (MQ=255)
                                                                                                                                                                                                   tgggTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATGTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGc                                           <  2:115596/248‑1 (MQ=255)
                                                                                                                                                                                                       tttCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGATCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTGCCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGa                                       >  1:186863/1‑250 (MQ=255)
                                                                                                                                                                                                         tCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAAt                                     <  1:441974/250‑1 (MQ=255)
                                                                                                                                                                                                                                 tcAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGt            >  1:486214/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cagcagTTCAGCACGAGACTGAATTTTGCCTGGGTTGTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGt            <  2:293955/250‑1 (MQ=255)
                                                                                                                                                                                                                                   agcagTTCAGCACGAGCCTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTACAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTAAGTTCACGCGGGTTGAACGGttt          >  1:243880/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cagTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGAACTGCTCGCCATCAGGGACGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGAACCCAGATTCATGGTACGGGGCAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTgg        >  1:250004/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTTCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATg   <  2:6052/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCCTTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCCCGAATCGTCAGTTCACGCGGGTTGATCGGTTTGGGGATGt  >  2:230616/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
GGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGT  >  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: