Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,852,303 C→A 5.1% R139R (CGC→CGA)  pncA → nicotinamidase/pyrazinamidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,852,3030CA 5.1% 300.6 / 2.2 88R139R (CGC→CGA) pncAnicotinamidase/pyrazinamidase
Reads supporting (aligned to +/- strand):  ref base C (40/43);  new base A (3/2);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCGTGAATATCCAGCCCCAGGACAG  >  NC_000913/1852157‑1852449
                                                                                                                                                  |                                                                                                                                                  
aGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGTCCGTCGGCACAAAACCTCTCTCGATGACTGGTTAcgc                                                                                                                                                    >  4:16908/1‑147 (MQ=255)
  tGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGAt                                                                                                                                                 >  2:241523/1‑148 (MQ=255)
  tGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGAt                                                                                                                                                 >  2:462673/1‑148 (MQ=255)
   gAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGAt                                                                                                                                                 >  2:124046/1‑147 (MQ=255)
     aGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCAt                                                                                                                                              <  8:328488/148‑1 (MQ=255)
     aGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCAt                                                                                                                                              <  7:73907/148‑1 (MQ=255)
      ggCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCAt                                                                                                                                              >  2:407277/1‑147 (MQ=255)
        cgcACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGa                                                                                                                                            >  6:382397/1‑147 (MQ=255)
         gcACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGaa                                                                                                                                           >  8:164414/1‑147 (MQ=255)
          cacaATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAAt                                                                                                                                         >  7:239943/1‑148 (MQ=255)
          cacaATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAAt                                                                                                                                         >  7:199315/1‑148 (MQ=255)
           acaATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATc                                                                                                                                        <  5:323372/148‑1 (MQ=255)
            caATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCg                                                                                                                                       <  7:279681/148‑1 (MQ=255)
             aaTTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGa                                                                                                                                      <  7:164414/148‑1 (MQ=255)
             aaTTACATCCGGTACTGCACCAAAACGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAACCCTCTCTCGATGACTGGTTACGCGATCATGAAATCGa                                                                                                                                      <  6:62859/148‑1 (MQ=255)
              aTTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGTGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGAt                                                                                                                                     >  1:3307/1‑148 (MQ=255)
              aTTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGAt                                                                                                                                     >  1:180764/1‑148 (MQ=255)
                    tCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTg                                                                                                                               <  8:324285/148‑1 (MQ=255)
                     ccGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGAGATCATGAAATCGATGAATTGa                                                                                                                              >  3:306186/1‑148 (MQ=255)
                      cGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGa                                                                                                                              <  7:86323/147‑1 (MQ=255)
                      cGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGAt                                                                                                                             >  5:32799/1‑148 (MQ=255)
                      cGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGAt                                                                                                                             <  2:227343/148‑1 (MQ=255)
                           cTGCACCAAAAAGCGCTCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTAAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGtt                                                                                                                         <  1:124046/147‑1 (MQ=255)
                           cTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGtt                                                                                                                         <  7:165328/147‑1 (MQ=255)
                           cTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGtt                                                                                                                         <  1:28181/147‑1 (MQ=255)
                           cTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTa                                                                                                                        >  5:157918/1‑148 (MQ=255)
                           cTGCACCAAAAAGCGATCGCAGCGGTGTCCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGtt                                                                                                                         <  2:46504/147‑1 (MQ=255)
                             ccACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATg                                                                                                                      >  2:222656/2‑148 (MQ=255)
                              ctccAAAAAGCGATCTCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTTCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTAt                                                                                                                       <  7:451136/144‑1 (MQ=255)
                              cACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTAt                                                                                                                       <  6:157918/146‑1 (MQ=255)
                              cACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTAt                                                                                                                       <  7:156413/146‑1 (MQ=255)
                                  aaaaaGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCt                                                                                                                 >  1:229886/1‑148 (MQ=255)
                                   aaaaGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTg                                                                                                                <  5:449313/148‑1 (MQ=255)
                                          aTCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTg                                                                                                         <  2:379933/148‑1 (MQ=255)
                                          aTCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTg                                                                                                         <  1:407277/148‑1 (MQ=255)
                                          aTCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTg                                                                                                         <  2:229886/148‑1 (MQ=255)
                                           tCGCAGCGGTGTTCCATAAAGGTGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGa                                                                                                        <  2:3307/148‑1 (MQ=255)
                                            cGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGAc                                                                                                       >  5:393634/1‑148 (MQ=255)
                                            cGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGAc                                                                                                       >  6:70242/1‑148 (MQ=255)
                                             ggaGCGGTGTTGCTTATAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGATAAGCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGGCTGGCTACTGACt                                                                                                      <  8:357711/146‑1 (MQ=255)
                                                        ccATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGaa                                                                                            <  1:241523/147‑1 (MQ=255)
                                                        ccATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGAGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAg                                                                                           >  7:61095/1‑148 (MQ=255)
                                                           tAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGttt                                                                                        <  6:32799/148‑1 (MQ=255)
                                                           tAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGttt                                                                                        <  3:317949/148‑1 (MQ=255)
                                                                   aaaaTCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCt                                                                                >  2:321620/1‑148 (MQ=255)
                                                                      aTCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAACACCTCTCTCGATGACTGGTTACGAGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGCCTATTGCGTGAAGTTTACCGTGCTggc                                                                             >  6:405727/1‑147 (MQ=255)
                                                                         cTTTATTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGTCTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGAcgc                                                                          <  3:16908/148‑1 (MQ=255)
                                                                           ttAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGAt                                                                                                                             >  5:110691/1‑95 (MQ=255)
                                                                           ttAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGAt                                                                                                                             <  6:110691/95‑1 (MQ=255)
                                                                           ttAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGt                                                                        >  3:207388/1‑148 (MQ=255)
                                                                             aGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGtt                                                                       <  8:239943/147‑1 (MQ=255)
                                                                             aGTTGACAGTTACAGTGCCTTTTGTGATAACGGCCGTCGGCAGAAAACCTCTCTCTATGACTGGTTACGAGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTTGACGCGtt                                                                       <  4:306186/147‑1 (MQ=255)
                                                                               ttGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTa                                                                      >  2:343881/1‑146 (MQ=255)
                                                                                tGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTg                                                                                  >  5:389363/1‑133 (MQ=255)
                                                                                tGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAg                                                                   >  8:25949/1‑148 (MQ=255)
                                                                                tGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAg                                                                   >  3:69282/1‑148 (MQ=255)
                                                                                tGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATAATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTg                                                                                  <  6:389363/133‑1 (MQ=255)
                                                                                  aCAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTAt                                                                                                                       <  2:470475/94‑1 (MQ=255)
                                                                                  aCAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTAt                                                                                                                       >  1:470475/1‑94 (MQ=255)
                                                                                   cAGTTTCAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCTATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTa                                                                <  5:382397/148‑1 (MQ=255)
                                                                                   cAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTa                                                                <  2:180764/148‑1 (MQ=255)
                                                                                     gTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAgg                                                              <  6:393634/148‑1 (MQ=255)
                                                                                       tACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAACCCTCTCGCGATGACTGGTTACGCCATAAGGAAATCGATGAATTTAAATTAATGGGCCTGGCTACAGACTATGGGGTTAATTTTACCGTGCTGTACGCTTTAAAGTAAGGt                                                             >  4:276653/1‑147 (MQ=255)
                                                                                            tGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGAGGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAg                                                       <  4:64068/148‑1 (MQ=255)
                                                                                             gCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAgg                                                      >  5:112190/1‑148 (MQ=255)
                                                                                                ttttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTa                                                    >  4:433259/1‑147 (MQ=255)
                                                                                                 tttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGcc                                                                                                                  <  3:396433/84‑1 (MQ=255)
                                                                                                 tttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGcc                                                                                                                  >  4:396433/1‑84 (MQ=255)
                                                                                                 tttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTAt                                                                                                    >  8:114988/1‑98 (MQ=255)
                                                                                                 tttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTAc                                                                     >  6:208503/1‑129 (MQ=255)
                                                                                                 tttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTAc                                                                     <  5:208503/129‑1 (MQ=255)
                                                                                                 tttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTaaa                                                  >  2:367404/1‑148 (MQ=255)
                                                                                                 tttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGAGGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTAt                                                                                                    <  7:114988/98‑1 (MQ=255)
                                                                                                  ttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTaa                                                   <  7:25949/146‑1 (MQ=255)
                                                                                                      gATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTg                                              <  1:321620/147‑1 (MQ=255)
                                                                                                       aTAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGAGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGa                                             <  8:61095/147‑1 (MQ=255)
                                                                                                            ggCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTTAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCg                                       <  5:405727/148‑1 (MQ=255)
                                                                                                                     cAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTc                              >  2:252307/1‑148 (MQ=255)
                                                                                                                     cAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTCTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTc                              <  5:70242/148‑1 (MQ=255)
                                                                                                                       gAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGt                            >  5:279414/1‑148 (MQ=255)
                                                                                                                        aaaaCCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTg                           >  7:154823/1‑148 (MQ=255)
                                                                                                                        aaaaCCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTg                           >  3:250911/1‑148 (MQ=255)
                                                                                                                        aaaaCCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGATGGAAGAGTTACAGTTATGTTATAAGGTAAACTTTATTACCTATGGTTGTCGTg                           >  5:185229/1‑148 (MQ=255)
                                                                                                                          aaCCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGc                         >  8:133450/1‑148 (MQ=255)
                                                                                                                                     aTGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCg                        >  4:52769/1‑138 (MQ=255)
                                                                                                                                     aTGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCg                        <  3:52769/138‑1 (MQ=255)
                                                                                                                                      tGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCGTGAATATCCAg             >  6:32464/1‑148 (MQ=255)
                                                                                                                                        aCTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCGTGAATATCCAGcc           >  5:91706/1‑148 (MQ=255)
                                                                                                                                               aCGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCGGGAATATCCAGCCCCAGGa     <  4:69282/147‑1 (MQ=255)
                                                                                                                                                cgcgATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCGTGAATATCCAGCCCCAGGAc    <  4:250911/147‑1 (MQ=255)
                                                                                                                                                 gcgATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCGTGAATATCCAGCCCCAGGACAg  <  1:367404/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGATTACCGATGGTTGTCGTGGCGTGAATATCCAGCCCCAGGACAG  >  NC_000913/1852157‑1852449

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