Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 659,205 A→T 7.4% 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 7.4% 324.7 / 12.5 95intergenic (+55/+46)tatE/lipATatABCE protein translocation system subunit/lipoate synthase
Reads supporting (aligned to +/- strand):  ref base A (41/47);  new base T (2/5);  total (43/52)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

AGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGCT  >  NC_000913/659057‑659353
                                                                                                                                                    |                                                                                                                                                    
aGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCa                                                                                                                                                      >  2:826862/1‑149 (MQ=255)
   gTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaaaa                                                                                                                                                   >  1:826743/1‑149 (MQ=255)
      cAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGc                                                                                                                                                <  4:355959/149‑1 (MQ=255)
      cAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGc                                                                                                                                                <  2:826743/149‑1 (MQ=255)
      cAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGc                                                                                                                                                <  1:473918/149‑1 (MQ=255)
        agaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCgg                                                                                                                                              <  4:681737/149‑1 (MQ=255)
           agGCGATGAATGATGACGATGATGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGAACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGtg                                                                                                                                           <  4:200805/149‑2 (MQ=255)
            ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttt                                                                                                                                          <  4:741438/149‑1 (MQ=255)
            ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttt                                                                                                                                          <  3:706034/149‑1 (MQ=255)
            ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttt                                                                                                                                          >  1:891708/1‑149 (MQ=255)
             gCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGtttt                                                                                                                                         <  2:814843/149‑1 (MQ=255)
              cGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttttt                                                                                                                                        <  4:864001/149‑1 (MQ=255)
                  gaatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATc                                                                                                                                    <  2:247809/149‑1 (MQ=255)
                   aatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATca                                                                                                                                   >  3:246291/1‑149 (MQ=255)
                    atgatgACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcag                                                                                                                                  <  4:332341/149‑1 (MQ=255)
                       atgaCGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagaca                                                                                                                               <  3:577181/149‑1 (MQ=255)
                       atgaCGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCCTAAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCGGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagaca                                                                                                                               >  2:699358/1‑149 (MQ=255)
                              tgctgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaa                                                                                                                        <  1:296245/149‑1 (MQ=255)
                                ctgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagt                                                                                                                      >  4:516396/1‑149 (MQ=255)
                                 tgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagta                                                                                                                     >  2:54602/1‑149 (MQ=255)
                                 tgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagta                                                                                                                     >  2:273904/1‑149 (MQ=255)
                                       gAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagtaa                                                                                                                    <  4:621528/144‑1 (MQ=255)
                                       gAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagtaa                                                                                                                    >  3:621528/1‑144 (MQ=255)
                                        aaaaaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAAttatta                                                                                                              >  1:805016/1‑149 (MQ=255)
                                         aaaaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTAc                                                                                                             <  4:381944/149‑1 (MQ=255)
                                           aaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAg                                                                                                           >  4:375790/1‑149 (MQ=255)
                                             aGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                         >  2:190161/1‑149 (MQ=255)
                                             aGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                         <  1:310853/149‑1 (MQ=255)
                                             aGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                         <  2:348184/149‑1 (MQ=255)
                                              ggCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGAt                                                                                                        <  2:19438/149‑1 (MQ=255)
                                               gcgcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGCCAGATGTAAGTAATTATTACAGGAtt                                                                                                       >  1:106337/1‑149 (MQ=255)
                                                 gcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTAc                                                                                                     <  1:102635/149‑1 (MQ=255)
                                                     aCGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTaa                                                                                                 <  3:849217/149‑1 (MQ=255)
                                                        ttGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACtt                                                                                              >  4:855858/1‑149 (MQ=255)
                                                         tGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTa                                                                                                      >  3:320367/1‑140 (MQ=255)
                                                         tGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTa                                                                                                      <  4:320367/140‑1 (MQ=255)
                                                          gATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTcc                                                                                            >  1:25767/1‑149 (MQ=255)
                                                          gATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTcc                                                                                            <  3:258680/149‑1 (MQ=255)
                                                           aTCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCa                                                                                           >  2:752602/1‑149 (MQ=255)
                                                           aTCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCa                                                                                           <  4:55230/149‑1 (MQ=255)
                                                            tCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                                                                          >  1:797309/1‑149 (MQ=255)
                                                              ttCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATcc                                                                                        >  2:817094/1‑149 (MQ=255)
                                                               tCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATccc                                                                                       <  4:30783/149‑1 (MQ=255)
                                                                cAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCt                                                                                      <  2:302960/149‑1 (MQ=255)
                                                                 aGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCtt                                                                                     <  3:524310/149‑1 (MQ=255)
                                                                  ggCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCttt                                                                                    >  3:78609/1‑149 (MQ=255)
                                                                   gCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTc                                                                                   <  2:805016/149‑1 (MQ=255)
                                                                   gCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTc                                                                                   <  4:555796/149‑1 (MQ=255)
                                                                   gCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTc                                                                                   <  4:656667/149‑1 (MQ=255)
                                                                   gCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTc                                                                                   <  4:203697/149‑1 (MQ=255)
                                                                     tGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                                                                          >  3:140553/1‑140 (MQ=255)
                                                                     tGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                                                                          <  4:140553/140‑1 (MQ=255)
                                                                       aaaaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCt                                                                               <  4:78609/149‑1 (MQ=255)
                                                                             tctctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAAt                                                                         >  2:219527/1‑149 (MQ=255)
                                                                              ctctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATc                                                                        <  4:518341/149‑1 (MQ=255)
                                                                              ctctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATc                                                                        <  1:283339/149‑1 (MQ=255)
                                                                               tctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCg                                                                       >  1:797379/1‑149 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTAc                                                                                                             <  3:266205/109‑1 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTAc                                                                                                             >  4:266205/1‑109 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTa                                                                                                  <  1:557629/120‑1 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTa                                                                                                  >  2:557629/1‑120 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                                                                          <  3:80775/128‑1 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                                                                          >  4:80775/1‑128 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGc                                                                     <  1:826862/149‑1 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGc                                                                     >  1:382601/1‑149 (MQ=255)
                                                                                 tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGc                                                                     <  3:312324/149‑1 (MQ=255)
                                                                                         agTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaag                                                             <  3:627825/149‑1 (MQ=255)
                                                                                          gTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaa                                                              <  1:302655/147‑1 (MQ=255)
                                                                                          gTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaa                                                              >  2:302655/1‑147 (MQ=255)
                                                                                               gTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGc                                                       <  1:173094/149‑1 (MQ=255)
                                                                                                tGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCg                                                      <  3:837813/149‑1 (MQ=255)
                                                                                                  gCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGa                                                    <  4:23456/149‑1 (MQ=255)
                                                                                                  gCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGa                                                    <  2:891708/149‑1 (MQ=255)
                                                                                                  gCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGAAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGa                                                    >  4:595832/1‑149 (MQ=255)
                                                                                                    gAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACa                                                  >  3:314424/1‑149 (MQ=255)
                                                                                                    gAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACa                                                  >  1:789723/1‑149 (MQ=255)
                                                                                                      gCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACaaa                                                <  4:256047/149‑1 (MQ=255)
                                                                                                      gCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACaaa                                                <  1:631641/149‑1 (MQ=255)
                                                                                                       ccggACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAACTCGGCGTGGTAAGAAGAGCGGACAAAc                                               >  2:790296/3‑149 (MQ=255)
                                                                                                       cAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAAc                                               <  4:589139/149‑1 (MQ=255)
                                                                                                          gACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGa                                            >  2:635595/1‑149 (MQ=255)
                                                                                                           aCGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGAc                                           <  4:78974/149‑1 (MQ=255)
                                                                                                                cccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCAt                                      >  3:129835/1‑149 (MQ=255)
                                                                                                                 ccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATg                                     >  3:682246/1‑149 (MQ=255)
                                                                                                                        atCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATg                              <  4:246291/149‑1 (MQ=255)
                                                                                                                         tCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATgg                             >  1:557275/1‑149 (MQ=255)
                                                                                                                         tCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATgg                             >  2:879981/1‑149 (MQ=255)
                                                                                                                           ttGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGt                           >  2:571351/1‑149 (MQ=255)
                                                                                                                                   gataaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCa                   >  2:68861/4‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCg                                                      >  4:715180/1‑111 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCg                                                      <  3:715180/111‑1 (MQ=255)
                                                                                                                                         aaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGATAAAGCCCATCGCCa             >  3:184444/1‑149 (MQ=255)
                                                                                                                                         aaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCa             >  3:610619/1‑149 (MQ=255)
                                                                                                                                          aaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAg            >  4:695952/1‑149 (MQ=255)
                                                                                                                                             ccGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAgcgc         >  4:530266/1‑149 (MQ=255)
                                                                                                                                             ccGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAgcgc         >  3:621695/1‑149 (MQ=255)
                                                                                                                                                   cAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGc   <  1:273904/149‑1 (MQ=255)
                                                                                                                                                    aaaaaGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGCt  <  1:235720/149‑1 (MQ=255)
                                                                                                                                                    aaaaaGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGCt  >  1:160522/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGCT  >  NC_000913/659057‑659353

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