Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 15,338 A→T intergenic (+40/‑107) 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,3380AT93.8% 182.8 / ‑4.7 66intergenic (+40/‑107)dnaJ/insL1chaperone protein DnaJ/IS186/IS421 family transposase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (18/43);  minor base G (0/4);  total (19/47)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTC  >  NC_000913/15092‑15580
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
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)
                                                                       ccTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTAAACGGCAAACTTAAGGGttt                                                                                                                                                                          >  1:349855/1‑250 (MQ=255)
                                                                       ccTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgt                                                                                                                                                                         <  1:175811/251‑1 (MQ=255)
                                                                        cTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAATCCCCAAAACGCTAACTTTAAGGGtgtg                                                                                                                                                                         >  1:450319/1‑248 (MQ=255)
                                                                          gAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCGCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGt                                                                                                                                                                      <  2:47701/251‑1 (MQ=255)
                                                                             cgagaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                                                                                                                                                   <  1:152779/245‑1 (MQ=255)
                                                                             cgagaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCAAAAAGCCTGCCCGTGGGCAGGCTTGGGTAACATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCTCTCCCATCAACTCTAACATAAGGGGTGTgggttt                                                                                                                                                                   >  2:498974/7‑247 (MQ=255)
                                                                               agaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                                                                                 <  1:159834/247‑1 (MQ=255)
                                                                                gaagCAGAAACATCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTGGACGACCTGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATAGGCGAGCACCTGAAAAGTGGCGGGGATCACGCCCATAAGCGCTAACTGAAGGGTAGTGGTATTACg                                                                                                                                                                <  2:422747/248‑1 (MQ=255)
                                                                                  agCAGAAACAGCTGCTGCCAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGACGATATGCGAGCACCTGAAAAGTGGCGGGGGTACCACCCCTAAGCGCTTACTTTAGGGTTGGGGTATTAAGcc                                                                                                                                                              >  1:189465/2‑251 (MQ=255)
                                                                                     aGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                                                                                                                                           <  1:10040/251‑1 (MQ=255)
                                                                                       aaaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACACCCCTAAACCCCAACATAAAGGTTGTggggttacaccccgt                                                                                                                                                          >  1:380401/1‑237 (MQ=255)
                                                                                        aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                                                                                                                                                        <  2:198576/251‑1 (MQ=255)
                                                                                        aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGACACTCCCACAAACGCAAACTTACAGGTTGGGGTAtttatccatgta                                                                                                                                                         >  1:454305/1‑240 (MQ=255)
                                                                                        aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                                                                                                                                                        <  1:120610/251‑1 (MQ=255)
                                                                                            gctgctGCAAGAGCTGCAAGAAAGCTGCGGTTGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGATCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCGTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGTGTTGTGTTATTACGCCTGATATGAtt                                                                                                                                                    <  2:555209/251‑1 (MQ=255)
                                                                                               gctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCTCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTGGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                                   <  1:586745/249‑1 (MQ=255)
                                                                                               gctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                                   <  2:13169/249‑1 (MQ=255)
                                                                                                ctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGGGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTGACACCCCTAAACGCTAACTTAAGGGTTGTGGTTTTTCGACCCATATGATTTatt                                                                                                                                                >  1:125552/1‑249 (MQ=255)
                                                                                                   caagagctgcaagaAAGCTGCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCGTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCGGACCCGCTAACCTCCCCAAAAGCCGGCCCGTGGGCAGGCTTGGGTAAAAGTAGGGTGCGTTGAAGATAGGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAGAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                                                                                                                                             <  1:330087/251‑1 (MQ=255)
                                                                                                      gagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGGAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                                                                          <  1:566777/251‑1 (MQ=255)
                                                                                                           gcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                                                                       <  2:450319/249‑1 (MQ=255)
                                                                                                                  aGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct                                                                                                                              <  2:188528/251‑1 (MQ=255)
                                                                                                                   gCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                                                                             <  2:503815/251‑1 (MQ=255)
                                                                                                                      tCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTTGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                                                                                           <  1:505929/250‑1 (MQ=255)
                                                                                                                       cGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACGCTCACg                                                                                                                          <  2:4021/250‑1 (MQ=255)
                                                                                                                             cccAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                                                                   <  1:249522/251‑1 (MQ=255)
                                                                                                                                  ccttggttcacaACCGCCCGCGCTCAAAGAGCCTCTGTGATGGTGTGAAGAAGTTTTTTGACGACCTCACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCACAAGCGCTCACTTAACGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                                                                              <  2:125552/243‑1 (MQ=37)
                                                                                                                                  ccGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACCTCCATTACCGCCAACTTAAAGGTTGTGGGATTTCGCCCGATTTTATTTTACGTGCCGATGGATTTCTTTCACATTAACTgttc                                                                                                                >  2:14455/1‑246 (MQ=255)
                                                                                                                                        aGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTc                                                                                                        <  1:624504/251‑1 (MQ=255)
                                                                                                                                           acaacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                                                                     <  1:119982/251‑1 (MQ=255)
                                                                                                                                            caacaGCCCGCGCTAAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAAGTAGGGTGCGTTGAAGATATGCGAGCTCCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                                                                                    <  2:4485/251‑1 (MQ=255)
                                                                                                                                                      cgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCATAAGCGCTAAATTAAAGGGTGTGGGATTACGCCTGATATGATTTTAAGTGCCGGTGAATTATTCTCACGATAACTGGTCAGCAAATCTGGCCCCTATTGGt                                                                                          >  1:244270/1‑251 (MQ=255)
                                                                                                                                                                      ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCTCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACTTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTg                                                                           <  1:71049/250‑1 (MQ=255)
                                                                                                                                                                                 aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCa                                                               <  2:523059/251‑1 (MQ=255)
                                                                                                                                                                                 aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGCTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTAAATCCCAAAAATGCCAACTTAAAGGCTTGGGGGTTACCCCTGCTATGATTTTACGCTCCCAAGAAATATTCTCACGATAACAGGGTACTAATTCTTGCCCAAATTGGTAAGCACgtagaaataggtaattctgcc                                                               >  1:95824/1‑236 (MQ=255)
                                                                                                                                                                                        tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt                                                         <  1:603409/250‑1 (MQ=255)
                                                                                                                                                                                                     gACCCTCTAACCTCACCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTGACTTAAGGGTTGTGGTATTACGCCTGATATGAGTTAGCGTGCCGATGAATTTCTCTCACGATAACTGGTCAGCAATTCTGGCCCATAGTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTa                                             <  1:620925/249‑1 (MQ=18)
                                                                                                                                                                                                                          gCCTGCCCGTGTGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGTGATCACTCCCATAAGCGCTAACTTAAGGGTTGGGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGa                      <  1:611921/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGCAAAGTGGCGGGGATAACTCCCATAAACGCTCACTTTAGGGTTGGGGGATTAAGCCCGATATGATTTAAAGTGCCGCTGAATTAATCTCACGATAACTGGTCAGCAATTCTGGCCCATATTTGTAAACCCGCAGAACTGGATAATTCCGCACCTAATGACGGGGCTCTAACCCGCCGCCGCGAAAATCGGGAtgctgccac             >  1:200007/1‑248 (MQ=17)
                                                                                                                                                                                                                                   ggggCAGTCTTGGGTAAAAGTAGGGTGCGTTGAAGATATGCTAGCATCTGAAAAGTGGCGTGGATCACGCCCATAAGCGCGAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAAc             <  2:511982/250‑1 (MQ=14)
                                                                                                                                                                                                                                   ggggCAGGCGTGGGGAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACTTGCCGATGAATGACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCTTGATGCTGCAAc             <  2:301891/250‑1 (MQ=17)
                                                                                                                                                                                                                                              gggTAAAATTAGGGGGCGTGGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTc  <  2:349855/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                                                                                                                                            <  1:95999/111‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTC  >  NC_000913/15092‑15580

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: