Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,869,637 C→A L87L (CTC→CTA)  cbrA → colicin M resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,869,6370CA90.0% 183.7 / 2.7 70L87L (CTC→CTA) cbrAcolicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base A (39/24);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGAT  >  NC_000913/3869396‑3869880
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cgagtattGGCCGCGGCCCGTCAGGATCGGCGTTGGCACGAAAGGTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGGGGTACTGATGGTGTCAGCACACCCTGCGGCGGTCTGCTGGCACCTGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑AGCGCTAACCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAAc                                                                                                                                                                                                                                             <  2:420606/247‑1 (MQ=255)
cGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGTACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAAc                                                                                                                                                                                                                                             <  1:483945/251‑1 (MQ=255)
 gATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCCCCGGACGCGCCGCGGTCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGGG‑ATCGCTAATCCGCAGATTTTTAGCGTCAAAACCGTCGGCGTCGCCGCCTCGCTCACGCGTAACt                                                                                                                                                                                                                                            >  2:5554/1‑251 (MQ=255)
     attagCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAAACCTGCGGCGGTCTGCTGGCACCGGACGAGCAACGGTCTTTTATTCGCGGTGGGCTGACGCTTCCTGTC‑GAGGGG‑ATCGCTAAACCGCAGAGTTTCAGCGTCAAAACCGTCGACGACGCAGCATCGCGAACGCGTAACTAcc                                                                                                                                                                                                                                         >  2:517929/5‑250 (MQ=255)
       tGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGc                                                                                                                                                                                                                                      >  1:603838/1‑251 (MQ=255)
           cTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCCGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTTCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGGCGCCGCATAGCTAACGCGTAACTACCAGCGa                                                                                                                                                                                                                                    >  1:612418/1‑249 (MQ=255)
            ccGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGc                                                                                                                                                                                                                                 <  1:498263/250‑1 (MQ=255)
                    gCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCaa                                                                                                                                                                                                                         <  1:5554/251‑1 (MQ=255)
                      aGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAata                                                                                                                                                                                                                       <  1:150311/251‑1 (MQ=255)
                          ccGGCGTTGGCACCAAAGTTAGCAGGCAAAATGCAGGTTATCGCGCTGGATAAAAAGCACCAGTGTGGTACTTAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGATGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTaa                                                                                                                                                                                                                  <  2:682983/250‑1 (MQ=255)
                            ggCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATc                                                                                                                                                                                                                <  1:546000/251‑1 (MQ=255)
                              cGTTGGCACGCAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGAGACGCGCTCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑CTCGCTTCTCAGCAGCTTTTCAGCGTCAAAACCGTCGCCGACGCCGCCTAGATAAAGCGTAACTACACGCGAAGCTATATCAATA‑TTAATcac                                                                                                                                                                                                              >  2:209320/1‑249 (MQ=255)
                               gTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGCAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGcc                                                                                                                                                                                                             >  1:680719/1‑251 (MQ=255)
                                 gggCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAATGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCa                                                                                                                                                                                                            <  2:224114/249‑1 (MQ=255)
                                   gCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGCTAAAAAGCACCAGTGTCCTACTGAAGGTTTCCGCAAACCCTGCGGAGGTCTGCTGGCCCCGGACGCGCCGCGCTCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTTAGCGTCCAAACCGCCGACGTCGCCGCATCGGTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATgg                                                                                                                                                                                                         >  1:336347/1‑250 (MQ=255)
                                       gAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTc                                                                                                                                                                                                     <  1:115044/251‑1 (MQ=255)
                                       gAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTAGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGACATGCtttt                                                                                                                                                                                                     >  2:527954/1‑250 (MQ=255)
                                        aacgTTAGCAGGCAAAATGCAGGCGATCGCGCGGGATAACAAGCACCAGTGAGGTATTGAAGGTTTCAGCAAACCCTGCGGCGATCTGCTGGCACCGGACTCGCAGCGATCTTTTATTCGCGATGGGCTCACGCGTCCTGTC‑GATGTG‑ATCGCTCATCCGCAGATTTTCAGCGTCAAAACCGTCGACACCGCCGAATCGCTCACGCGACACTACCAGCAACGCTCTCTCACTA‑TTAATTGCAATGCTTTCg                                                                                                                                                                                                    >  2:487313/4‑251 (MQ=255)
                                           gTTAGCCGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTAAGCAAACCCTGCGGCGGTCTGCTGGCACCGGAAGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAACACCGTCGACGTCCCCGCATCGCTAACGCGGAACTAACAGCGACGCTATATCAATA‑TTAATCGCCATGCTTTCGACt                                                                                                                                                                                                 >  2:302389/1‑251 (MQ=255)
                                                    cAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATg                                                                                                                                                                                         >  1:451335/1‑250 (MQ=255)
                                                    cAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTACTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTTCCAGCGCAGCTATATCAATA‑TTACACGCACTGCTTTCGACTTGGGGATgt                                                                                                                                                                                        >  2:70273/1‑250 (MQ=255)
                                                     aacaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGAACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCCCCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGATAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCCTCGCTCACGCGGAAATACCAGCGCAGCTATATCAATA‑TTACTCGCCCTGCTTTCAGCTTGTGGATgc                                                                                                                                                                                        >  2:405337/4‑249 (MQ=255)
                                                               ggATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATc                                                                                                                                                                             <  2:182225/250‑1 (MQ=255)
                                                                              aaaaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTg                                                                                                                                                                >  1:170726/1‑249 (MQ=255)
                                                                                aaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGt                                                                                                                                                            >  1:571591/1‑251 (MQ=255)
                                                                                  gCACCAGAGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGAGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGAGTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCCATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCACTGATCCCCGCCAGCGTTGAGGttt                                                                                                                                                          >  1:484993/1‑251 (MQ=255)
                                                                                  cacccaggtgggTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCTGCGATCTTTTATTCGCGAGGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCATATTTCCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAAGA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTAGAGGttg                                                                                                                                                          <  2:447226/241‑2 (MQ=255)
                                                                                    aCCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTAt                                                                                                                                                        >  1:378440/1‑251 (MQ=255)
                                                                                          gtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGAc                                                                                                                                                  >  1:170489/1‑251 (MQ=255)
                                                                                                   aaGGCTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGCTCCTTTATTCGCGATGGGCTGACGCTTCCTGCC‑CATGTG‑ATCGCTCATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGAATCGCTCCCCCGTCAATACCAGCGAAGCTATAACAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGCTTATCACGACAGCCTGTGc                                                                                                                                         >  1:36179/1‑251 (MQ=255)
                                                                                                           agcaaccccggggggggTCTGCTGGCCCCGGACGCGATGCGTTCTTTTATTCCCGATGGGCTGCCGCTTCC‑GTCCGATGTG‑ATCCCTAATCCGCAGATTTTCAGCGTCAAAACCGTAGACGTCGCCGCATCGCTACCGCGTAACTCCCAGCGAAGTTATATCAATA‑TTAATCGCATTGCTTTCGACTTGTGGATGAAATCGCCGATCACCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAAt                                                                                                                                 <  2:150353/236‑1 (MQ=255)
                                                                                                              aaaCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCg                                                                                                                                                                         <  1:697246/208‑1 (MQ=255)
                                                                                                              aaaCCCTGCGGCGGTCTGCTGGCAACGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGAGAAGCTATATCAATA‑TTAATCGCCCTGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTg                                                                                                                              >  1:441032/1‑251 (MQ=255)
                                                                                                                       ggcggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGCTGAAATCGCAGATCCCAGCCAGCGTTGAGGTTTATCACGCCAGCCTGTGACGAAAAATATGGCGTggga                                                                                                                      >  2:563475/1‑247 (MQ=255)
                                                                                                                          ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatgt                                                                                                                  >  2:679570/1‑250 (MQ=255)
                                                                                                                          ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatg                                                                                                                   >  1:174010/1‑250 (MQ=255)
                                                                                                                                 tgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATgg                                                                                                           >  1:531817/1‑251 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCtt                                                                                                                                                                                                       >  2:533879/1‑152 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAACACTGGCGTGCGGATGCTAAATGGCATGTc                                                                                                     >  2:278096/1‑250 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑GTAATCGCCATGCtt                                                                                                                                                                                                       <  1:533879/152‑1 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCATTTATCTGCGATGGGCTGATGCTTCCGGGC‑GA‑GTGGATCCCTAATCCGCATATTTTCAGCGACAAAACCGTCGACGCCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCa                                                                                                    <  2:129870/251‑1 (MQ=255)
                                                                                                                                         gACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAt                                                                                                   >  1:173259/1‑251 (MQ=255)
                                                                                                                                                aGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGTTAAATGGCATGTAATTTTTCgg                                                                                            >  2:502779/1‑250 (MQ=255)
                                                                                                                                                  cGATTTTTTATTCGCGATGGGGTGACGCTTCCTGTC‑GATTGG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACGACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTGTCGACTTGTGGAGGAAATCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGc                                                                                          <  1:405337/251‑1 (MQ=255)
                                                                                                                                                    aTCTTTTATTCGCGAAGGGCTGACGCTTCCTGTC‑GAAGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGACGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATTTCATTTTTCGTGccg                                                                                        >  2:390094/1‑249 (MQ=255)
                                                                                                                                                        tttATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGc                                                                                          <  1:521299/245‑1 (MQ=255)
                                                                                                                                                        tttATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACg                                                                                     <  1:319493/250‑1 (MQ=255)
                                                                                                                                                              ggggaggggCGGACGCTCCCTGTT‑GATGTG‑ACCGCTAAACCGCAGATTTTCAGCGACAAACCCGCGGAGGGCGCCCCATCGCTAACGCGTATCTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGTTTTCGATTGGGGGATGAAAGCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCAGGTCATTTTTCGTGCAGACGGCTGGGa                                                                              <  2:628186/245‑1 (MQ=255)
                                                                                                                                                                             cTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgc                                                               >  1:612867/1‑251 (MQ=255)
                                                                                                                                                                               cccTGTC‑GTTTTG‑ATCGCTAAGCCTCAGATTGTCAACTTCAAAACCGCAGCCGTCGCCCCACCGTTACCGCGCAACTACCAGAGAAG‑TTTAACAATAATTAATCGCCATGCTTCCGACTTTGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGAGGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgccc                                                             <  2:496321/250‑1 (MQ=255)
                                                                                                                                                                                  tGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCACAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTAACAGCGACGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGACATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAACCTGTGCCGCAAAATCTGGCGTAAGGAGGATAAATGGCCTGTCATTTTTCGTGCAGACGGCTGGGAgaaccataatcccgccagct                                                          >  2:179713/1‑232 (MQ=255)
                                                                                                                                                                                  tGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTAGTGCAGACGGCTGGGAGCAACATATTAccgcccgtt                                                          >  2:385607/1‑249 (MQ=255)
                                                                                                                                                                                      ggATGTG‑ACCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATc                                                       <  1:385607/250‑1 (MQ=255)
                                                                                                                                                                                       gATGTG‑ATCGCTACTCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTAGTGCAGACGGCTGGGAGCAACATATTACCCACCGCTATCt                                                      >  2:509042/1‑251 (MQ=255)
                                                                                                                                                                                           tg‑ATCGCTAATCCGCAGCTTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGCCAGCCTCTGCCGAACAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGCCGGCTGGGCGCCACATATTACCTCCCGCTATCTGGTc                                                  >  2:363070/1‑251 (MQ=255)
                                                                                                                                                                                            g‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCg                                                 >  1:75880/1‑251 (MQ=255)
                                                                                                                                                                                            g‑ATCGCTAACCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTACCTACCAGCGAAGCTATACCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTc                                                  <  1:140280/250‑1 (MQ=255)
                                                                                                                                                                                                cGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTg                                              >  1:506264/1‑251 (MQ=255)
                                                                                                                                                                                                   tAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTgacg                                           >  2:485194/1‑248 (MQ=255)
                                                                                                                                                                                                    aaTCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGCGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGGGCCgc                                          >  2:129903/1‑250 (MQ=255)
                                                                                                                                                                                                      tCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATg                                        >  2:12515/1‑251 (MQ=255)
                                                                                                                                                                                                        cGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGGTATATCAATA‑TTAATCGCCAGGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgc                                      <  1:658894/251‑1 (MQ=255)
                                                                                                                                                                                                        cGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgc                                      >  1:14406/1‑251 (MQ=255)
                                                                                                                                                                                                             aTTTTCAGCGTCAAAACCGTCGCCGTCGCCGCTTCGCTAACGCGTCACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGCCAGCCTGTGCCGACAAATCTGGCGTGCGGATGATAAATGGCATGTCATTTTTCGTGCATCCGGCTGGGAGCAACATATTACCACCCGCTATCTGGTCGGTGCCGATggagcca                                  >  2:312521/1‑245 (MQ=255)
                                                                                                                                                                                                                ttCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                       >  2:511828/1‑114 (MQ=255)
                                                                                                                                                                                                                ttCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                       <  1:511828/114‑1 (MQ=255)
                                                                                                                                                                                                                         aaaaCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTgcg                     >  1:510470/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACAtc             <  1:679570/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cATCGCTACCGCGTAACTACCAGCAAAGCTATAATAATA‑TGAACCTCAATGCTTTCGATTGGTGGATAAAATCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGat  <  2:173259/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGat  <  1:517929/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTC‑GATGTG‑ATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATA‑TTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGAT  >  NC_000913/3869396‑3869880

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