Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,182,929 A→T 100% E562V (GAA→GTA)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,9290AT95.1% 337.6 / 5.2 101E562V (GAA→GTA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base T (51/45);  total (53/48)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 103 total aligned reads displayed.

TTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGC  >  NC_000913/4182784‑4183067
                                                                                                                                                 |                                                                                                                                          
tttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa                                                                                                                                          >  3:16225/1‑148 (MQ=255)
 ttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGtaa                                                                                                                                          <  6:211378/147‑1 (MQ=255)
     ggACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                                                                                                                     >  1:152005/1‑148 (MQ=255)
      gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCt                                                                                                                                     >  3:268838/1‑147 (MQ=255)
      gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCt                                                                                                                                     >  4:263986/1‑147 (MQ=255)
          aGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAgg                                                                                                                                >  7:47121/1‑148 (MQ=255)
             acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGt                                                                                                                               >  8:189122/1‑146 (MQ=255)
                acCCGCTGTCTGATATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   <  5:94590/139‑1 (MQ=255)
                acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGaa                                                                                                                          >  3:385768/1‑148 (MQ=255)
                 cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGATACGTACACACGACTCAATACGGTCGCGTATTTCAAATCGTAACCCCTGAAgagatatg                                                                                                                          >  6:94590/1‑140 (MQ=255)
                 cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGaa                                                                                                                          >  1:323886/1‑147 (MQ=255)
                        tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGACCAGGCGGTCTGACCCGGGAACTTTCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCgg                                                                                                                    <  5:430556/146‑1 (MQ=255)
                         cTGAGATTACGCACAAACGTCGTATCTCCGCACACGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGCCGTACAACAGACTCACTACCGTCGCGTATGTCCAATCGTAACCCCTGAAGTTCCGAACATCGGTCt                                                                                                                 >  5:223790/1‑148 (MQ=255)
                          tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCTCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCt                                                                                                                 >  8:142406/1‑147 (MQ=255)
                            agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGACTGACCCGTGAAAGAGCATGCTTCGAAGTTCGACACATACACCCGACTCACTACGGTCGAGTATGTCCAATCGTAACCCATGAAGGTCCGAACATCGGtatta                                                                                                               >  2:66355/1‑143 (MQ=255)
                              aTTACGCACAACCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGATCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAGGTCCAATCGTAACCCCGGAAGGTCCGACCATCGGTCGGagc                                                                                                             >  8:161331/1‑145 (MQ=255)
                               ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGCACGTGCAGGCTTCGAAGTTCGAGCCGGACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCaa                                                                                                           <  2:131233/148‑1 (MQ=255)
                                tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCTCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   >  8:160475/1‑123 (MQ=255)
                                tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   <  7:160475/123‑1 (MQ=255)
                                tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAAc                                                                                                          >  1:59015/1‑148 (MQ=255)
                                    cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAATATCgg                                                                                                                    >  3:188084/1‑134 (MQ=255)
                                    cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAATATCgg                                                                                                                    <  4:188084/134‑1 (MQ=255)
                                      caAACGTCGTATCTCCGCACTCGGCACATGCGGTCTGACACGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCCTATGTCCAATCGAAACCCCATAAGCTCCGAACAACTGTCTGATCAACTCTCTg                                                                                                    >  5:325660/1‑148 (MQ=255)
                                        aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   >  7:234352/1‑115 (MQ=255)
                                        aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                   <  8:234352/115‑1 (MQ=255)
                                        aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                                                                  <  4:180691/148‑1 (MQ=255)
                                          cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCg                                                                                                <  8:175955/148‑1 (MQ=255)
                                               tATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAc                                                                                           <  8:47121/148‑1 (MQ=255)
                                                aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCAGTGTACg                                                                                          >  4:70305/1‑148 (MQ=255)
                                                  ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                                                                         <  2:323886/147‑1 (MQ=255)
                                                    ccGCACTCGGCCCAGGCGGTCTTACACGTGAACGTGCAGGCTTCGAAGTTCGAGACGAACACCCGACTCACTACGGTCGCGTATGTCAAATCGTAACCCCTGAAGGTCCGAACATCTGTCTGCACAACTCTCTGTCCGTGTACGcaca                                                                                      >  7:261680/1‑148 (MQ=255)
                                                    ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                                                                        >  3:275212/1‑146 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  1:243684/146‑1 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  2:59015/146‑1 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  5:336745/146‑1 (MQ=255)
                                                     cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACACGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                                                                         <  3:263986/144‑1 (MQ=255)
                                                       cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                                                                                   >  7:13418/1‑148 (MQ=255)
                                                           cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                               >  7:385781/1‑148 (MQ=255)
                                                            ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                                                                <  1:384706/146‑1 (MQ=255)
                                                            ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                                                                <  7:204339/146‑1 (MQ=255)
                                                            ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACg                                                                              <  5:176547/148‑1 (MQ=255)
                                                             gctcAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                                                             <  7:142406/145‑1 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                                                            >  8:272703/1‑147 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                                                           <  4:385768/148‑1 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                                                           <  4:16225/148‑1 (MQ=255)
                                                               ccAGGCGGTCTGACCCGTGAACGTGCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCTTATGTCCAATCGTAACCCCTTAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                                                           <  7:384037/148‑1 (MQ=255)
                                                                cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                                                          <  7:189122/148‑1 (MQ=255)
                                                                cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                                                          >  4:342652/1‑148 (MQ=255)
                                                                  ggCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                                                        >  8:34779/1‑148 (MQ=255)
                                                                   gCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTACCCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                                                       >  7:432086/1‑148 (MQ=255)
                                                                      gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                                                    >  6:54693/1‑148 (MQ=255)
                                                                       tCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            <  3:116766/107‑1 (MQ=255)
                                                                       tCTGACCCGAGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            >  4:116766/1‑107 (MQ=255)
                                                                       gcTGACCCGTGAACGTGCAGGCTTCGAAGTCCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                                                                   <  8:385781/147‑1 (MQ=255)
                                                                        cTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                                                                  >  1:419717/1‑148 (MQ=255)
                                                                          gACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  3:314385/1‑148 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                                                                                 <  3:165095/130‑1 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                                                                                 >  4:165095/1‑130 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                                                               <  2:17625/148‑1 (MQ=255)
                                                                           aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                                                               <  4:268838/148‑1 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            <  5:41380/102‑1 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCa                                                                                                            >  6:41380/1‑102 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  3:23911/1‑146 (MQ=255)
                                                                            cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  2:68618/1‑146 (MQ=255)
                                                                                    cGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                                                       >  5:283253/1‑147 (MQ=255)
                                                                                       gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                                                             <  3:142656/106‑1 (MQ=255)
                                                                                       gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                                                             >  4:142656/1‑106 (MQ=255)
                                                                                           gCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa                                               <  2:152005/148‑1 (MQ=255)
                                                                                           gCTTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGACCATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa                                               <  3:338975/148‑1 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                >  4:129635/1‑130 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                <  3:129635/130‑1 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              >  4:186218/1‑148 (MQ=255)
                                                                                            cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              <  4:275212/148‑1 (MQ=255)
                                                                                              tCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt                                            >  8:121316/1‑148 (MQ=255)
                                                                                              tCGAAGTTAGATACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCTTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCATGAGACTCCGTATCGTaaa                                              >  1:311373/1‑146 (MQ=255)
                                                                                                gAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa                                          <  4:314385/148‑1 (MQ=255)
                                                                                                    ttCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGa                                      >  2:121037/1‑148 (MQ=255)
                                                                                                       gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACgg                                   <  3:70305/148‑1 (MQ=255)
                                                                                                            gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGtgttg                              <  3:342652/148‑1 (MQ=255)
                                                                                                             tACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              >  2:329676/1‑131 (MQ=255)
                                                                                                             tACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              <  1:329676/131‑1 (MQ=255)
                                                                                                               cacCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTa                            <  6:283253/147‑1 (MQ=255)
                                                                                                                  ccGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGGACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATAGTAAAGTGACCGACCGTGTTGTAACt                         >  5:210887/1‑147 (MQ=255)
                                                                                                                   cGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              <  7:325637/125‑1 (MQ=255)
                                                                                                                   cGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                              >  8:325637/1‑125 (MQ=255)
                                                                                                                   cGACTCACAACGGTCGCGAATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCACCTCTCTCTCCGTGTACGCACACACGAACGAATACGGGTTCCTTGCGACTCCGTATCGTAAATTGCCCCACGGTGTTGTAACTGa                       >  6:181742/1‑148 (MQ=255)
                                                                                                                         aCTACGGTCGCGTATGTACAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAt                 <  7:392348/148‑1 (MQ=255)
                                                                                                                          cTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                >  8:106435/1‑148 (MQ=255)
                                                                                                                          cTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                <  3:186218/148‑1 (MQ=255)
                                                                                                                           tACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTaa                           <  7:6287/136‑1 (MQ=255)
                                                                                                                           tACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                >  4:299667/1‑147 (MQ=255)
                                                                                                                             cGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTc               <  5:54693/146‑1 (MQ=255)
                                                                                                                             cGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTc               <  4:178048/146‑1 (MQ=255)
                                                                                                                                 cgcgTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTACCGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAcc         >  5:112802/1‑148 (MQ=255)
                                                                                                                                      aTGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       >  6:422406/1‑65 (MQ=255)
                                                                                                                                      aTGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                       <  5:422406/65‑1 (MQ=255)
                                                                                                                                       tGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTg   <  1:68618/148‑1 (MQ=255)
                                                                                                                                         tCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGc  >  2:195561/1‑147 (MQ=255)
                                                                                                                                          ccAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAc          >  5:361221/1‑138 (MQ=255)
                                                                                                                                          ccAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAc          <  6:361221/138‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
TTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGC  >  NC_000913/4182784‑4183067

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: