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,3280CT97.6% 124.3 / NA 41intergenic (+30/‑117)dnaJ/insL1chaperone protein DnaJ/IS186/IS421 family transposase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/32);  minor base A (0/1);  total (8/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATG‑CCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGC  >  NC_000913/15091‑15574
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
aGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAACGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTaaaatt                                                                                                                                                                                                                                                >  1:48785/1‑247 (MQ=255)
   cGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAgggg                                                                                                                                                                                                                                           >  1:412761/1‑250 (MQ=255)
      cAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGc                                                                                                                                                                                                                                         <  1:72484/250‑1 (MQ=255)
        aggcTGTCAGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGATCTTCTTTTATGGTGTGAAGAAGTTTTTTGACGTCCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGTCTTGGGTAAAATTAGGGTGCGtt                                                                                                                                                                                                                                      <  2:430518/248‑1 (MQ=255)
           cTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCACGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGa                                                                                                                                                                                                                                    >  1:430310/1‑250 (MQ=255)
                    ggggCGCACAGGGTGATTTGCGGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGa                                                                                                                                                                                                                          <  2:22421/249‑1 (MQ=255)
                                gTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGCTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTgaaa                                                                                                                                                                                                                >  1:52118/1‑246 (MQ=255)
                                gTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTgaaa                                                                                                                                                                                                                >  2:175414/1‑246 (MQ=255)
                                          gtgCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGAt                                                                                                                                                                                                    <  2:106283/251‑1 (MQ=255)
                                            gCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCa                                                                                                                                                                                                  <  1:175414/251‑1 (MQ=255)
                                                   tGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCCAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACGccccc                                                                                                                                                                                           >  1:58386/1‑249 (MQ=255)
                                                        gcgAAACCCCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCTCTCCCATAAgcg                                                                                                                                                                                      <  2:159008/250‑1 (MQ=255)
                                                           aaaCACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTa                                                                                                                                                                                   <  1:436699/251‑1 (MQ=255)
                                                             acacCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAAc                                                                                                                                                                                 <  2:58386/251‑1 (MQ=255)
                                                                                   agCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGATCCGCTAACCACCCCAAAAGCCTGCCCGGGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGcc                                                                                                                                                          <  1:202386/250‑1 (MQ=255)
                                                                                      aGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCTCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGa                                                                                                                                                       <  1:363802/251‑1 (MQ=255)
                                                                                       gAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCGTTGAGGGTGTGAAGAAGTTTTTTGACGACCTGCCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGat                                                                                                                                                      <  1:41993/251‑1 (MQ=255)
                                                                                          aCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATg                                                                                                                                                   <  1:124476/251‑1 (MQ=255)
                                                                                                       gagctgcaagaAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                                                                      <  1:122538/251‑1 (MQ=255)
                                                                                                                  aaGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTtg                                                                                                                                                                      <  2:375799/208‑1 (MQ=255)
                                                                                                                   aGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                               <  2:409402/230‑1 (MQ=255)
                                                                                                                   aGCTTCGGGGGCCCAGCCGGTGAGCACAACAGCCCGCTCCTAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct                                                                                                                          <  2:179613/251‑1 (MQ=38)
                                                                                                                          gTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa                                                                                                                   <  1:48435/251‑1 (MQ=255)
                                                                                                                            ggCCCAACCGGTGATCATTACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTGTTTTGACGACCTGACCCGTTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCATGGGTAAAATGAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCTCTCCCATAAGCTCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCCCGATAAc                                                                                                                 <  2:363208/251‑1 (MQ=37)
                                                                                                                                           cacaACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTGTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGTGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGACAACTGGTC‑AGCAATTCTg                                                                                                 <  2:280715/251‑1 (MQ=255)
                                                                                                                                              aacaGCCCGCGCGCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGTCCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGcc                                                                                              <  2:81234/251‑1 (MQ=255)
                                                                                                                                               aaaGCCCGCGGTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGccc                                                                                             <  2:338520/249‑1 (MQ=255)
                                                                                                                                                 aGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCat                                                                                           <  1:235756/251‑1 (MQ=255)
                                                                                                                                                      gcgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTgg                                                                                      <  1:263975/251‑1 (MQ=255)
                                                                                                                                                             aaGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGccc                                                                               <  1:265387/251‑1 (MQ=255)
                                                                                                                                                                         gATGGTGTGAAGAAGTTTTTTGACGAACTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                                                                   <  2:435859/251‑1 (MQ=255)
                                                                                                                                                                               gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCg                                                             <  2:88269/251‑1 (MQ=255)
                                                                                                                                                                                  aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCa                                                          <  2:403648/251‑1 (MQ=255)
                                                                                                                                                                                   agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc                                                         <  1:504419/251‑1 (MQ=255)
                                                                                                                                                                                                     tGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTGACTCCCCTAAACGCTAACTTTAAGGGT‑GTGGGGTTACGCCTGGTATGATTTTACGCTCCGCTTAATTACTTTCCCGATAAATCGTCCACCAAT‑CTTGCCCACATTGTTAAGCCCGACGAAATGGATATCTCTGCCCGTAATGACCGGGGcctt                                        >  2:115043/1‑247 (MQ=255)
                                                                                                                                                                                                      gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATG‑CCGGGGCTCTAAc                                     <  1:54307/251‑1 (MQ=255)
                                                                                                                                                                                                             tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGACACTCCCCTAAACGCTAACTT‑AAAGGTTGTGGGGTTAAGCCTGGTATGATTTTAAATGCCGGTGAATTACTCTCACGATAAATGGTC‑AGCAAATCTTGCCCATATTTGTAAACCCGAAGAACTgg                                                                     >  1:437557/1‑213 (MQ=255)
                                                                                                                                                                                                             tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                                                                     <  2:437557/213‑1 (MQ=255)
                                                                                                                                                                                                                     ccAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATG‑CCGGGGCTCTAACCCGCCGCCGCGAAAt                      <  1:100641/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCTGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATG‑CCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCaa        <  2:128810/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aGGCTTGGGTAAAATGAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATG‑CCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGc  <  1:201439/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
AGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTC‑AGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATG‑CCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGC  >  NC_000913/15091‑15574

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