Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 659,205 A→T 17.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,2050AT17.4% 414.2 / 78.7 180intergenic (+55/+46)tatE/lipATatABCE protein translocation system subunit/lipoate synthase
Reads supporting (aligned to +/- strand):  ref base A (82/65);  new base T (23/8);  total (106/74)
Fisher's exact test for biased strand distribution p-value = 7.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 185 total aligned reads displayed.

GCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGAC  >  NC_000913/659106‑659304
                                                                                                   |                                                                                                   
 naGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaa                                                                                                    >  1:82357/2‑100 (MQ=255)
       ttGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGACGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCg                                                                                              <  1:3091037/100‑1 (MQ=255)
        tGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCgg                                                                                             <  2:1298813/100‑1 (MQ=255)
        tGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCgg                                                                                             <  1:3636541/100‑1 (MQ=255)
        tGAGCTTCGGGATGAAAAGCGATCACATAAAGAGCCACCGGGAGAGCAGGACGCTTCCTCAATATCTAAAAAAAAAAAAAAAACCCGCTTCAAAAAGc                                                                                               <  1:2482525/98‑1 (MQ=255)
         gATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGGCGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCggg                                                                                            >  2:2175016/1‑100 (MQ=255)
            cTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCggg                                                                                            >  1:3368152/1‑97 (MQ=255)
             ttCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtttt                                                                                        >  2:4008260/1‑100 (MQ=255)
              tCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttttt                                                                                       <  1:3859923/100‑1 (MQ=255)
              tCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtttt                                                                                        >  2:2926895/1‑99 (MQ=255)
                aGGCTGAAAAGCTCTCTCATAAAGAGTGACGTTGCTACCAGGACTCTCCCTCAATCTCTTGTTCGCTACAAACAACCCGCTTCAAAAAGCGGGtttt                                                                                        <  1:1778526/97‑1 (MQ=255)
                aGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTa                                                                                     >  1:3416020/1‑100 (MQ=255)
                 ggatGAAAAGCTCTCTCATAAAGAGAGACGTGGGCAGAAGGACGCTCCCTCAATATCTTATTCGAAACAACAAACGCGCTTCAAAAAGCGGGATTGTTa                                                                                     >  1:4282632/4‑99 (MQ=255)
                 ggCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttttttg                                                                                     <  1:2432013/99‑2 (MQ=37)
                 ggCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTAt                                                                                    <  2:1766402/100‑1 (MQ=255)
                  gCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATc                                                                                   <  2:161418/100‑1 (MQ=255)
                   cTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttttttgtat                                                                                  <  2:1824848/100‑5 (MQ=37)
                    tGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcag                                                                                 >  2:2748973/1‑100 (MQ=255)
                     gAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGTTACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcag                                                                                 <  1:666242/99‑1 (MQ=255)
                           ctctctCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATg                                                                          <  2:3154368/100‑1 (MQ=255)
                           ctctctCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATg                                                                          >  2:170408/1‑100 (MQ=255)
                             ctctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgta                                                                        >  2:2349638/1‑100 (MQ=255)
                             ctctcATAAAGAGTGACGTGGCGAGCAGGACGCCCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgta                                                                        <  2:473992/100‑1 (MQ=255)
                             ctctcATAAAGAGTGACGTGGCGAGCAGAACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgta                                                                        <  2:877335/100‑1 (MQ=255)
                               ctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaag                                                                      <  2:1423324/100‑1 (MQ=255)
                               ctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaag                                                                      >  2:2614129/1‑100 (MQ=255)
                                tcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagta                                                                    <  2:4031409/100‑1 (MQ=255)
                                 cATAAAGAGTGACTTGTCAAGAAGTACGCACACTCAAGATCTTGTTCGATACAAAAAACCCGCGTCAAAAAGCGGGTTTTTTATCAGACAGATgtaag                                                                      <  1:374178/98‑1 (MQ=255)
                                  aTAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagtaa                                                                   >  2:550328/1‑100 (MQ=255)
                                  aTAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagtaa                                                                   >  1:2105961/1‑100 (MQ=255)
                                  aTAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagta                                                                    >  2:2939631/1‑99 (MQ=255)
                                      agagTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAattat                                                               >  2:3745713/1‑100 (MQ=255)
                                      agagTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAattat                                                               <  2:208850/100‑1 (MQ=255)
                                          tGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATAAAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACa                                                           <  2:37611/100‑1 (MQ=255)
                                           gACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAg                                                          <  2:2482617/100‑1 (MQ=255)
                                           gACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAg                                                          >  1:2545404/1‑100 (MQ=255)
                                            aCGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAgg                                                         <  2:381371/100‑1 (MQ=255)
                                             cGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                        >  2:123945/1‑100 (MQ=255)
                                                ggCGAGCAGGACGCTCGCTCAATATCTTGTTCGATACAAGAAACCCGCTTCAAAAAGGGGGTTTTTCCTCAGACAGATGTAAGTAATTATTACAGGATTa                                                     >  1:3318054/1‑100 (MQ=255)
                                                ggCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTa                                                     >  1:1897168/1‑100 (MQ=255)
                                                ggCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTa                                                     >  1:3220771/1‑100 (MQ=255)
                                                  cGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACt                                                   >  2:414022/1‑100 (MQ=255)
                                                  cGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACt                                                   >  2:4273227/1‑100 (MQ=255)
                                                  cGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACt                                                   >  2:1240167/1‑100 (MQ=255)
                                                   gAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACtt                                                  <  1:550326/100‑1 (MQ=255)
                                                   gAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATAATTACAgg                                                         >  1:4356334/1‑93 (MQ=255)
                                                    aGCAGGACGCTCCTTCAATATCTTGTTTGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTa                                                 <  1:159264/100‑1 (MQ=255)
                                                      cAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAAc                                               <  2:3416038/100‑1 (MQ=255)
                                                        ggACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACtt                                             <  2:1961942/100‑1 (MQ=255)
                                                             cTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                         >  1:4017529/1‑99 (MQ=255)
                                                              tCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATcc                                       <  2:524529/100‑1 (MQ=255)
                                                              tCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATc                                        <  2:1353696/99‑1 (MQ=255)
                                                                ccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCt                                     >  1:1687981/1‑100 (MQ=255)
                                                                ccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCt                                     >  2:1597462/1‑100 (MQ=255)
                                                                ccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGAAAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                         >  1:780172/1‑96 (MQ=255)
                                                                   cAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTc                                  >  2:1540366/1‑100 (MQ=255)
                                                                    aaTATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTATTTATTACAGGATTACTTAACTTCCATCCCTTTct                                 <  1:293746/100‑2 (MQ=255)
                                                                    aaTATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCg                                 <  1:1384996/100‑1 (MQ=255)
                                                                       atCTTCTTCGATACAAAAAAAAAAAAAAAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGcc                               <  1:1066847/99‑1 (MQ=255)
                                                                        tcgtGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTg                             <  1:4338535/97‑1 (MQ=255)
                                                                         cTTGTTCGTTACAAACAACCCTCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTg                             <  1:1656593/99‑1 (MQ=255)
                                                                          ttGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCa                           <  1:2505622/100‑1 (MQ=255)
                                                                          ttGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCa                           >  2:571473/1‑100 (MQ=255)
                                                                              tCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATc                       >  2:667794/1‑100 (MQ=255)
                                                                                gATACAAAAAACCCGCTTCCAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCgg                     <  1:891429/100‑1 (MQ=255)
                                                                                ctgataaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCgg                     >  1:1733585/6‑100 (MQ=25)
                                                                                 tgataaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGc                    >  2:2367538/5‑100 (MQ=37)
                                                                                 tgataaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGc                    >  1:3613187/5‑100 (MQ=37)
                                                                                  tACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCg                   >  2:2551076/1‑100 (MQ=255)
                                                                                  gataaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCg                   >  1:3739725/4‑100 (MQ=37)
                                                                                  gataaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCg                   >  2:2521908/4‑100 (MQ=37)
                                                                                  gataaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCg                   <  2:939433/97‑1 (MQ=255)
                                                                                    taaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTTCTTATTACAGGATTACTCAACTTACATCCATTTCGCTTGCAAATCGGCGTg                 >  1:597307/2‑100 (MQ=25)
                                                                                    taaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                 >  2:3896178/2‑100 (MQ=38)
                                                                                    taaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTAAACTTCAATCCCTTTAGCCTGCAAATCGGCGt                  >  1:4160333/2‑99 (MQ=25)
                                                                                    cAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                 <  2:765940/100‑1 (MQ=255)
                                                                                    cAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                 >  2:2255266/1‑100 (MQ=255)
                                                                                    cAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                 <  2:3488660/100‑1 (MQ=255)
                                                                                    cAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                 >  2:436371/1‑100 (MQ=255)
                                                                                    cAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                 <  2:3933460/100‑1 (MQ=255)
                                                                                     aaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTgg                >  2:3879321/1‑100 (MQ=38)
                                                                                     aaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTAGCCTGCAAATCGGCGt                  >  1:4271475/1‑98 (MQ=37)
                                                                                     aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTgg                <  2:3235691/100‑1 (MQ=255)
                                                                                      aaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGt               >  1:1814915/1‑100 (MQ=38)
                                                                                      aaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGt               >  2:3681631/1‑100 (MQ=255)
                                                                                       aaaaCCCGCTTTTTGAAGCGGGCTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGt               >  1:3322503/1‑99 (MQ=38)
                                                                                       aaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTAATACAGGATTACTTAACATCCATCCCTTGCGCCTACAAATCGGCGTGGTa              >  1:967889/1‑100 (MQ=255)
                                                                                       aaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAAATATTACAGGATTACTTAACTTCCATCACTTTCGCCTACAAATCGGCGTGGTa              >  1:2282789/1‑100 (MQ=255)
                                                                                           cccGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaagaa          >  2:353495/1‑100 (MQ=255)
                                                                                           ccagCTTCAAGAAGCGGGTAATGAATAAGACAGATGTAAGTAATAATTACAGGATTACTTAAATAAAATTCATATCACCTGCAAATGGGAGTGGTaagaa          >  1:4375952/4‑100 (MQ=255)
                                                                                             cGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAagaaga        >  2:1547083/1‑100 (MQ=255)
                                                                                             cGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAagaaga        <  2:3319120/100‑1 (MQ=255)
                                                                                             cGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAagaaga        >  2:106231/1‑100 (MQ=255)
                                                                                                 tCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTATTTAACTTCCATCCCTTTCTCCTGCAAATCGGCGTGGTAAGAACAGCg     >  1:1942215/1‑99 (MQ=255)
                                                                                                 tCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCgg    >  1:3782450/1‑100 (MQ=255)
                                                                                                  cAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGa   >  2:2513299/1‑100 (MQ=255)
                                                                                                  cAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGa   >  2:2455201/1‑100 (MQ=255)
                                                                                                   aaaaaGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGAc  <  2:3782474/100‑1 (MQ=255)
                                                                                                   aaaaaGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGAc  >  1:4240601/1‑100 (MQ=255)
                                                                                                   aaaaaGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGAc  >  1:1058643/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGAC  >  NC_000913/659106‑659304

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: