Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,188,953 G→A E1202K (GAA→AAA)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,188,9530GA98.4% 419.1 / ‑3.2 124E1202K (GAA→AAA) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (66/56);  total (67/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGA  >  NC_000913/4188808‑4189098
                                                                                                                                                 |                                                                                                                                                 
ggCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAAc                                                                                                                                                >  3:151286/1‑149 (MQ=255)
ggCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAAc                                                                                                                                                >  7:263729/1‑149 (MQ=255)
   aaTCGTGGCTGAAATGAGCGGTTTCGTTTCCTTCGGTAAAGAAACCTAAGGTATAGGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTAGGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtg                                                                                                                                             <  1:185181/149‑1 (MQ=255)
   aaTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGAGTCCGAAATGGCGGCCGCTCAACGGGGGCGAAGGTAACCgtg                                                                                                                                             >  6:254343/1‑149 (MQ=255)
    aTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtgt                                                                                                                                            <  5:218445/149‑1 (MQ=255)
    aTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtgt                                                                                                                                            <  5:358525/149‑1 (MQ=255)
    aTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtg                                                                                                                                             >  6:42895/1‑148 (MQ=255)
     tCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTa                                                                                                                                           <  6:97092/149‑1 (MQ=255)
       cTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAg                                                                                                                                          >  3:158202/1‑148 (MQ=255)
        tGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGaa                                                                                                                                        >  5:197098/1‑149 (MQ=255)
        tGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGaa                                                                                                                                        >  4:129360/1‑149 (MQ=255)
            tGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCCCCCCGGTAGACAGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtg                                                                                                                                    >  1:363696/1‑149 (MQ=255)
            tGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtg                                                                                                                                    >  2:324270/1‑149 (MQ=255)
             gAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGGCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTGCGAAGGTAAACGTGTAGAACgtgg                                                                                                                                   >  6:279204/1‑149 (MQ=255)
              aaaTCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtggt                                                                                                                                  >  8:241257/1‑149 (MQ=255)
              aaaTCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtg                                                                                                                                    <  4:151286/147‑1 (MQ=255)
              aaaTCAGCGGTAATGTTTCCTTCGGTAAATAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATACGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtgg                                                                                                                                   <  2:78366/148‑1 (MQ=255)
                 tCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGATATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGAc                                                                                                                               >  1:359738/1‑149 (MQ=255)
                      ggaaTCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCACGGTAGACCGTAGCGATCCGTGCGAAGAGATGATTCCGAAATGGCTTCACATCAACGTGTTCTAAGGTAAAAGTCTAGAACGTGCTGAAGTaaa                                                                                                                          >  8:207196/4‑148 (MQ=255)
                      ggTATCGTTTCCGTCGGTAAAGAAACCAAAGGGTCACGTCGGCGGGTTATGTCCCCGGGAGCCGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTCTCGTAAt                                                                                                                          >  3:382285/1‑149 (MQ=255)
                          tCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTcc                                                                                                                      <  1:324270/149‑1 (MQ=255)
                          tCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTcc                                                                                                                      >  7:297858/1‑149 (MQ=255)
                                  ttCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCGCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTcc                                                                                                              >  6:315364/1‑149 (MQ=255)
                                  ttCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTcc                                                                                                              >  7:338647/1‑149 (MQ=255)
                                    cGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTcc                                                                                                              <  4:368982/147‑1 (MQ=255)
                                       tAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAg                                                                                                         >  4:284875/1‑149 (MQ=255)
                                       tAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAg                                                                                                         >  8:190240/1‑149 (MQ=255)
                                              aCCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTTGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCa                                                                                                  <  6:92688/149‑1 (MQ=255)
                                               ccAAAGTTAAACGTCGTCTGGTTATCAACCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGCTGACGTAATTTCCGACGGTCAGGAAGCTCAGCAc                                                                                                 >  7:456435/1‑149 (MQ=255)
                                                cAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCAc                                                                                                 >  7:344933/1‑148 (MQ=255)
                                                cAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCAc                                                                                                 >  5:279514/1‑148 (MQ=255)
                                                  aaGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACATGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCg                                                                                                             >  7:325694/1‑134 (MQ=255)
                                                   aGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACa                                                                                             <  3:129360/149‑1 (MQ=255)
                                                     gTAAACGTCGTCTGGTTATCACCCGGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                           <  5:267620/149‑1 (MQ=255)
                                                     gTAAACGTCGTCTGGTTATCACCCCGGTCGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                           <  3:8955/149‑1 (MQ=255)
                                                     gTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                           <  3:2727/149‑1 (MQ=255)
                                                       aaaCGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACTAGGTGAACTAATTTACAAATATACGAAAGCGCCGCACGACATtct                                                                                         >  6:166976/1‑149 (MQ=255)
                                                       aaaCGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtct                                                                                         >  4:97686/1‑149 (MQ=255)
                                                         aCGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgc                                                                                       <  4:103019/149‑1 (MQ=255)
                                                           gtcgtcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgt                                                                                     >  8:244037/1‑149 (MQ=255)
                                                              gtcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgtct                                                                                   <  4:138884/148‑1 (MQ=255)
                                                               tcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgtctgc                                                                                 >  7:34903/1‑149 (MQ=255)
                                                               tcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgtctgc                                                                                 >  2:425768/1‑149 (MQ=255)
                                                                 tGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                               >  7:154419/1‑149 (MQ=255)
                                                                 tGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                               >  1:40224/1‑149 (MQ=255)
                                                                   gTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCgtgg                                                                             >  5:296672/1‑149 (MQ=255)
                                                                   gTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCgtg                                                                              <  1:28266/148‑1 (MQ=255)
                                                                    ttATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCgtggt                                                                            >  6:62299/1‑149 (MQ=255)
                                                                    ttATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTACGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCAGGAAGCGCCGAACGACAGTCTGCGTCTGCgtggt                                                                            >  6:419615/1‑149 (MQ=255)
                                                                       tCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGgtgt                                                                          <  3:97686/148‑1 (MQ=255)
                                                                       tCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGgtgt                                                                          <  6:125973/148‑1 (MQ=255)
                                                                         aCCCGGGTAGACGGTAGCGATCGTTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCa                                                                       <  4:382285/149‑1 (MQ=255)
                                                                         aCCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCa                                                                       <  5:42895/149‑1 (MQ=255)
                                                                           cccGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATg                                                                     >  8:396624/1‑149 (MQ=255)
                                                                           cccGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATg                                                                     >  2:398198/1‑149 (MQ=255)
                                                                               gTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGt                                                                 <  4:158202/149‑1 (MQ=255)
                                                                                    cGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTc                                                            <  7:241257/149‑1 (MQ=255)
                                                                                      gTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGt                                                          <  5:257662/149‑1 (MQ=255)
                                                                                       tAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGtt                                                         <  6:197098/149‑1 (MQ=255)
                                                                                         gCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTa                                                        <  8:297858/148‑1 (MQ=255)
                                                                                         gCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTa                                                        >  2:115871/1‑148 (MQ=255)
                                                                                           gATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGATCTCGGTGATCTCCTTAACATAAAAAAAAAAAAGCCGCACGACATTATTCGTCTGCGTGGGGTTCATGCTGTTACTCGTTACAt                                                     >  4:435200/1‑149 (MQ=255)
                                                                                             tCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCg                                                   >  6:282133/1‑149 (MQ=255)
                                                                                               cGTACGAAGAGATGATTCCGAAATGGCGTCCGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGtt                                                 <  7:92005/149‑1 (MQ=255)
                                                                                               cGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGTGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGt                                                  <  4:4593/148‑1 (MQ=255)
                                                                                                 tACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTaa                                               <  7:190240/149‑1 (MQ=255)
                                                                                                  aCGAAGAGATGATTCCGAAATGGCGTCAGCTCACCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGAAGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGGTAAc                                              >  3:222005/1‑149 (MQ=255)
                                                                                                   cgacgaGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAAc                                              <  1:100319/148‑1 (MQ=255)
                                                                                                   agAAGAGATGATTCCGAAATGGCGGCAGCTCACCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAAc                                              <  5:279204/147‑1 (MQ=255)
                                                                                                    gAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACTACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACg                                             >  4:274574/1‑148 (MQ=255)
                                                                                                    gAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGt                                                                 <  7:437398/128‑1 (MQ=255)
                                                                                                     aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACg                                             >  5:129388/1‑147 (MQ=255)
                                                                                                          atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACa                                      <  5:315364/149‑1 (MQ=255)
                                                                                                             atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGa                                   <  2:23285/149‑1 (MQ=255)
                                                                                                              ttCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGAc                                  <  2:363696/149‑1 (MQ=255)
                                                                                                                ccGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACg                                 >  3:148418/1‑148 (MQ=255)
                                                                                                                 cGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                               <  1:184508/101‑1 (MQ=255)
                                                                                                                 cGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGAAGGTCCGGAAGCGCCGCACGACATtctgagtctgctt                                                                               >  2:184508/1‑99 (MQ=255)
                                                                                                                    aaTGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATAc                            <  2:401336/149‑1 (MQ=255)
                                                                                                                     aTGGCGTCAGCTCAACGTGGTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATAcc                           <  7:207196/149‑1 (MQ=255)
                                                                                                                        gCGTCAGCTCAAGGTGTTCGAATGTAAACGTGTAGAACGTTGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTTGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGt                         <  8:344933/148‑1 (MQ=255)
                                                                                                                        gCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGt                         <  5:245629/148‑1 (MQ=255)
                                                                                                                         cGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCt                       <  1:445810/149‑1 (MQ=255)
                                                                                                                           tCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGCCGGTCCGGAAGCGCCGCACGCCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                      >  5:364108/1‑148 (MQ=255)
                                                                                                                           tCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                      >  1:404764/1‑148 (MQ=255)
                                                                                                                              gCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAgg                  <  5:62299/149‑1 (MQ=255)
                                                                                                                               cTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                 >  1:106871/1‑149 (MQ=255)
                                                                                                                               cTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                 <  2:29169/149‑1 (MQ=255)
                                                                                                                                  aaCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt              >  8:389538/1‑149 (MQ=255)
                                                                                                                                   aCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt             <  8:34903/149‑1 (MQ=255)
                                                                                                                                   aCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt              >  6:247137/1‑148 (MQ=255)
                                                                                                                                      tgtTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCTCAAGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGAAGGGCGtt             >  3:57895/1‑146 (MQ=255)
                                                                                                                                       gtTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAg                                     <  3:116240/121‑1 (MQ=255)
                                                                                                                                        ttCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt        >  1:299038/1‑149 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTaa                                       <  2:52706/115‑2 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTaa                                       >  1:52706/1‑114 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGAGTCTTCGGGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa     <  3:274574/149‑1 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGAAGTAATTTCCGAAAGAAAGTAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa     >  4:97648/1‑149 (MQ=255)
                                                                                                                                             aGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTc                                                            <  6:263151/92‑1 (MQ=255)
                                                                                                                                              ggTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTCCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGa  <  2:450431/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGA  >  NC_000913/4188808‑4189098

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