Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,722,608 T→C 5.1% S886P (TCA→CCA)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,722,6080TC 5.1% 376.3 / 7.8 118S886P (TCA→CCA) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (50/62);  new base C (3/3);  total (53/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 118 total aligned reads displayed.

TCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACA  >  NC_000913/2722463‑2722755
                                                                                                                                                 |                                                                                                                                                   
tCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATTGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGTATCGTTGGGCTTACTCTAAATCTTGCCCAGCGCTAGGAATCAt                                                                                                                                                  <  7:8422/149‑1 (MQ=255)
tCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                                                                                                                                                  <  7:90921/149‑1 (MQ=255)
tCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                                                                                                                                                  <  6:41251/149‑1 (MQ=255)
  aCGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                                                                                                                                                  >  7:243339/1‑147 (MQ=255)
  aCGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                                                                                                                                                  >  3:99250/1‑147 (MQ=255)
  aCGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGa                                                                                                                                                <  3:20631/149‑1 (MQ=255)
  aCGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGa                                                                                                                                                <  8:435527/149‑1 (MQ=255)
       cTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAg                                                                                                                                           <  8:107028/149‑1 (MQ=255)
        tACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAgg                                                                                                                                          >  7:159254/1‑149 (MQ=255)
        tACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAgg                                                                                                                                          >  5:398742/1‑149 (MQ=255)
          cAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTAGAAGAGGTGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTa                                                                                                                                        >  6:531689/1‑149 (MQ=255)
           aaCGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCa                                                                                                                                                   >  6:109099/1‑137 (MQ=255)
           aaCGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCa                                                                                                                                                   <  5:109099/137‑1 (MQ=255)
           aaCGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTAc                                                                                                                                       >  6:469629/1‑149 (MQ=255)
           aaCGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTAc                                                                                                                                       <  4:99250/149‑1 (MQ=255)
            aCGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACt                                                                                                                                      >  7:329829/1‑149 (MQ=255)
             cGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTg                                                                                                                                     <  6:4722/149‑1 (MQ=255)
                 tGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGaaa                                                                                                                                 >  5:225340/1‑149 (MQ=255)
                   aaTGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGaaa                                                                                                                                 <  8:443509/147‑1 (MQ=255)
                    aTGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGt                                                                                                                              <  1:205952/149‑1 (MQ=255)
                     tGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt                                                                                                                             >  1:79470/1‑149 (MQ=255)
                        tATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGa                                                                                                                           <  6:398742/148‑1 (MQ=255)
                        tATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAc                                                                                                                          <  5:81181/149‑1 (MQ=255)
                        tATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAc                                                                                                                          <  6:225340/149‑1 (MQ=255)
                         aTTACGATGCAAAACTATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATAGTTGGGCTTACGCTAAATCTTGCCCAGAGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAc                                                                                                                          >  4:430596/1‑148 (MQ=255)
                          ttACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCa                                                                                                                        <  2:228357/149‑1 (MQ=255)
                             cGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACt                                                                                                                      <  8:159254/148‑1 (MQ=255)
                                 gCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAc                                                                                                                          >  1:332683/1‑140 (MQ=255)
                                 gCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAc                                                                                                                          <  2:332683/140‑1 (MQ=255)
                                  ccAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATCTCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGAGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATc                                                                                                                 <  6:144006/148‑1 (MQ=255)
                                  ccAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATc                                                                                                                 <  7:502442/148‑1 (MQ=255)
                                  ccAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCg                                                                                                                >  1:72933/1‑149 (MQ=255)
                                   cAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACt                                                                                                                      <  4:440700/142‑1 (MQ=255)
                                    aaaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCa                                                                                                                        <  3:148108/139‑1 (MQ=255)
                                    aaaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCgg                                                                                                               >  7:99529/1‑148 (MQ=255)
                                    aaaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCgg                                                                                                               >  6:229198/1‑148 (MQ=255)
                                     aaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAg                                                                                                                                           <  7:363559/119‑1 (MQ=255)
                                     aaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAg                                                                                                                                           >  8:363559/1‑119 (MQ=255)
                                     aaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTgg                                                                                                                                    <  4:165926/126‑1 (MQ=255)
                                     aaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTgg                                                                                                                                    >  3:165926/1‑126 (MQ=255)
                                     aaCAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGa                                                                                                             <  2:79470/149‑1 (MQ=255)
                                         aTCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTa                                                                                                                                        >  1:350628/1‑118 (MQ=255)
                                         aTCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTa                                                                                                                                        <  2:350628/118‑1 (MQ=255)
                                          tCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGgt                                                                                                        >  4:425867/1‑149 (MQ=255)
                                             tGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGgtgt                                                                                                      >  4:348435/1‑148 (MQ=255)
                                             tGGCATGGTGGCGCTAGCCAGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGtt                                                                                                     <  8:243339/149‑1 (MQ=255)
                                                cATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGttatt                                                                                                  <  8:201913/149‑1 (MQ=255)
                                                cATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGttat                                                                                                   >  8:429576/1‑148 (MQ=255)
                                                 aTGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGttatta                                                                                                 <  8:474738/149‑1 (MQ=255)
                                                    gtggCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTattatt                                                                                               <  2:72933/148‑1 (MQ=255)
                                                    gtggCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTattatt                                                                                               <  3:240995/148‑1 (MQ=255)
                                                     tggCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGc                                                                                             <  8:329829/149‑1 (MQ=255)
                                                           tAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAAt                                                                                                                  >  5:16890/1‑122 (MQ=255)
                                                           tAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAAt                                                                                                                  <  6:16890/122‑1 (MQ=255)
                                                           tAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTattatt                                                                                               >  2:151245/1‑141 (MQ=255)
                                                           tAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTattatt                                                                                               <  1:151245/141‑1 (MQ=255)
                                                             gCCCGCAAGCTCGGGTTTAACGTTGATATGCAGCTCGAAGAGGGGATCGTTGGTCTTACTATAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTTGATATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCtt                                                                                       <  3:211810/147‑1 (MQ=255)
                                                              cccGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATg                                                                                    <  3:425867/149‑1 (MQ=255)
                                                              cccGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATg                                                                                    <  7:429576/149‑1 (MQ=255)
                                                               ccGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgt                                                                                   >  5:46580/1‑149 (MQ=255)
                                                               ccGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgt                                                                                   >  8:252822/1‑149 (MQ=255)
                                                               ccGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgt                                                                                   >  8:238914/1‑149 (MQ=255)
                                                                  cAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATgtcgt                                                                                >  4:474809/1‑149 (MQ=255)
                                                                    aGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCt                                                                              >  4:107237/1‑149 (MQ=255)
                                                                       tCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAACCATGAGTAAGGTACTGGAAAt                                                                                                                                <  4:181088/96‑1 (MQ=255)
                                                                       tCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAACCATGAGTAAGGTACTGGAAAt                                                                                                                                >  3:188390/1‑96 (MQ=255)
                                                                       tCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAACCATGAGTAAGGTACTGGAAAt                                                                                                                                >  3:181088/1‑96 (MQ=255)
                                                                       tCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAACCATGAGTAAGGTACTGGAAAt                                                                                                                                <  4:188390/96‑1 (MQ=255)
                                                                       tCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAACCATGAGTAAGGTACTGGAAAt                                                                                                                                >  5:349994/1‑96 (MQ=255)
                                                                         gggTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCAtt                                                                         <  3:348435/149‑1 (MQ=255)
                                                                            tttAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGa                                                                                                             >  7:43485/1‑110 (MQ=255)
                                                                            tttAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGa                                                                                                             <  8:43485/110‑1 (MQ=255)
                                                                            tttAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGa                                                                                                             >  5:75750/1‑110 (MQ=255)
                                                                                  gTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTg                                                                        <  3:61534/141‑1 (MQ=255)
                                                                                   ttGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAg                                                                                                                                           >  7:190355/1‑73 (MQ=255)
                                                                                   ttGATATCCAGCTCGAAGAGGGGATAGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGa                                                               >  5:156473/1‑149 (MQ=255)
                                                                                         tCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCtt                                                          >  4:158015/1‑148 (MQ=255)
                                                                                          acAGCTCGAAGAGGGGATCGATGGGCTTCCGCTAAATCTTGACCAGCGCGAGCAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCTGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCa                                                        <  6:334741/148‑1 (MQ=255)
                                                                                           cAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCaa                                                       <  3:107237/149‑1 (MQ=255)
                                                                                            aGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAAt                                                      <  7:128289/149‑1 (MQ=255)
                                                                                            aGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAAt                                                      <  3:158015/149‑1 (MQ=255)
                                                                                              cTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGa                                                    <  5:229198/149‑1 (MQ=255)
                                                                                                 gAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGaa                                                   >  3:442393/1‑147 (MQ=255)
                                                                                                 gAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAAc                                                  <  7:252822/148‑1 (MQ=255)
                                                                                                       gggATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACagag                                              >  3:481527/1‑146 (MQ=255)
                                                                                                          aTCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATg                                                                                    <  6:291775/105‑1 (MQ=255)
                                                                                                          aTCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATg                                                                                    >  5:291775/1‑105 (MQ=255)
                                                                                                             gTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTTGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATg                                      <  5:531689/148‑1 (MQ=255)
                                                                                                                 ggCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACtg                                 <  3:289104/149‑1 (MQ=255)
                                                                                                                 ggCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACt                                  <  6:46580/148‑1 (MQ=255)
                                                                                                                       cGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGAtgt                           <  4:288173/149‑1 (MQ=255)
                                                                                                                        gCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTTAtgtt                          >  1:424316/1‑149 (MQ=255)
                                                                                                                          tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGgtg                                                                                                       <  6:120928/70‑1 (MQ=255)
                                                                                                                          tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGgtg                                                                                                       >  5:120928/1‑70 (MQ=255)
                                                                                                                          tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGAtgttgt                        >  4:344834/1‑149 (MQ=255)
                                                                                                                          tAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGAtgttg                         >  1:363911/1‑148 (MQ=255)
                                                                                                                                    ccAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAAttt                >  3:90807/1‑147 (MQ=255)
                                                                                                                                    ccAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTaa              <  3:344834/149‑1 (MQ=255)
                                                                                                                                            aGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGtaat       >  3:333785/1‑148 (MQ=255)
                                                                                                                                                aTCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACa  <  7:364648/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACCACTTAATCGGGACTGGTGTTATTATTGCCCGCTTATGTCGTCTGCATTGCACAGAGGACCCTTCAATGAACAGAGAAGAAATGCACTGTGATGTTGTCAAAATTTAAGCGTAATAAACA  >  NC_000913/2722463‑2722755

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