Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 15,328 C→T intergenic (+30/‑117) dnaJ → / → insL1 chaperone protein DnaJ/IS186/IS421 family transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,3280CT95.5% 196.5 / NA 66intergenic (+30/‑117)dnaJ/insL1chaperone protein DnaJ/IS186/IS421 family transposase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/44);  minor base G (0/3);  total (19/47)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAAC  >  NC_000913/15083‑15569
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cgcgGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCGGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCTCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGt                                                                                                                                                                                                                                                 <  1:70834/251‑1 (MQ=255)
       aaGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTaaaatta                                                                                                                                                                                                                                          <  2:4330/251‑4 (MQ=255)
         ggCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCGGAACGAGAAGCAGAAACAGGTGCGGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGGCCTGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCGGGCGGGGGTAAAATTAgg                                                                                                                                                                                                                                        <  2:259091/251‑1 (MQ=255)
          gCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAggg                                                                                                                                                                                                                                       <  1:422946/251‑1 (MQ=255)
          gCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTGTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAggg                                                                                                                                                                                                                                       <  1:568102/251‑1 (MQ=255)
            gTCAAGTCTGTCCGCGGTGGCGCACAGGGTGCTTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAACAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTg                                                                                                                                                                                                                                     >  2:563507/1‑251 (MQ=255)
            gTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTCTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTg                                                                                                                                                                                                                                     >  1:396251/1‑251 (MQ=255)
            gTCAAGTCTGTACGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACTGCCCGCGCTCAAAGAGCTTCTTTGAATGTGTGAAGAAGTGTGGTGACGACCTGGCCCGCTAACCGCGCCAAAAGGGTGCCCGTGGGCAGTCTTGGGTAAAAGTAGGGTg                                                                                                                                                                                                                                     <  2:102665/251‑1 (MQ=255)
               aaGTCTGTCCGCGGTGGCGCACAGGGTGATTTACTGTACCGCGTTGTCATCGCAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCATCATTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGACGTTTTTTGACGACCGGACCCGCTAACCTCCCCAACAGCCTGCACGTGGGCAGGATTGGGTAATACTAGGGAGCGt                                                                                                                                                                                                                                  >  2:604094/1‑251 (MQ=255)
               aaGTCTGTCCGCGGAGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAACCACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCACAAGCCTGCCCGTGGGCCGGCTTGGGTAAAATTAGGGTGCGt                                                                                                                                                                                                                                  >  2:604097/1‑251 (MQ=255)
                aGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTCTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGtt                                                                                                                                                                                                                                 <  2:396251/251‑1 (MQ=255)
                aGTCTGTCCGCGGGGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGtt                                                                                                                                                                                                                                 <  2:147750/251‑1 (MQ=255)
                   cggTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTTGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTTGGCAGGCTTGGGTAAAATTAGGGTGCGTTGaa                                                                                                                                                                                                                              <  2:397683/249‑1 (MQ=255)
                       ccGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGata                                                                                                                                                                                                                          <  2:459200/251‑1 (MQ=255)
                       ccGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGata                                                                                                                                                                                                                          <  1:438340/251‑1 (MQ=255)
                          cgGTTGCGCACAGGTGGATTGGCTGTGCCGCGTTGTCGTCGAAACACCGGGAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGGAAAATTAGGGTGCGTTGAAGATATGc                                                                                                                                                                                                                       <  2:604744/251‑1 (MQ=255)
                                     aGGGTGATGTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTgaa                                                                                                                                                                                                            <  1:148046/251‑3 (MQ=255)
                                               gCTGTGCCGCGTGGTCGTCGAAACCCCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATGAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCggg                                                                                                                                                                                                  <  2:456894/251‑1 (MQ=255)
                                                    gCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAGAGTGGCGGGGATc                                                                                                                                                                                              <  1:236502/250‑1 (MQ=255)
                                                     ccGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCTCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCgggggtacc                                                                                                                                                                                            >  2:112392/1‑246 (MQ=255)
                                                       gcgTTGTCGTCGAAACACCGGTAGGCCGGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTc                                                                                                                                                                                          <  2:28400/251‑1 (MQ=255)
                                                          ttGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAATccc                                                                                                                                                                                        >  1:364168/1‑250 (MQ=255)
                                                            gtcgtcGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                                                                                                                                                                      <  2:237356/250‑1 (MQ=255)
                                                                 cgAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACcccccttagagc                                                                                                                                                                                >  2:428753/1‑243 (MQ=255)
                                                                   aaaCACCGGTAGGCCTGAACGAGCAGCAGAAACAGCTGCTGCAAAAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCggaggtcaatccccaaacagcc                                                                                                                                                                               >  2:212975/1‑230 (MQ=255)
                                                                                ccTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTAAACGGCAAACTT‑AA‑GGGttt                                                                                                                                                               >  1:349855/1‑250 (MQ=255)
                                                                                ccTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTtgt                                                                                                                                                              <  1:175811/251‑1 (MQ=255)
                                                                                 cTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAATCCCCAAAACG‑CTAACTTTAA‑GGGtgtg                                                                                                                                                              >  1:450319/1‑248 (MQ=255)
                                                                                   gAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCGCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGt                                                                                                                                                           <  2:47701/251‑1 (MQ=255)
                                                                                      cgagaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTAtt                                                                                                                                                        <  1:152779/245‑1 (MQ=255)
                                                                                      cgagaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCAAAAAGCCTGCCCGTGGGCAGGCTTGGGTAACATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCTCTCCCATCAACT‑CTAACAT‑AA‑GGGGTGTgggttt                                                                                                                                                        >  2:498974/7‑247 (MQ=255)
                                                                                        agaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTAc                                                                                                                                                      <  1:159834/247‑1 (MQ=255)
                                                                                         gaagCAGAAACATCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTGGACGACCTGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATAGGCGAGCACCTGAAAAGTGGCGGGGATCACGCCCATAAGCG‑CTAACTG‑AA‑GGGTAGTGGTATTACg                                                                                                                                                     <  2:422747/248‑1 (MQ=255)
                                                                                           agCAGAAACAGCTGCTGCCAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGACGATATGCGAGCACCTGAAAAGTGGCGGGGGTACCACCCCTAAGCG‑CTTACTT‑TA‑GGGTTGGGGTATTAAGcc                                                                                                                                                   >  1:189465/2‑251 (MQ=255)
                                                                                              aGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGa                                                                                                                                                <  1:10040/251‑1 (MQ=255)
                                                                                                aaaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACACCCCTAAACC‑CCAACAT‑AA‑AGGTTGTggggttacaccccgt                                                                                                                                               >  1:380401/1‑237 (MQ=255)
                                                                                                 aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGatat                                                                                                                                             <  2:198576/251‑1 (MQ=255)
                                                                                                 aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGACACTCCCACAAACG‑CAAACTT‑AC‑AGGTTGGGGTAtttatccatgta                                                                                                                                              >  1:454305/1‑240 (MQ=255)
                                                                                                 aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGatat                                                                                                                                             <  1:120610/251‑1 (MQ=255)
                                                                                                     gctgctGCAAGAGCTGCAAGAAAGCTGCGGTTGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGATCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCGTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GTGTTGTGTTATTACGCCTGATATGAtt                                                                                                                                         <  2:555209/251‑1 (MQ=255)
                                                                                                        gctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCTCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTGGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                        <  1:586745/249‑1 (MQ=255)
                                                                                                        gctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                        <  2:13169/249‑1 (MQ=255)
                                                                                                         ctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGGGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTGACACCCCTAAACG‑CTAACTT‑AA‑GGGTTGTGGTTTTTCGACCCATATGATTTatt                                                                                                                                     >  1:125552/1‑249 (MQ=255)
                                                                                                            caagagctgcaagaAAGCTGCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCGTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCGGACCCGCTAACCTCCCCAAAAGCCGGCCCGTGGGCAGGCTTGGGTAAAAGTAGGGTGCGTTGAAGATAGGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAGAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTg                                                                                                                                  <  1:330087/251‑1 (MQ=255)
                                                                                                               gagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGGAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                                                               <  1:566777/251‑1 (MQ=255)
                                                                                                                    gcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                                                            <  2:450319/249‑1 (MQ=255)
                                                                                                                           aGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct                                                                                                                   <  2:188528/251‑1 (MQ=255)
                                                                                                                            gCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                                                                  <  2:503815/251‑1 (MQ=255)
                                                                                                                               tCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTTGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                                                                                <  1:505929/250‑1 (MQ=255)
                                                                                                                                cGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACGCTCACg                                                                                                               <  2:4021/250‑1 (MQ=255)
                                                                                                                                      cccAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                                                        <  1:249522/251‑1 (MQ=255)
                                                                                                                                           ccttggttcacaACCGCCCGCGCTCAAAGAGCCTCTGTGATGGTGTGAAGAAGTTTTTTGACGACCTCACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCACAAGCG‑CTCACTT‑AA‑CGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                                                                   <  2:125552/243‑1 (MQ=37)
                                                                                                                                           ccGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACCTCCATTACCG‑CCAACTT‑AA‑AGGTTGTGGGATTTCGCCCGATTTTATTTTACGTGCCGATGGATTTCTTTCACATTAACTgttc                                                                                                     >  2:14455/1‑246 (MQ=255)
                                                                                                                                                 aGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTc                                                                                             <  1:624504/251‑1 (MQ=255)
                                                                                                                                                    acaacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                                                          <  1:119982/251‑1 (MQ=255)
                                                                                                                                                     caacaGCCCGCGCTAAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAAGTAGGGTGCGTTGAAGATATGCGAGCTCCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                                                                         <  2:4485/251‑1 (MQ=255)
                                                                                                                                                               cgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCATAAGCG‑CTAAATT‑AA‑AGGGTGTGGGATTACGCCTGATATGATTTTAAGTGCCGGTGAATTATTCTCACGATAACTGGTCAGCAAATCTGGCCCCTATTGGt                                                                               >  1:244270/1‑251 (MQ=255)
                                                                                                                                                                               ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCT‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACTTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTg                                                                <  1:71049/250‑1 (MQ=255)
                                                                                                                                                                                          aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCa                                                    <  2:523059/251‑1 (MQ=255)
                                                                                                                                                                                          aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGCTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTAAATCCCAAAAATG‑CCAACTT‑AAAGGCTTG‑GGGGTTACCCCTGCTATGATTTTACGCTCCCAAGAAATATTCTCACGATAACAGGGTACTAATTCTTGCCCAAATTGGTAAGCACgtagaaataggtaattctgcc                                                    >  1:95824/1‑236 (MQ=255)
                                                                                                                                                                                                 tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt                                              <  1:603409/250‑1 (MQ=255)
                                                                                                                                                                                                              gACCCTCTAACCTCACCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTGACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGAGTTAGCGTGCCGATGAATTTCTCTCACGATAACTGGTCAGCAATTCTGGCCCATAGTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTa                                  <  1:620925/249‑1 (MQ=18)
                                                                                                                                                                                                                                   gCCTGCCCGTGTGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGTGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGGGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGa           <  1:611921/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGCAAAGTGGCGGGGATAACTCCCATAAACG‑CTCACTT‑TA‑GGGTTGGGGGATTAAGCCCGATATGATTTAAAGTGCCGCTGAATTAATCTCACGATAACTGGTCAGCAATTCTGGCCCATATTTGTAAACCCGCAGAACTGGATAATTCCGCACCTAATGACGGGGCTCTAACCCGCCGCCGCGAAAATCGGGAtgctgccac  >  1:200007/1‑248 (MQ=17)
                                                                                                                                                                                                                                            ggggCAGTCTTGGGTAAAAGTAGGGTGCGTTGAAGATATGCTAGCATCTGAAAAGTGGCGTGGATCACGCCCATAAGCG‑CGAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAAc  <  2:511982/250‑1 (MQ=14)
                                                                                                                                                                                                                                            ggggCAGGCGTGGGGAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACTTGCCGATGAATGACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCTTGATGCTGCAAc  <  2:301891/250‑1 (MQ=17)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAAC  >  NC_000913/15083‑15569

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