Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,608 G→A C233C (TGC→TGT 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,6080GA94.1% 142.9 / NA 51C233C (TGC→TGTarcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/16);  minor base C (3/0);  total (35/16)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCG  >  NC_000913/4639365‑4639851
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
cAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGaa                                                                                                                                                                                                                                                <  2:458393/249‑1 (MQ=255)
cAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAg                                                                                                                                                                                                                                               <  1:123981/250‑1 (MQ=255)
        ccGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTACTCTTCCAGATCACCACAGAAGCGATAACCt                                                                                                                                                                                                                                      >  1:518856/1‑251 (MQ=255)
                 cgcTGTTGTTGGGTACCACTGACCTGGTGGTGGCGATCGCCGCAATGATGTATGAGACCCGAAATATCGTCTGGTATGCGTTTTCAATGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGCCGCTCACTTACGTATCTGGAAATCAGGTTGAAAACTAAAAACGGCGCTAACAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTACTCTTCCGGATCACACCAGAAGCGATacactccacagtgg                                                                                                                                                                                                                             >  2:567781/1‑238 (MQ=255)
                              aaCCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGCAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGTTGGTGGCgatg                                                                                                                                                                                                                 >  2:179229/1‑250 (MQ=255)
                              aaCCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCgatg                                                                                                                                                                                                                 >  1:539101/1‑250 (MQ=255)
                                ccAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGAtt                                                                                                                                                                                                              <  1:527758/251‑1 (MQ=255)
                                    ggtgCTGGTGGTGGGGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTtccg                                                                                                                                                                                                          <  2:426933/250‑1 (MQ=255)
                                      tgCTGGTGGGGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGCCTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGGGGCGATGATTtccggc                                                                                                                                                                                                        >  2:532058/1‑251 (MQ=255)
                                                ggCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTtccggcgtatccggc                                                                                                                                                                                               <  2:85577/250‑1 (MQ=255)
                                                   gcTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCCCAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTtccggcgtatccggcgt                                                                                                                                                                                             >  1:52521/1‑249 (MQ=255)
                                                      ggatggCATGATGGTTGTGAGTCGTAATATCGACTGGTATGGGTGGTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAAGTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGTGCGAAAAAGCGCCGGTTTTTTTGTCGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGTGATAACCGTGAGCGTGGATGGTGGCGAAGGTGTCCGGCCTAGCCGGCGTAGAtt                                                                                                                                                                                        <  2:242157/247‑1 (MQ=255)
                                                                   tttGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACg                                                                                                                                                                           >  1:287072/1‑251 (MQ=255)
                                                                        gactcGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAc                                                                                                                                                                      <  2:281787/247‑1 (MQ=255)
                                                                         actcGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcg                                                                                                                                                                     <  2:266957/248‑1 (MQ=255)
                                                                                                 ttttCACTGCCGAAAATGAAATCCAGTAAAGAAGTTACAATGGACGCTCAAGTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTATAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAACGATAAACTTCACCGTGGATGTTGGAGATGATTTCCGGCGTATCCGGCGTAGATTCGAACTGTTTACGAATACTGCGGATCGTCACGTGTACAGtgcg                                                                                                                                             >  2:620052/1‑248 (MQ=255)
                                                                                                  tttCACTGCCGAACCTGAAAGCCAGTAAAGAAGTTGCCATGGACGATCAGTTACGTATCTGGACATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACgg                                                                                                                                            >  2:90315/1‑251 (MQ=255)
                                                                                                                       ccAGTAAAGACGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTAGTGTGGTTTCCGCTCACgg                                                                                                                       >  2:257849/1‑251 (MQ=255)
                                                                                                                           tAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCCCGTCTACAGTAAGGTCGTGTGGTTTCAACTCACAGCCgg                                                                                                                   >  2:354842/1‑251 (MQ=255)
                                                                                                                            aaaGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGt                                                                                                                  >  1:302947/1‑251 (MQ=255)
                                                                                                                             aagaagTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTc                                                                                                                 <  1:591048/251‑1 (MQ=255)
                                                                                                                                         tggACGATCACTTACGTATCTGGGAATAAGGTTGAAAAATACAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGCTGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTAACGTACGGTAGAGTGGTTTCCGCTCACGGCCGGACAtttttattagc                                                                                                     >  2:360348/2‑244 (MQ=255)
                                                                                                                                            acgatcaattACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTcagcagg                                                                                                  >  2:217217/9‑250 (MQ=255)
                                                                                                                                             cgatcaattACGTATCTGGAAATAAGGTTTAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGAGCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGtt                                                                                                 <  1:362985/244‑1 (MQ=255)
                                                                                                                                             cgatcaattACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGtt                                                                                                 >  1:367522/8‑251 (MQ=255)
                                                                                                                                                  aattACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCTGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGc                                                                                             >  1:282906/3‑250 (MQ=255)
                                                                                                                                                   attACGTATCTGGAAGTAAGGTTGAAATATAAAAACGGCGCTCAAAAGCGCCGTTTTTTTTGACGGAGGTAAAGACGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCATCACGTCTACAGTACGGTAGTGTGGTTTCAGCTCACGGACGGTCAtttcattccgtcgttcagcac                                                                                           >  2:424203/2‑233 (MQ=255)
                                                                                                                                                    ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCAcc                                                                                                                                                                                                                                 >  1:505482/1‑116 (MQ=255)
                                                                                                                                                    ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCAcc                                                                                                                                                                                                                                 <  2:505482/116‑1 (MQ=255)
                                                                                                                                                                 aaTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTg                                                                             <  1:471639/251‑1 (MQ=255)
                                                                                                                                                                 aaTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGGTTTTTTTGACGGGGGTAAAACCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGGCTACAGTACGGGCGTGTGGTTTCAGCTCACTGCCGGTCATTTTCTTGCGCAGTTCCGCACAGGCCTGCATTTTg                                                                             >  2:160288/1‑251 (MQ=255)
                                                                                                                                                                           aaaaaTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGAAGGAGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCCGCTCACGGCCGGGCATTTTCTTCAGCAGTTCAGCACGAGACTGCATTTTGCCTGGGtttt                                                                   >  2:558927/1‑251 (MQ=255)
                                                                                                                                                                             aaaTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGAGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCCACAGTACGGTCGTGTGGGTTCAGCTCACGGCCGGGCACTTTCTTCCACCGTTCCGCGAGCGACTGAATTTTGCCTGGGttttt                                                                  >  2:40207/1‑249 (MQ=255)
                                                                                                                                                                             aaaTAAAAACGGCGCTAAAAAACGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGAGGCGATGATTTCCGGCGAATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATAGTCCCGTCTACAGTACAGTCGGGTGGTTTCAGCTCACGGCCGGTCATTTTCTGCAGCAGTTAAGCACGAGACTGCATTTTGCCTGGGTTTTc                                                                  >  1:355895/1‑250 (MQ=255)
                                                                                                                                                                                    aaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAgg                                                          >  2:452962/1‑250 (MQ=255)
                                                                                                                                                                                      cGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCATCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACCGAAGTg                                                         >  2:479992/1‑250 (MQ=255)
                                                                                                                                                                                            tAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCCCAGAAGTGAAGCATg                                                  >  2:159006/1‑251 (MQ=255)
                                                                                                                                                                                              aaaaGCGCCGTCTTTCTTGACGATGGTAAAGCCGATTAAACTTCCAGATCACCACAGAAGCGATACCATTCAACGTGGATGGTGGCGATGATTTTCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATCCGGCGGATCGTCACGGCTACAGTAAGGTTGAGTGGTTTAAGCGAACGGCCGGGCATTTTCTTCCGCCGTTCAGCCCGAGACAGAATTTTGCGTGGGATTTCACAGAAGGGAAGCAGGgc                                                >  2:520479/1‑251 (MQ=255)
                                                                                                                                                                                               aaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGgcc                                               >  2:188763/1‑250 (MQ=255)
                                                                                                                                                                                                               tGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTAAGGTCGTGTGGTTTCAGCTCACGGCCGGGCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCCCTGcgcg                               >  2:182313/1‑251 (MQ=255)
                                                                                                                                                                                                                        tAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTa                                                                                                                                               <  1:326516/130‑1 (MQ=255)
                                                                                                                                                                                                                        tAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAAACTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTa                                                                                                                                               >  2:326516/1‑130 (MQ=255)
                                                                                                                                                                                                                         aaaGCCGATTAATCTTCCAGATCACCACCGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTa                     <  2:332873/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACCGAAGTGAAGCCTGGCGCGGAACTCGCTGCGCGGCCGCTTGTa                     >  2:522579/1‑250 (MQ=255)
                                                                                                                                                                                                                           aGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACt                   <  1:479992/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGcc            >  1:362214/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaTCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTAAGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGAATGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAAATCGCTGCGCGGCAGCTTGTACTGCTCGcc            >  2:12451/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ccAGAGCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGATTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAgg      <  1:341032/250‑1 (MQ=255)
                                                                                                                                                                                                                                         acAGAGCACCCCAGAAGCGATAACCTTCACCGTGGCTGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGAAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTAATGCTCGCCATCAggg     >  2:418165/2‑251 (MQ=255)
                                                                                                                                                                                                                                          cAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGGTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGc    >  1:369327/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aTCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCAGGCGTATCCGGCGTAGCTTCGAAATGTTTACGAATACGGCGGATCGTCCCGTCTACAGTACGGTTGTGTGGATTCAGCTCCAGGCCGGTCATTTTCTTCAGCTGTTAAGAACGAGAATGAATTTTGCCTGGGTTTTAACAGGAGTGAAGCATGGCGCGGCACTCGCTGGCCGGCAGGTTGTACTGCTCGCCATCAGGGCCg  >  2:426478/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
CAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCG  >  NC_000913/4639365‑4639851

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