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.

TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  NC_000913/3334435‑3334969
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                          
ggacggggaccaaggCACCCCAGAACAGGCCCATCTGTTCCGTACCGCCATCGAACAGGGTAACTGTCGCCATCTTCCGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcg                                                                                                                                                                                                                                                                       <  4:342541/261‑1 (MQ=255)
     aataacgaacacgacCGGAAGAGACACAGAGGGTCCGTAGCGCAATCGGACAGGGTAACGGTCGCCATCGTCGGGANCCGGTGCGGGAAGCAATGAACGCTTGTGGTTCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGGG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGa                                                                                                                                                                                                                                                                  <  3:163307/261‑1 (MQ=255)
                    cAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACATGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGACTCGCGCACATCGATGTTCAACGCGATCGTTAATGCCCTCCCCTCACCCCTCGCCGC‑CCGCGGGATGGTTCCAGTAAGCTc                                                                                                                                                                                                                                                  >  4:384669/1‑275 (MQ=255)
                          gggCCCAGATGACACGTACGGCCATCGAACAGGGTAACGGCCGCCATCTTCTGGAACCGGTGCTGGAAGCAAGGAAGGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGATAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAGCCCGCGCAGC‑GCGCGGGATGGTTCCAGGAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:483120/274‑1 (MQ=255)
                            gCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGACGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                                                                                                                                                      >  1:163308/1‑247 (MQ=255)
                                                                             tGGAACCGGTTCTGCAAGCAATGAACGCTGGTGGATCTCTTGAATTGACGCGTCAGCGTGCCGAGGAAGAAGCATACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGAACCCATCGGCCTCGCGAACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtat                                                                                                                                                                                         <  4:167921/275‑1 (MQ=255)
                                                                             gtgAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTGGAATGGAAGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtat                                                                                                                                                                                         <  2:586664/273‑1 (MQ=255)
                                                                                  acgggtagggAAGAAATTGACGTGTGTGGGTGTCTTGAATGGACGCGTCTGCGTGGCGAGGAATAAGTAGGCAAAGGCATCGCAGCGTTAGCGGTGACTCCTGG‑CACGCCATGGCGGGGAGCAGTCATCGGACTCGCGCAAATCGCTGGTCAACGCGATAGTTAATCCCCTCCCCTAATCCCGCGGAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTatatat                                                                                                                                                                                     <  4:112053/266‑1 (MQ=255)
                                                                                    ggTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                                >  2:640080/1‑213 (MQ=255)
                                                                                        cTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                            <  4:640078/213‑1 (MQ=255)
                                                                                               cAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcaca                                                                                                                                                                       <  2:368721/275‑1 (MQ=255)
                                                                                                aaTGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                                                                                                                                            >  2:168010/1‑189 (MQ=255)
                                                                                                                   ctTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  2:647980/1‑224 (MQ=255)
                                                                                                                       aaTGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  4:647978/224‑1 (MQ=255)
                                                                                                                              gcgTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                                                                                                                                                      >  8:13452/1‑149 (MQ=255)
                                                                                                                               actCAGCTGGCCAAGGAAAGAGCGGCCAAAGGCCTAGCCGCGTTACAGGTG‑CGCCCGGAAACCCCTGGGCGAGACGCACTCATCGGCCTGGCGCACATGCCGGGTCAACGCGATCGTTAATCCACACACCTCATCACACGCTGC‑GCGCGGGATGGTTCCAGTAATTTCCATAAAAACAGTTATTCAGCTCTGTTAAACCTGAAAAATAGGGATGTGCATCTTTATAAATTCTGAATATTCACACTCGTTAAAGGAACTTTTTAGAGCAATAGGcc                                                                                                                                       <  4:538669/273‑1 (MQ=255)
                                                                                                                                 tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCggg                                                                                                                                                                                                                                                                   >  8:626353/1‑149 (MQ=255)
                                                                                                                                  cAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                              >  4:599981/1‑169 (MQ=255)
                                                                                                                                     cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGATga                                                                                                                                                                                                                                                               <  8:78118/149‑2 (MQ=255)
                                                                                                                                      gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaaaa                                                                                                                                                                                                                                          <  2:599983/169‑1 (MQ=255)
                                                                                                                                          cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTcc                                                                                                                                                                                                                                                           >  5:91249/1‑148 (MQ=255)
                                                                                                                                           gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCa                                                                                                                                                                                                                                                          <  6:383705/148‑1 (MQ=255)
                                                                                                                                           gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                                                                                                                                                         <  7:872842/149‑1 (MQ=255)
                                                                                                                                              gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      >  7:792259/1‑149 (MQ=255)
                                                                                                                                              gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTa                                                                                                                                                                                                                                                       <  5:840507/148‑1 (MQ=255)
                                                                                                                                               aagaagCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                      <  6:91249/148‑1 (MQ=255)
                                                                                                                                                          aaaGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                                                                                                                                           <  5:773637/148‑1 (MQ=255)
                                                                                                                                                                aTCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                                                                                                                                                    >  8:490233/1‑149 (MQ=255)
                                                                                                                                                                 tCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                   >  8:440617/1‑149 (MQ=255)
                                                                                                                                                                  cGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCta                                                                                                                                                                                                                          >  1:33377/1‑157 (MQ=255)
                                                                                                                                                                      gCGTTACAGGTG‑CTCCCGNACACCCNTTGNCGAGAAGNACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAac                                                                                                                                                                                                                      <  3:33376/157‑1 (MQ=255)
                                                                                                                                                                       cGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGc                                                                                                                                                                                                                              >  6:490932/1‑148 (MQ=255)
                                                                                                                                                                        gTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                                                                                                                                                      >  2:656931/1‑107 (MQ=255)
                                                                                                                                                                         ttACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                                                                                                                                            <  6:64788/148‑1 (MQ=255)
                                                                                                                                                                          tACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                                                                                                                                                  <  2:84978/141‑1 (MQ=255)
                                                                                                                                                                            cAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGa                                                                                                                                                                                                                                                                  <  4:656929/107‑1 (MQ=255)
                                                                                                                                                                              ggTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                                       <  6:34227/148‑1 (MQ=255)
                                                                                                                                                                              ggTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCA‑‑‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAacac                                                                                                                                                                                                                    <  8:184883/149‑1 (MQ=255)
                                                                                                                                                                                   cTCCCGGACACCCCTTGGCGAGAAGAACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGa                                                                                     <  2:492689/274‑1 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                                 >  8:42784/1‑149 (MQ=255)
                                                                                                                                                                                    tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGC‑GCTGCGCGGGa                                                                                                                                                                                                                                                                  <  7:481627/100‑1 (MQ=255)
                                                                                                                                                                                      ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATg                                                                                                                                                                                                                                                                <  6:676817/100‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                                                                                                                                            <  8:582848/98‑1 (MQ=255)
                                                                                                                                                                                            acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                                                                                                                                            >  7:582848/1‑98 (MQ=255)
                                                                                                                                                                                             ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGATg                                                                                                                                                                                                                                                                <  5:690784/93‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      <  6:557603/148‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                                                                                                                                     <  8:83828/149‑1 (MQ=255)
                                                                                                                                                                                                cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCt                                                                                                                                                                                             >  1:1501/1‑157 (MQ=255)
                                                                                                                                                                                                  tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgt                                                                                                                                                                                                   <  7:523569/149‑1 (MQ=255)
                                                                                                                                                                                                  tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                                                                                                                                                                    <  5:462126/148‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcg                                                                                                                                                                                                                                                                       <  2:118825/79‑1 (MQ=255)
                                                                                                                                                                                                    gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:667252/1‑144 (MQ=255)
                                                                                                                                                                                                      gagaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                                                                                                                                            <  2:154758/73‑1 (MQ=255)
                                                                                                                                                                                                        gaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:667254/144‑1 (MQ=255)
                                                                                                                                                                                                          aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                                                                                                                                                      >  2:226709/1‑74 (MQ=255)
                                                                                                                                                                                                          aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAACCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCtt                                                                                                                                                                                            <  6:651944/148‑1 (MQ=255)
                                                                                                                                                                                                             aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                                                                                                                                         <  7:416409/132‑1 (MQ=255)
                                                                                                                                                                                                             aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGAGATGTGCATCTTtata                                                                                                                                                                                        >  7:573657/1‑149 (MQ=255)
                                                                                                                                                                                                              cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGa                                                                                                                                                                                                                                                                  <  4:226708/74‑1 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  >  6:618439/1‑137 (MQ=255)
                                                                                                                                                                                                               tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  5:618439/137‑1 (MQ=255)
                                                                                                                                                                                                                 aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  8:531183/1‑149 (MQ=255)
                                                                                                                                                                                                                  tCGGCGTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  7:149089/149‑1 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                    >  6:168541/1‑148 (MQ=255)
                                                                                                                                                                                                                  tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                   <  8:798121/149‑1 (MQ=255)
                                                                                                                                                                                                                    ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                                                                                                                                                 <  8:87641/149‑1 (MQ=255)
                                                                                                                                                                                                                      ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGTCTGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTTCATCTTTTNATATTCTGAATATTCACAATCTTAACCGGAACTTTTTAGAGTAATAGGCCATCTGGAGAATAGTGATTCTCGAANGAAGCAGTGTCTTGAATGAAAGCAAATAAANNAAGAAGGAAGGAAAAAGGGtaaaacattt                                                 >  3:88229/1‑268 (MQ=255)
                                                                                                                                                                                                                          gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAatat                                                                                                                                                                            <  5:157714/148‑1 (MQ=255)
                                                                                                                                                                                                                           cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                                                                                                                                            >  2:29526/1‑52 (MQ=255)
                                                                                                                                                                                                                           cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                      >  4:743242/1‑121 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcg                                                                                                                                                                                                                                                                       <  4:29526/52‑1 (MQ=255)
                                                                                                                                                                                                                               caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                  <  2:743244/121‑1 (MQ=255)
                                                                                                                                                                                                                                aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGCTGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTCATCCTCGACAGAAGAAGTTTTCGGAATGAAAGTGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTATTTTCATTGAATgg                                       >  4:18998/1‑275 (MQ=255)
                                                                                                                                                                                                                                 tCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAACTTCCATAAAAACACTTATTGAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCGCtctc                                                                                                                                                                    <  8:573657/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGAGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCa                                     <  1:88230/275‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGct                                                                                                                                                                                                                             >  4:447259/1‑89 (MQ=255)
                                                                                                                                                                                                                                        tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAgg                                                                                                                                         >  4:800805/1‑169 (MQ=255)
                                                                                                                                                                                                                                          aaCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg                                                                                                                                                             <  7:88329/147‑1 (MQ=255)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCggg                                                                                                                                                                                                                                                                   <  6:284221/44‑1 (MQ=18)
                                                                                                                                                                                                                                           aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCggg                                                                                                                                                                                                                                                                   >  5:284221/1‑44 (MQ=18)
                                                                                                                                                                                                                                             gcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                                                                                                                                                >  5:739974/1‑141 (MQ=255)
                                                                                                                                                                                                                                                 tCGTTAATACCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                    >  7:635756/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:403580/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                   <  8:635756/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTaga                                                                                                                                                 <  7:820312/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                                                                                                                                                  >  5:34046/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     tAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTagag                                                                                                                                                <  7:312982/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        tcccctccccTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCaa                                                                                                                                             >  8:279176/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ccctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAAt                                                                                                                                            >  5:602340/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                           >  5:413044/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                           >  5:668523/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                                                                                                                                                        <  5:649137/118‑1 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                                                                                                                                                 >  7:243461/1‑125 (MQ=255)
                                                                                                                                                                                                                                                            ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                                                                                                                                          >  6:820800/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCa                                                                                                                                              >  2:385288/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                   cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                    >  5:544009/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                   cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa                                                                                                                                   >  6:737285/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     tCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCAATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAAATTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCCTCCCGCAGACATCATTGCGGGttttggacgtaagg  >  4:582359/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGAt                                                                                                                      >  2:787312/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                          cgcAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                   >  2:452560/1‑45 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                          
TGCGATGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGG  >  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: ‑