Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 659,205 A→T 5.3% intergenic (+55/+46) tatE → / ← lipA TatABCE protein translocation system subunit/lipoate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913659,2050AT 5.3% 311.0 / 4.5 104intergenic (+55/+46)tatE/lipATatABCE protein translocation system subunit/lipoate synthase
Reads supporting (aligned to +/- strand):  ref base A (49/48);  new base T (2/4);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 109 total aligned reads displayed.

AAGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGC  >  NC_000913/659056‑659352
                                                                                                                                                     |                                                                                                                                                   
aaGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTTGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCa                                                                                                                                                     >  2:1616509/1‑150 (MQ=255)
aaGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGCCGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCa                                                                                                                                                     >  2:620499/1‑150 (MQ=255)
aaGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCa                                                                                                                                                     >  2:567661/1‑150 (MQ=255)
    gTTCAAGAAGTCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGAAGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATAGCTTGTTCGATACAAAAAACCCGCTTCaaaaa                                                                                                                                                 >  1:296829/1‑150 (MQ=255)
       cAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCg                                                                                                                                              >  1:1524930/1‑150 (MQ=255)
          gaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGGGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGt                                                                                                                                           <  1:834051/150‑1 (MQ=255)
           aagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtt                                                                                                                                          >  1:309632/1‑150 (MQ=255)
            agGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttt                                                                                                                                         <  2:1408527/150‑1 (MQ=255)
            agGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttt                                                                                                                                         <  2:1412605/150‑1 (MQ=255)
              gCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttttt                                                                                                                                       <  1:535562/150‑1 (MQ=255)
              gCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttttt                                                                                                                                       >  1:1074433/1‑150 (MQ=255)
              gCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttttt                                                                                                                                       >  1:1080075/1‑150 (MQ=255)
               cGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTGGATATCTTGTTCGATACAAAAAACCAGCTTCAAAAAGCGGGtttttt                                                                                                                                      >  1:824760/1‑150 (MQ=255)
                 atgaatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttttttgt                                                                                                                                    <  2:1423427/150‑3 (MQ=255)
                   taatgaTGACGATGGTGCGGCGATAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAATGAGTGACGTTGCGAGCTGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAATGCGGGTTTTTTATca                                                                                                                                  <  2:1164763/149‑1 (MQ=255)
                   taatgaTGACGATGCTGCGGCGATAAAAGGCGCAGACGTTGATCTTAATGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGATCAGGACTCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATca                                                                                                                                  <  2:1164770/149‑1 (MQ=255)
                   gaatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATca                                                                                                                                  >  2:1024198/1‑150 (MQ=255)
                   gaatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATca                                                                                                                                  <  1:827203/150‑1 (MQ=255)
                    aatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcag                                                                                                                                 <  2:1217128/150‑1 (MQ=255)
                      tgatgaCGATGCTGCGTCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCGGCTTCAAAAAGCGGGTTTTTTATcagac                                                                                                                               <  1:1085688/150‑1 (MQ=255)
                            cGATGCTGCGGCGAAAAAAGGCGCAGCCGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgt                                                                                                                         >  1:833517/1‑150 (MQ=255)
                             gATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgta                                                                                                                        <  2:1520307/150‑1 (MQ=255)
                               tgctgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaag                                                                                                                      >  2:1509351/1‑150 (MQ=255)
                                  tgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagacag                                                                                                                             >  2:932772/1‑140 (MQ=255)
                                  tgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagacag                                                                                                                             <  1:932772/140‑1 (MQ=255)
                                  tacggcgAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagac                                                                                                                               >  2:971650/3‑138 (MQ=255)
                                  tacggcgAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagac                                                                                                                               <  1:971650/136‑1 (MQ=255)
                                             aaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGCGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTATTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                        >  1:67153/1‑150 (MQ=255)
                                                gcgcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTTACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTa                                                                                                     <  2:123830/150‑1 (MQ=255)
                                                gcgcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTa                                                                                                     >  1:1251627/1‑150 (MQ=255)
                                                 cgcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTAc                                                                                                    >  2:756802/1‑150 (MQ=255)
                                                  gcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACt                                                                                                   <  2:527526/150‑1 (MQ=255)
                                                  gcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACt                                                                                                   <  1:1024198/150‑1 (MQ=255)
                                                    aGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTa                                                                                                 <  2:146603/150‑1 (MQ=255)
                                                    aGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTa                                                                                                 <  1:1659021/150‑1 (MQ=255)
                                                     gACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCggg                                                                                                                                            <  1:337895/106‑1 (MQ=255)
                                                     gACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCggg                                                                                                                                            >  2:337895/1‑106 (MQ=255)
                                                      aCGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAAc                                                                                               <  1:811251/150‑1 (MQ=255)
                                                          tGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTcc                                                                                           >  1:1209235/1‑150 (MQ=255)
                                                           gATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCa                                                                                          >  1:696926/1‑150 (MQ=255)
                                                           gATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCa                                                                                          <  2:1524930/150‑1 (MQ=255)
                                                            aTCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                                                                         >  1:512914/1‑150 (MQ=255)
                                                              cTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATcc                                                                                       <  2:309632/150‑1 (MQ=255)
                                                                tCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCt                                                                                     >  1:1344504/1‑150 (MQ=255)
                                                                tCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCt                                                                                     <  1:1509351/150‑1 (MQ=255)
                                                                 cAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGAtt                                                                                                      >  2:281275/1‑132 (MQ=255)
                                                                   ggCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTc                                                                                  <  1:1188345/150‑1 (MQ=255)
                                                                    gCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCg                                                                                 <  1:259049/150‑1 (MQ=255)
                                                                    gCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCg                                                                                 <  2:891243/150‑1 (MQ=255)
                                                                     cTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGc                                                                                >  2:384330/1‑150 (MQ=255)
                                                                     cTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATAACAGGATTACTTAACTTCCATCCCTTTCGc                                                                                <  1:815396/150‑1 (MQ=255)
                                                                     cTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGAAGATAAATAAATATATTGTTCGATACAAAAAACCAGCTTAAAAAAGATTTTTTTTTATAAGACAGATGTAAGTAATTATTACAGGATTACTAAACAACCAAcccaaaccc                                                                                >  1:1459762/1‑144 (MQ=255)
                                                                      tGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGcc                                                                               <  2:691989/150‑1 (MQ=255)
                                                                       gtaaaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTACAAAAAGCGGGATTTTTATCAGACAGATGTAATTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCt                                                                              >  2:157177/3‑150 (MQ=255)
                                                                       gAAAAGCTCTCTCATAAAGAGTGACGTGGCGATCAGGACGCTCCCTCAATATCGTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATGATTACAGGATTACTTACCTTCCATCCCATTCGCCt                                                                              >  1:292863/1‑150 (MQ=255)
                                                                       gAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCt                                                                              <  1:1192817/150‑1 (MQ=255)
                                                                       gAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCACGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCCCTTCAAAAAGCGGGTTTTTTATAAGACAGATGTTAGTAATTATTACAGGATTACTTAACATCCATCCCTTTCGcat                                                                              >  2:1597262/1‑148 (MQ=255)
                                                                        aaaaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTg                                                                             >  1:1201925/1‑150 (MQ=255)
                                                                         aaaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGGTCCGTCAATATCTTGTTCGATACAAAAAAACCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTACTTCCCTTACGCCTGc                                                                            >  2:298409/1‑150 (MQ=255)
                                                                            gCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCaaa                                                                         >  2:21601/1‑150 (MQ=255)
                                                                             ctctctCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAAt                                                                        <  2:1209235/150‑1 (MQ=255)
                                                                               ctctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCg                                                                      >  1:1040987/1‑150 (MQ=255)
                                                                                tctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCgg                                                                     <  2:1344504/150‑1 (MQ=255)
                                                                                   aaTAAAGAGTGACGTGGCGAGCAGGACTCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGt                                                                  <  2:1295019/149‑1 (MQ=255)
                                                                                    aTAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                                                                 <  1:908429/150‑1 (MQ=255)
                                                                                     tAAAGAGTGACGTGGCGATCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTgg                                                                <  1:317864/150‑1 (MQ=255)
                                                                                     tAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTgg                                                                >  1:905492/1‑150 (MQ=255)
                                                                                       aaGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTa                                                              <  1:384330/150‑1 (MQ=255)
                                                                                         gagTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaag                                                            >  2:453070/1‑150 (MQ=255)
                                                                                             gACGTGTCGAGCAGTAAGCACCCACAATATCTAGTTCGAATCCAACAACCCGCTTCCAAAAGCGTGTTTTATATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACATCCATCCCTTTCGCCTGCAACTCGGCGTGGGAagaaga                                                        <  2:1617076/150‑1 (MQ=255)
                                                                                               cGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGc                                                      >  1:159093/1‑150 (MQ=255)
                                                                                                gTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAATAAGCGGGTTTTATATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCg                                                     >  1:591280/1‑150 (MQ=255)
                                                                                                    cGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACa                                                 <  2:1251627/150‑1 (MQ=255)
                                                                                                     taGCTGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAATAGCGGGTTTTTTTTCAGACTGTTGTAAGTAATTATTACAGGATTAATTAACTTCCAAACCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACaa                                                <  2:943151/149‑1 (MQ=255)
                                                                                                         aGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACgg                                            <  1:578426/150‑1 (MQ=255)
                                                                                                           gACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCTAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGAc                                          <  2:1074433/150‑1 (MQ=255)
                                                                                                           gACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGAc                                          <  2:833517/150‑1 (MQ=255)
                                                                                                           gACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGAc                                          <  2:1080075/150‑1 (MQ=255)
                                                                                                               cTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCa                                      >  2:1245604/1‑150 (MQ=255)
                                                                                                                 cccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGTCTGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGAGAAACGGACCGCATg                                    >  2:124091/1‑150 (MQ=255)
                                                                                                                  ccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTACATCCCTTTCGCCTGCAAATCCGCGTGGTAAGAAGAGAGGACAAACGGACCGCATgc                                   >  2:326947/1‑150 (MQ=255)
                                                                                                                      aaTATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATgcagca                               <  2:159093/150‑1 (MQ=255)
                                                                                                                        tatCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATg                             <  1:21601/150‑1 (MQ=255)
                                                                                                                        tatCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATg                             <  2:901219/150‑1 (MQ=255)
                                                                                                                          tCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATggg                           >  2:961332/1‑150 (MQ=255)
                                                                                                                                   tgataaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCa                  >  1:1422478/5‑150 (MQ=255)
                                                                                                                                   aTACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAATAAGAGCGGACAAACGGAACGCATGCAGAATGGGTAAAGccaa                  >  2:376083/1‑148 (MQ=255)
                                                                                                                                       aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGc              <  1:567661/150‑1 (MQ=255)
                                                                                                                                       aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGc              <  1:620499/150‑1 (MQ=255)
                                                                                                                                       aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGc              >  2:807286/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGcc             >  2:562852/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCa            >  1:1276223/1‑150 (MQ=255)
                                                                                                                                          aaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAg           >  1:1544003/1‑150 (MQ=255)
                                                                                                                                            aCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAgcg         >  1:11208/1‑150 (MQ=255)
                                                                                                                                                gCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCAATCGCCAGCGCTTc     >  2:1306936/1‑150 (MQ=255)
                                                                                                                                                 cTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCg    <  1:453070/150‑1 (MQ=255)
                                                                                                                                                  ttCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAAACGGCGTGGTAATAAGAGCGGACAAAAGGACCGCATTCAGCAAGGGTAAAGCCAATCGCCAGCTCTTCgg   >  2:162821/1‑150 (MQ=255)
                                                                                                                                                   tCAAAAATCGGGTTTTTTATCAGACAGATGTAAGTATTTATTTGAGGATTACTTAACTTCCCTCCCTTTCGCCTGCAAATCGGCGTGGTACGCAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGc  <  1:162821/150‑1 (MQ=255)
                                                                                                                                                   tCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGc  <  1:617839/150‑1 (MQ=255)
                                                                                                                                                   tCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGc  >  2:1493340/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGC  >  NC_000913/659056‑659352

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