Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,182,934 C→A 41.9% P564T (CCT→ACT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,9340CA41.9% 25.6 / 57.2 64P564T (CCT→ACT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base A (16/11);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAA  >  NC_000913/4182789‑4183082
                                                                                                                                                 |                                                                                                                                                    
ggACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACt                                                                                                                                                    >  6:95739/1‑148 (MQ=255)
 gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa                                                                                                                                                  <  3:355965/149‑1 (MQ=255)
    cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAgg                                                                                                                                               >  5:252402/1‑149 (MQ=255)
            cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                                                                                                                                        >  1:225188/1‑148 (MQ=255)
                 tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCTGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGt                                                                                                                                  >  8:8743/1‑149 (MQ=255)
                   tCTGATATTACGCACAAACGTCGTATCTCAGCACGCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCAAATCGAAACCCCTGAAGGTCCGAACATCGGTc                                                                                                                                 <  5:228112/148‑1 (MQ=255)
                   tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGaa                                                                                                                                         <  3:198199/140‑1 (MQ=255)
                   tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGaa                                                                                                                                         >  4:198199/1‑140 (MQ=255)
                     tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGa                                                                                                                              >  7:319594/1‑149 (MQ=255)
                     tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGa                                                                                                                              >  2:45437/1‑149 (MQ=255)
                       agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAt                                                                                                                             <  3:160632/148‑1 (MQ=255)
                       agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGGCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGa                                                                                                                              <  4:74227/147‑1 (MQ=255)
                          ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAATTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGTTCTGATCAAc                                                                                                                         <  2:58254/149‑1 (MQ=255)
                          ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAAc                                                                                                                         <  7:341107/149‑1 (MQ=255)
                              gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTGCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctc                                                                                                                     >  6:155418/1‑149 (MQ=255)
                              gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAcccc                                                                                                                                                     <  8:238608/117‑1 (MQ=255)
                              gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAcccc                                                                                                                                                     >  7:238608/1‑117 (MQ=255)
                              gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctc                                                                                                                     >  6:82060/1‑149 (MQ=255)
                                acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct                                                                                                                    >  8:350035/1‑148 (MQ=255)
                                acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTg                                                                                                                   >  4:262927/1‑149 (MQ=255)
                                acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTg                                                                                                                   >  7:250354/1‑149 (MQ=255)
                                acaAACGTCGTATCTCCGCACTAGGCCCAGGCGGTCTGACCCGGGATCGTGCAGGCTTCGAAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTACAATAGAAAGAACTGAAAGTCCGAACATCGGTCGGATAACCTCTCTg                                                                                                                   >  6:258810/1‑149 (MQ=255)
                                     cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt                                                                                                              <  2:225188/149‑1 (MQ=255)
                                        cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTa                                                                                                           >  3:165503/1‑149 (MQ=255)
                                        cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTa                                                                                                           >  3:69606/1‑149 (MQ=255)
                                         gtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAc                                                                                                          >  2:197331/1‑149 (MQ=255)
                                                 gCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                                                                                                  <  5:82060/149‑1 (MQ=255)
                                                       ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                              <  5:357954/147‑1 (MQ=255)
                                                       ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                              <  6:252402/147‑1 (MQ=255)
                                                       ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                              <  8:250354/147‑1 (MQ=255)
                                                         cccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                              >  1:149429/1‑145 (MQ=255)
                                                         cccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                              <  2:149429/145‑1 (MQ=255)
                                                              gCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                                                                                     >  6:95730/1‑149 (MQ=255)
                                                              gCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACACCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                                                                                     <  7:350035/149‑1 (MQ=255)
                                                                 gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                                                                                  <  5:95739/149‑1 (MQ=255)
                                                                        ccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgaga                                                                           <  3:262927/149‑1 (MQ=255)
                                                                             aaCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                                                                      >  6:87525/1‑149 (MQ=255)
                                                                                      gCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                                             <  1:45437/149‑1 (MQ=255)
                                                                                           gAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAc                                                        >  1:256732/1‑149 (MQ=255)
                                                                                            aaGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc                                                       <  8:319594/149‑1 (MQ=255)
                                                                                             aGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGCCTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc                                                       >  6:239253/1‑148 (MQ=255)
                                                                                                  gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGgt                                                 >  2:42858/1‑149 (MQ=255)
                                                                                                   agaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGgtg                                                >  7:344441/1‑149 (MQ=255)
                                                                                                    gtcGTACACCCGACTCTCTACGGTCGCGTATGTCCAATCGAAACCCCTGTAGGTCCGAACTTCTGTCTTTTCAACTCTCTGTCCGTTTACGcaca                                                                                                     <  2:369516/93‑1 (MQ=255)
                                                                                                    gaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGgtgt                                               >  6:352032/1‑149 (MQ=255)
                                                                                                    gaCGTACACCCGACACACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGACTGATCAACTCTCTGTCCGTGTAAGcaca                                                                                                     >  1:369516/1‑95 (MQ=255)
                                                                                                          cacCCGACTCACTGCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAAc                                         <  5:155418/149‑1 (MQ=255)
                                                                                                             ccGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGCCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGa                                      >  7:270654/1‑149 (MQ=255)
                                                                                                             ccGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGa                                      <  7:233314/149‑1 (MQ=255)
                                                                                                               gACTCTCTACGGTCGCGTTTGTGCAATCGGACCCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGAc                                     >  4:97754/1‑148 (MQ=255)
                                                                                                                    aCTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAt                                >  3:98694/1‑148 (MQ=255)
                                                                                                                      tACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCa                             >  6:60067/1‑149 (MQ=255)
                                                                                                                      tACGGTCGCGTATGTCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCa                             >  1:26112/1‑149 (MQ=255)
                                                                                                                         gggcgcgTCTGTCCAATCGAAACCACTGAAGGTCCGAACATATGGCTGATCACCGCTCTGTCCGTGTACGCACAGACTACCGACTCCGTCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACGGACGAAATTCCCTa                          <  5:258809/146‑1 (MQ=255)
                                                                                                                         gggcgcgTATGTCCAATCGAAACCACTGAATCTCCGACCATATGTCTGATCAACTCTCTCTCCGTGTACGCACAGACTAACGAATCCGGCTTCCATGAGACTCCGTATAGTACAGTGACCGACGGTGTTGTAACTGACGAAATTCACTa                          <  5:258810/146‑1 (MQ=255)
                                                                                                                         ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTa                          <  7:8743/149‑1 (MQ=255)
                                                                                                                           tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAcc                        >  8:23139/1‑149 (MQ=255)
                                                                                                                             tcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACg                                                   >  6:131113/2‑120 (MQ=255)
                                                                                                                             gcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACg                                                   <  5:131113/120‑1 (MQ=255)
                                                                                                                               gTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTc                    <  4:69606/149‑1 (MQ=255)
                                                                                                                                    tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTa               >  3:351071/1‑149 (MQ=255)
                                                                                                                                    tCCAATCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTa               <  5:87525/149‑1 (MQ=255)
                                                                                                                                       aaTCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATc             >  4:28316/1‑148 (MQ=255)
                                                                                                                                       aaTCGAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCAGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATc             >  6:282422/1‑148 (MQ=255)
                                                                                                                                           gAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCgaaga        <  1:42858/149‑1 (MQ=255)
                                                                                                                                           gAAACCACTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCgaag         <  1:197331/148‑1 (MQ=255)
                                                                                                                                                 ccTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCaa  <  2:256732/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAA  >  NC_000913/4182789‑4183082

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