Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,524 A→T 6.3% K807M (AAG→ATG)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5240AT 6.3% 786.0 / 6.0 248K807M (AAG→ATG) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (99/124);  new base T (2/21);  total (101/147)
Fisher's exact test for biased strand distribution p-value = 6.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

TGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTC  >  NC_000913/164375‑164672
                                                                                                                                                     |                                                                                                                                                     
 gCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAg                                                                                                                                                      >  2:4757958/1‑150 (MQ=255)
  cccAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGt                                                                                                                                                     <  2:3367608/150‑1 (MQ=255)
  cccAAAGGCCATAACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCAGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGt                                                                                                                                                     >  1:5655994/1‑150 (MQ=255)
   ccAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTa                                                                                                                                                    <  1:1366866/150‑1 (MQ=255)
   ccAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTa                                                                                                                                                    >  2:795446/1‑150 (MQ=255)
    cAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAt                                                                                                                                                   <  1:927855/150‑1 (MQ=255)
    cAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAt                                                                                                                                                   >  2:3250576/1‑150 (MQ=255)
     aaaGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAtt                                                                                                                                                  <  2:1490827/150‑1 (MQ=255)
         gCCATTACAACTCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTa                                                                                                                                              <  1:1006197/150‑1 (MQ=255)
          ccATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTaa                                                                                                                                             <  1:112670/150‑1 (MQ=255)
          ccATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTaa                                                                                                                                             <  2:2480872/150‑1 (MQ=255)
             ttACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggt                                                                                                                                          >  1:3640935/1‑150 (MQ=255)
               aCAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggttg                                                                                                                                        <  2:3583182/150‑1 (MQ=255)
                 aaaTCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggttgg‑t                                                                                                                                     >  1:3054323/1‑150 (MQ=255)
                 aaaTCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggttgg‑t                                                                                                                                     <  1:986674/150‑1 (MQ=255)
                  aaTCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAggttgg‑tt                                                                                                                                    >  1:5029531/1‑150 (MQ=255)
                    tCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑tttt                                                                                                                                  <  1:2351454/150‑1 (MQ=255)
                    tCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑tttt                                                                                                                                  >  1:643176/1‑150 (MQ=255)
                     cacaCGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtc                                                                                                                                 <  1:4287242/150‑1 (MQ=255)
                     cacaCGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtc                                                                                                                                 <  1:4556358/150‑1 (MQ=255)
                      acacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtct                                                                                                                                <  1:2784599/150‑1 (MQ=255)
                      acacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtct                                                                                                                                <  2:1436202/150‑1 (MQ=255)
                      acacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtct                                                                                                                                <  1:4148478/150‑1 (MQ=255)
                      acacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtct                                                                                                                                <  1:1252350/150‑1 (MQ=255)
                      acacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtct                                                                                                                                <  1:5234333/150‑1 (MQ=255)
                      acacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtct                                                                                                                                <  1:5150790/150‑1 (MQ=255)
                       cacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtctc                                                                                                                               <  2:439620/150‑1 (MQ=255)
                       cacGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtctc                                                                                                                               <  1:5269925/150‑1 (MQ=255)
                        acGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtctcc                                                                                                                              >  1:335523/1‑150 (MQ=255)
                        acGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtctcc                                                                                                                              <  2:5581457/150‑1 (MQ=255)
                          gagaTTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtctccct                                                                                                                            <  1:4949811/150‑1 (MQ=255)
                          gagaTTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtctccct                                                                                                                            >  2:1801808/1‑150 (MQ=255)
                                  aGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCtgt                                                                                                                    <  2:945381/150‑1 (MQ=255)
                                   gCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTtct                                                                                                                                <  2:5021973/137‑1 (MQ=255)
                                   gCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCtgtg                                                                                                                   <  2:2302499/150‑1 (MQ=255)
                                    cGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTgg                                                                                                                  <  2:5160449/150‑1 (MQ=255)
                                    cGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTgg                                                                                                                  >  1:4379959/1‑150 (MQ=255)
                                    cGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAAACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTgg                                                                                                                  >  1:4857213/1‑150 (MQ=255)
                                     gACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTggg                                                                                                                 >  2:5327692/1‑150 (MQ=255)
                                     gACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTggg                                                                                                                 >  2:5327814/1‑150 (MQ=255)
                                      aCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCACACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGga                                                                                                                <  2:4052659/150‑1 (MQ=255)
                                      aCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGga                                                                                                                <  1:4757958/150‑1 (MQ=255)
                                       cTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGgag                                                                                                               <  1:1953666/150‑1 (MQ=255)
                                        ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGgaga                                                                                                              >  2:3216370/1‑150 (MQ=255)
                                        ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGgaga                                                                                                              <  1:4648254/150‑1 (MQ=255)
                                        tgcTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGgaga                                                                                                              <  2:5658460/148‑1 (MQ=255)
                                         tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGagag                                                                                                             <  2:2767325/150‑1 (MQ=255)
                                           tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                           >  2:5740297/1‑150 (MQ=255)
                                           tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                           >  2:2165673/1‑150 (MQ=255)
                                           tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCACAGGGAGAggg                                                                                                           <  1:4709432/150‑1 (MQ=255)
                                            ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                          >  2:1897792/1‑150 (MQ=255)
                                            ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCACAGGGAGAggga                                                                                                          <  2:844373/150‑2 (MQ=255)
                                             gAAAGGAGCGTAGCGTGAGGTGCAAAAAGATATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACTGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGTGAGAGGGcc                                                                                                         >  2:4096012/1‑150 (MQ=255)
                                              aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCGCCCTCTCCCACAGGGAGAgggaga                                                                                                        <  1:2428675/150‑4 (MQ=255)
                                               aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                       <  2:1373298/150‑1 (MQ=255)
                                                  gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                    >  1:3996347/1‑150 (MQ=255)
                                                    gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGATATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                  <  2:5802904/150‑1 (MQ=255)
                                                     cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                 <  2:1911111/150‑1 (MQ=255)
                                                     cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                 <  2:604008/150‑1 (MQ=255)
                                                     cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                 >  1:1594763/1‑150 (MQ=255)
                                                     cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                 <  2:1887413/150‑1 (MQ=255)
                                                     cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                 >  1:5805946/1‑150 (MQ=255)
                                                       tACCGTGAGGTGCAAAAAGTGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGG‑TTTTATCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                               <  2:293580/150‑1 (MQ=255)
                                                       tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCCACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAggg                                                                                               <  2:487794/150‑1 (MQ=255)
                                                          cGTGAGGTGCAAAAAGAGATGATAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                            <  2:2972408/150‑1 (MQ=255)
                                                          cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                            <  2:1860334/150‑1 (MQ=255)
                                                          cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                            >  2:3178212/1‑150 (MQ=255)
                                                           gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                           <  2:2742876/150‑1 (MQ=255)
                                                            tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                          >  1:2325166/1‑150 (MQ=255)
                                                                   cAAAATGAGATGAAAGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGAGGGACGAAAAAGTTTTCTTAAGGTTGT‑TTTTCTCCCTCTCCCTGTGGTTGAGGGCGGGGGTGAGGGCATCAGCGCGCa                                                                                   <  2:4297572/150‑1 (MQ=38)
                                                                      aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                >  1:5491807/1‑150 (MQ=255)
                                                                        agtgaTTAATGGGCGTTAGCCCAAACATGTCTGGCGGGACGACCGGGAAAATACTGCACCGACGCTACGGACGAAAATGTATTCGTATGGTTGG‑TTTTCTCCCTCTCCCTTTGGTAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTtc                                                                              <  2:3106498/147‑1 (MQ=255)
                                                                          agaTTAATGGGCGTTATCCCAAACATGTCTGGCGGGCCGCCCCGGAAAATACGGCACCGAGGCGACGGACGAAAATGTTTTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTTTGGTTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                            <  2:4777920/150‑1 (MQ=255)
                                                                           tttgattgggCTTTATCCCAAACTTGTGTGGCCGTACGACCCGGAAAATACTGCACGGACGCTAGGGACTAAAAATTTTTCTTATGTTTGT‑TTTTGACCCTCTCCCTGTGGTAGTGGGCCGGGGTTAGGGCATCAGCGCGCACGTtcacc                                                                           <  2:1224463/143‑1 (MQ=255)
                                                                              gAATGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGCCGGACGAAAAAGTATTCGTATGGTTGG‑TTTTCGCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccctc                                                                        <  2:921085/150‑1 (MQ=255)
                                                                              gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGtgacccac                                                                        >  2:1329106/1‑148 (MQ=255)
                                                                              gAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGCCAAATACGCCACCGACGCGACGCACCAAAAAGTATTCGTATGGTTGT‑TTTCCACCCACGCCGTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                        <  2:1611227/150‑1 (MQ=255)
                                                                                   ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTgcgccccccccca                                                                   >  1:3574960/1‑149 (MQ=38)
                                                                                      gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTTGGCCCGCCCCCGAAc                                                                >  1:4405240/1‑150 (MQ=38)
                                                                                          cccaaaaCATGTCGGGCCGGACGACCCGGCAAATCCGGCACCGCCGCGCCGGACCAAAAAGTATTCTTATGTTTGT‑TTTTCCCCCCCACCCTGTGGTTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                            <  2:5397475/146‑1 (MQ=255)
                                                                                            ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCactc                                                          >  1:1300100/1‑150 (MQ=255)
                                                                                            ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCCAAACCactc                                                          >  1:1860633/1‑150 (MQ=37)
                                                                                             cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGCCCCGCCCCCCAAcccacacc                                                         >  1:2553135/1‑145 (MQ=37)
                                                                                                  tggtggggccggccGCCCCGGAAAATACTGCCCCGCCGGGCCGTACTAAAATTTA‑TCGGTAGGTTGGTTTTTCTCCCCCCCCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                    <  1:5458848/144‑1 (MQ=255)
                                                                                                    gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCtgt                                                                                                                    >  1:4476037/1‑84 (MQ=255)
                                                                                                      cgggCCGGACGACCCGGCAAAGACGGCACCGACGCGACGGCCGAAAAAGTATTCGTACGGGTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                <  2:3840335/148‑1 (MQ=38)
                                                                                                           cGGACGACCCGGCAATTACGGCACCGACGCGACGGACGAAAATGTATTCGTATGGTTGG‑TTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                           <  2:2392270/150‑1 (MQ=39)
                                                                                                             gactccccGGAAAATCCGGCACCGACGCGCGGGACGAAAATTAATTGGTATGTTTGT‑TTTTCTCCCTCTCCCTGTGGGTTTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                         <  2:5515901/145‑1 (MQ=25)
                                                                                                                gccccGCCAAATCCTGCACCGACGCGCCGGACGAAAAAGTTTTCGTATGGGTGG‑TTGTCCCCCGCCCCCTGGGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                      <  2:4033996/148‑6 (MQ=255)
                                                                                                                  cccgcaaaaTACTCCCCCTACGCCACTGAAGAAAATTTTTTCTTACGTTTTT‑TTTTCCCCCTCCCCGTGTGGTTGTGGGCGGGTGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                    <  2:259604/144‑1 (MQ=255)
                                                                                                                   ccGGCAAAAACTGCACCGCCGCGACGCACTAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCCCCCTGTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                   <  2:4974867/150‑1 (MQ=255)
                                                                                                                      tcAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGT‑TTTTCTCCCACTCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                <  2:1981728/149‑1 (MQ=255)
                                                                                                                         aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCCACCAAACCCACAACCACCACGGGGAGGGGTCCGAACGAGccc                             >  1:1414439/1‑148 (MQ=21)
                                                                                                                                cACCGCGGCGACGGACTAAAACTTTTTGGTATGTTTGG‑TTTCCCCCCGCCCGCGGGGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                      <  1:2209324/150‑1 (MQ=255)
                                                                                                                                cACCGCCGCGCCGCACCAAAACTTATTCTAATGGTTGT‑TTTTCCCCCTCCCGCTGTGGTTGAGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                      <  1:3158459/150‑1 (MQ=255)
                                                                                                                                    gccgcgacgGACGAAAATGTATTGGTATGTTTGT‑TTTTCTCCCGCTCCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                  <  2:1690040/148‑1 (MQ=255)
                                                                                                                                            ggacgAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCACACCCAAAACCAAACCCACAAGGGAGAGGGGCCCGAAAGAGCCCAAAATTTTAGAAataaaat          >  2:5681089/1‑146 (MQ=21)
                                                                                                                                             gacgAAAATGTAGTCGTAAGGTTGG‑TTTTCGCCCTCTCCCTGTGGGTGAGGGCCGGGGGGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc         <  2:3617162/150‑1 (MQ=255)
                                                                                                                                               ctaaaaaGTATTCGTATGTTTGT‑TTTACTCCCCCTCCGTGTGGTTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt       <  1:1871481/148‑1 (MQ=255)
                                                                                                                                                    aaGTATTCGTTAGGTTGG‑TTTTCTCCCTCGCCCTGTGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAATTTTTGAGAGATATCTTCTTCTGTc  <  2:3910642/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTC  >  NC_000913/164375‑164672

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