Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,704 G→C 5.6% intergenic (+23/‑205) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7040GC 5.6% 667.1 / 17.8 193intergenic (+23/‑205)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (82/98);  new base C (5/6);  total (88/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 207 total aligned reads displayed.

TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  NC_000913/3334435‑3334969
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                          
ggacggggaccaaggCACCCCAGAACAGGCCCATCTGTTCCGTACCGCCATCGAACAGGGTAACTGTCGCCATCTTCCGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  4:342541/261‑1 (MQ=255)
     aataacgaacacgacCGGAAGAGACACAGAGGGTCCGTAGCGCAATCGGACAGGGTAACGGTCGCCATCGTCGGGANCCGGTGCGGGAAGCAATGAACGCTTGTGGTTCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGGGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                                                                                                                                                                  <  3:163307/261‑1 (MQ=255)
                    cAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACATGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGACTCGCGCACATCGATGTTCAACGCGATCGTTAATGCCCTCCCCTCACCCCTCGCCGCCCGCGGGATGGTTCCAGTAAGCTc                                                                                                                                                                                                                                                  >  4:384669/1‑275 (MQ=255)
                          gggCCCAGATGACACGTACGGCCATCGAACAGGGTAACGGCCGCCATCTTCTGGAACCGGTGCTGGAAGCAAGGAAGGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGATAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAGCCCGCGCAGCGCGCGGGATGGTTCCAGGAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:483120/274‑1 (MQ=255)
                            gCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGACGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  1:163308/1‑247 (MQ=255)
                                                                             tGGAACCGGTTCTGCAAGCAATGAACGCTGGTGGATCTCTTGAATTGACGCGTCAGCGTGCCGAGGAAGAAGCATACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGAACCCATCGGCCTCGCGAACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtat                                                                                                                                                                                         <  4:167921/275‑1 (MQ=255)
                                                                             gtgAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTGGAATGGAAGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtat                                                                                                                                                                                         <  2:586664/273‑1 (MQ=255)
                                                                                  acgggtagggAAGAAATTGACGTGTGTGGGTGTCTTGAATGGACGCGTCTGCGTGGCGAGGAATAAGTAGGCAAAGGCATCGCAGCGTTAGCGGTGACTCCTGGCACGCCATGGCGGGGAGCAGTCATCGGACTCGCGCAAATCGCTGGTCAACGCGATAGTTAATCCCCTCCCCTAATCCCGCGGAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTatatat                                                                                                                                                                                     <  4:112053/266‑1 (MQ=255)
                                                                                    ggTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                                >  2:640080/1‑213 (MQ=255)
                                                                                        cTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:640078/213‑1 (MQ=255)
                                                                                               cAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcaca                                                                                                                                                                       <  2:368721/275‑1 (MQ=255)
                                                                                                aaTGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                                                                                                                                                                            >  2:168010/1‑189 (MQ=255)
                                                                                                                   ctTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  2:647980/1‑224 (MQ=255)
                                                                                                                       aaTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  4:647978/224‑1 (MQ=255)
                                                                                                                              gcgTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  8:13452/1‑149 (MQ=255)
                                                                                                                               actCAGCTGGCCAAGGAAAGAGCGGCCAAAGGCCTAGCCGCGTTACAGGTGCGCCCGGAAACCCCTGGGCGAGACGCACTCATCGGCCTGGCGCACATGCCGGGTCAACGCGATCGTTAATCCACACACCTCATCACACGCTGCGCGCGGGATGGTTCCAGTAATTTCCATAAAAACAGTTATTCAGCTCTGTTAAACCTGAAAAATAGGGATGTGCATCTTTATAAATTCTGAATATTCACACTCGTTAAAGGAACTTTTTAGAGCAATAGGcc                                                                                                                                       <  4:538669/273‑1 (MQ=255)
                                                                                                                                 tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   >  8:626353/1‑149 (MQ=255)
                                                                                                                                  cAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                              >  4:599981/1‑169 (MQ=255)
                                                                                                                                     cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATga                                                                                                                                                                                                                                                               <  8:78118/149‑2 (MQ=255)
                                                                                                                                      gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaaa                                                                                                                                                                                                                                          <  2:599983/169‑1 (MQ=255)
                                                                                                                                          cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTcc                                                                                                                                                                                                                                                           >  5:91249/1‑148 (MQ=255)
                                                                                                                                           gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCa                                                                                                                                                                                                                                                          <  6:383705/148‑1 (MQ=255)
                                                                                                                                           gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAg                                                                                                                                                                                                                                                         <  7:872842/149‑1 (MQ=255)
                                                                                                                                              gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      >  7:792259/1‑149 (MQ=255)
                                                                                                                                              gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTa                                                                                                                                                                                                                                                       <  5:840507/148‑1 (MQ=255)
                                                                                                                                               aagaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      <  6:91249/148‑1 (MQ=255)
                                                                                                                                                          aaaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                                                                                                                                           <  5:773637/148‑1 (MQ=255)
                                                                                                                                                                aTCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                                                                                                                                                    >  8:490233/1‑149 (MQ=255)
                                                                                                                                                                 tCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                   >  8:440617/1‑149 (MQ=255)
                                                                                                                                                                  cGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCta                                                                                                                                                                                                                          >  1:33377/1‑157 (MQ=255)
                                                                                                                                                                      gCGTTACAGGTGCTCCCGNACACCCNTTGNCGAGAAGNACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAac                                                                                                                                                                                                                      <  3:33376/157‑1 (MQ=255)
                                                                                                                                                                       cGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGc                                                                                                                                                                                                                              >  6:490932/1‑148 (MQ=255)
                                                                                                                                                                        gTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  2:656931/1‑107 (MQ=255)
                                                                                                                                                                         ttACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                                                                                                                                            <  6:64788/148‑1 (MQ=255)
                                                                                                                                                                          tACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                                                                                                                                                  <  2:84978/141‑1 (MQ=255)
                                                                                                                                                                            cAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                                                                                                                                                                  <  4:656929/107‑1 (MQ=255)
                                                                                                                                                                              ggTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                                       <  6:34227/148‑1 (MQ=255)
                                                                                                                                                                              ggTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCA‑‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAacac                                                                                                                                                                                                                    <  8:184883/149‑1 (MQ=255)
                                                                                                                                                                                   cTCCCGGACACCCCTTGGCGAGAAGAACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGa                                                                                     <  2:492689/274‑1 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                                 >  8:42784/1‑149 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGa                                                                                                                                                                                                                                                                  <  7:481627/100‑1 (MQ=255)
                                                                                                                                                                                      ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATg                                                                                                                                                                                                                                                                <  6:676817/100‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                                                                                                                                                                            <  8:582848/98‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTc                                                                                                                                                                                                                                                            >  7:582848/1‑98 (MQ=255)
                                                                                                                                                                                             ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATg                                                                                                                                                                                                                                                                <  5:690784/93‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      <  6:557603/148‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                                                                                                                                     <  8:83828/149‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCt                                                                                                                                                                                             >  1:1501/1‑157 (MQ=255)
                                                                                                                                                                                                  tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgt                                                                                                                                                                                                   <  7:523569/149‑1 (MQ=255)
                                                                                                                                                                                                  tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                                                                                                                                                                    <  5:462126/148‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  2:118825/79‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:667252/1‑144 (MQ=255)
                                                                                                                                                                                                      gagaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                                                                                                                                            <  2:154758/73‑1 (MQ=255)
                                                                                                                                                                                                        gaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:667254/144‑1 (MQ=255)
                                                                                                                                                                                                          aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcgc                                                                                                                                                                                                                                                                      >  2:226709/1‑74 (MQ=255)
                                                                                                                                                                                                          aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAACCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCtt                                                                                                                                                                                            <  6:651944/148‑1 (MQ=255)
                                                                                                                                                                                                             aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                                                                                                                                         <  7:416409/132‑1 (MQ=255)
                                                                                                                                                                                                             aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGAGATGTGCATCTTtata                                                                                                                                                                                        >  7:573657/1‑149 (MQ=255)
                                                                                                                                                                                                              cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGa                                                                                                                                                                                                                                                                  <  4:226708/74‑1 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  >  6:618439/1‑137 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  5:618439/137‑1 (MQ=255)
                                                                                                                                                                                                                 aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  8:531183/1‑149 (MQ=255)
                                                                                                                                                                                                                  tCGGCGTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  7:149089/149‑1 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  6:168541/1‑148 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  8:798121/149‑1 (MQ=255)
                                                                                                                                                                                                                    ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                                                                                                                                                 <  8:87641/149‑1 (MQ=255)
                                                                                                                                                                                                                      ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGTCTGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTTCATCTTTTNATATTCTGAATATTCACAATCTTAACCGGAACTTTTTAGAGTAATAGGCCATCTGGAGAATAGTGATTCTCGAANGAAGCAGTGTCTTGAATGAAAGCAAATAAANNAAGAAGGAAGGAAAAAGGGtaaaacattt                                                 >  3:88229/1‑268 (MQ=255)
                                                                                                                                                                                                                          gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAatat                                                                                                                                                                            <  5:157714/148‑1 (MQ=255)
                                                                                                                                                                                                                           cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                                                                                                                                            >  2:29526/1‑52 (MQ=255)
                                                                                                                                                                                                                           cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:743242/1‑121 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcgcg                                                                                                                                                                                                                                                                       <  4:29526/52‑1 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:743244/121‑1 (MQ=255)
                                                                                                                                                                                                                                aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGCTGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTCATCCTCGACAGAAGAAGTTTTCGGAATGAAAGTGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTATTTTCATTGAATgg                                       >  4:18998/1‑275 (MQ=255)
                                                                                                                                                                                                                                 tCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAACTTCCATAAAAACACTTATTGAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCGCtctc                                                                                                                                                                    <  8:573657/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGAGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCa                                     <  1:88230/275‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGct                                                                                                                                                                                                                             >  4:447259/1‑89 (MQ=255)
                                                                                                                                                                                                                                        tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg                                                                                                                                         >  4:800805/1‑169 (MQ=255)
                                                                                                                                                                                                                                          aaCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg                                                                                                                                                             <  7:88329/147‑1 (MQ=255)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   <  6:284221/44‑1 (MQ=18)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCggg                                                                                                                                                                                                                                                                   >  5:284221/1‑44 (MQ=18)
                                                                                                                                                                                                                                             gcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                                                                                                                                                >  5:739974/1‑141 (MQ=255)
                                                                                                                                                                                                                                                 tCGTTAATACCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                    >  7:635756/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:403580/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:635756/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTaga                                                                                                                                                 <  7:820312/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                                                                                                                                                  >  5:34046/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     tAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTagag                                                                                                                                                <  7:312982/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        tcccctccccTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCaa                                                                                                                                             >  8:279176/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ccctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAAt                                                                                                                                            >  5:602340/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           cctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                           >  5:413044/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           cctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                           >  5:668523/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                                                                                                                                                        <  5:649137/118‑1 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                                                                                                                                                 >  7:243461/1‑125 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                                                                                                                                          >  6:820800/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             tcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCa                                                                                                                                              >  2:385288/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                   cATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                    >  5:544009/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                   cATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa                                                                                                                                   >  6:737285/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     tCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCAATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAAATTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCCTCCCGCAGACATCATTGCGGGttttggacgtaagg  >  4:582359/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGAt                                                                                                                      >  2:787312/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                          cgcAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                   >  2:452560/1‑45 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                          
TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  NC_000913/3334435‑3334969

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