Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,721,243 G→A 56.1% G431R (GGA→AGA)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,721,2430GA56.1% 45.0 / 155.4 123G431R (GGA→AGA) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base A (36/33);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GCAAGTATTAATTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTAC  >  NC_000913/2721101‑2721388
                                                                                                                                              |                                                                                                                                                 
gcAAGTATTAATTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCa                                                                                                                                             >  7:3371/1‑149 (MQ=255)
 cAAGTATTACTTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGGTTATTCAGCGAAGCCGGGCGGCCGACCTACCGTACCCCGGAAAGAACCAt                                                                                                                                            <  4:366047/149‑1 (MQ=255)
  aaGTATTAATTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATc                                                                                                                                           <  8:3371/149‑1 (MQ=255)
    gTATTAATTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCAc                                                                                                                                         >  1:140832/1‑149 (MQ=255)
     tATTAATTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACt                                                                                                                                        >  4:340341/1‑149 (MQ=255)
       ttAATTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGc                                                                                                                                      <  2:268937/149‑1 (MQ=255)
         aaTTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTTACGAACTGGTGCTGCTAGCACTCCTCGCAAGAGGCACTACGTTTATTCAGCGAAGCCGGGATGCCTAACAACCGTACCCCGGAAGGACCAATCACTGCtt                                                                                                                                    >  8:469171/1‑149 (MQ=255)
          aTTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCttt                                                                                                                                   >  2:72773/1‑149 (MQ=255)
          aTTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCtta                                                                                                                                   <  7:37148/149‑2 (MQ=255)
             gAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTa                                                                                                                                 <  7:469171/148‑1 (MQ=255)
              aaGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTAt                                                                                                                                <  3:279067/148‑1 (MQ=255)
               aGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATg                                                                                                                               >  1:324902/1‑148 (MQ=255)
                    tAAAGCATCATCCCCGCAGCAATTATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATAtg                                                                                                                         <  6:60856/149‑1 (MQ=255)
                    tAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATAtg                                                                                                                         >  3:273443/1‑149 (MQ=255)
                    tAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATAtg                                                                                                                         >  4:305243/1‑149 (MQ=255)
                    tAAAGAATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATAtg                                                                                                                         <  8:484894/149‑1 (MQ=255)
                       aGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGc                                                                                                                              <  3:384811/141‑1 (MQ=255)
                       aGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGc                                                                                                                              >  4:384811/1‑141 (MQ=255)
                       aGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGAAGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATAtggtg                                                                                                                      >  1:333621/1‑149 (MQ=255)
                        gCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATAtggtgg                                                                                                                     <  2:140832/149‑1 (MQ=255)
                         catcatCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGa                                                                                                                    <  8:142443/149‑1 (MQ=255)
                         catcatCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGa                                                                                                                    <  4:352105/149‑1 (MQ=255)
                         catcatCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCAt                                                                                                                                            >  3:433193/1‑125 (MQ=255)
                         catcatCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCAt                                                                                                                                            <  4:433193/125‑1 (MQ=255)
                                 ccGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATg                                                                                                                               >  2:372384/1‑130 (MQ=255)
                                 ccGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATg                                                                                                                               <  1:372384/130‑1 (MQ=255)
                                    cagcaAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTa                                                                                                         <  8:219473/149‑1 (MQ=255)
                                      gcaAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAAt                                                                                                       <  3:165937/149‑1 (MQ=255)
                                        aaaTATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCa                                                                                                     <  7:172486/149‑1 (MQ=255)
                                         aaTATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAAt                                                                                                       >  8:153022/1‑146 (MQ=255)
                                         aaTATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAAt                                                                                                       >  2:172252/1‑146 (MQ=255)
                                           tatGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGaa                                                                                                  <  1:72773/149‑1 (MQ=255)
                                           tatGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGaa                                                                                                  <  3:340341/149‑1 (MQ=255)
                                                ctctTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAAc                                                                                             >  7:294773/1‑149 (MQ=255)
                                                   tttGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTAt                                                                                                                                <  5:142172/111‑1 (MQ=255)
                                                   tttGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTAt                                                                                                                                >  6:142172/1‑111 (MQ=255)
                                                   tttGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTAc                                                                                          >  4:251976/1‑149 (MQ=255)
                                                        tGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGa                                                                                      >  8:108873/1‑148 (MQ=255)
                                                        tGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTACCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGa                                                                                      >  8:192172/1‑148 (MQ=255)
                                                        tGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATCCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGa                                                                                      >  5:324740/1‑148 (MQ=255)
                                                         gACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGaaa                                                                                    >  4:265092/1‑149 (MQ=255)
                                                         gACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGaaa                                                                                    >  8:105555/1‑149 (MQ=255)
                                                          aCGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCACGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGaa                                                                                     <  4:229102/147‑1 (MQ=255)
                                                            gAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGc                                                                                 >  4:76418/1‑149 (MQ=255)
                                                              aCTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTTCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCg                                                                               <  2:333621/149‑1 (MQ=255)
                                                              aCTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCg                                                                               <  3:251976/149‑1 (MQ=255)
                                                                  gTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGt                                                                           >  6:16034/1‑149 (MQ=255)
                                                                      ggcgAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGcc                                                                       <  1:166067/149‑1 (MQ=255)
                                                                        cgAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCa                                                                     <  3:265092/149‑1 (MQ=255)
                                                                               tcctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAAt                                                                >  1:98863/1‑147 (MQ=255)
                                                                               tcctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCt                                                              >  3:211168/1‑149 (MQ=255)
                                                                                 ctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCt                                                              >  5:224647/1‑147 (MQ=255)
                                                                                 ctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGa                                                            >  5:159814/1‑149 (MQ=255)
                                                                                 ctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCt                                                              <  7:105555/147‑1 (MQ=255)
                                                                                 ctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGa                                                            >  8:259212/1‑149 (MQ=255)
                                                                                 ctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGa                                                            >  2:135559/1‑149 (MQ=255)
                                                                                 ctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGa                                                            >  2:416717/1‑149 (MQ=255)
                                                                                 ctcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGa                                                            >  7:487430/1‑149 (MQ=255)
                                                                                  tcGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGAc                                                           >  4:402654/1‑149 (MQ=255)
                                                                                   cGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCAGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACt                                                          >  1:140728/1‑149 (MQ=255)
                                                                                     cAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTc                                                        <  3:76418/149‑1 (MQ=255)
                                                                                      aaGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTcc                                                       <  8:7196/149‑1 (MQ=255)
                                                                                       agagGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCa                                                      <  8:487430/149‑1 (MQ=255)
                                                                                       agagGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCa                                                      <  8:294773/149‑1 (MQ=255)
                                                                                        gagGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACt                                                          <  3:213352/144‑1 (MQ=255)
                                                                                        gagGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACt                                                          >  4:213352/1‑144 (MQ=255)
                                                                                        gagGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCaa                                                     <  6:159872/149‑1 (MQ=255)
                                                                                             acgacgTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATAtggtg                                                                                                                      <  3:394374/79‑1 (MQ=255)
                                                                                             acgacgTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATAtggtg                                                                                                                      >  4:394374/1‑79 (MQ=255)
                                                                                               gacgTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCg                                                >  2:418638/1‑147 (MQ=255)
                                                                                                  gTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGa                                                                                      >  2:191654/1‑106 (MQ=255)
                                                                                                   tttATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAg                                          >  1:103689/1‑149 (MQ=255)
                                                                                                      aTTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATAc                                                  >  4:164093/1‑138 (MQ=255)
                                                                                                      aTTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATAc                                                  <  3:164093/138‑1 (MQ=255)
                                                                                                      aTTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAgcg                                        >  8:345832/1‑148 (MQ=255)
                                                                                                       ttCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAgcgc                                       <  7:459342/148‑1 (MQ=255)
                                                                                                          aGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCAtct                                   >  7:386940/1‑149 (MQ=255)
                                                                                                          aGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCAtct                                   >  8:311099/1‑149 (MQ=255)
                                                                                                             gAAGCCGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCAtcttct                                <  7:345832/149‑1 (MQ=255)
                                                                                                               aGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGt                              <  6:159814/149‑1 (MQ=255)
                                                                                                                gCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGtt                             <  4:42592/149‑1 (MQ=255)
                                                                                                                gCCGGGCTACCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGtt                             <  7:192172/149‑1 (MQ=255)
                                                                                                                  cGGGCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGc                           <  2:324902/149‑1 (MQ=255)
                                                                                                                   gggCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGGAATCAGAAGCAACTACGCGAACCGCCGGCGTTGCCCCGCCAGCTGACGTCCAGTACCGCAGAAGCGCATCTACTGTTGCa                          >  2:408613/1‑149 (MQ=255)
                                                                                                                    ggCTGCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCaa                         <  1:132931/149‑1 (MQ=255)
                                                                                                                        gCCGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAgg                     >  4:54734/1‑149 (MQ=255)
                                                                                                                         ccGACCTACCGTACCCCGGAAAGAACCATCACTGCTTTTCTGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGc                    <  1:408613/149‑1 (MQ=255)
                                                                                                                            aCCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTc                                                        <  1:68754/110‑1 (MQ=255)
                                                                                                                              cTACCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTg               >  2:250812/1‑149 (MQ=255)
                                                                                                                                aCCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGcc             >  5:204809/1‑149 (MQ=255)
                                                                                                                                aCCGTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGcc             >  8:227377/1‑149 (MQ=255)
                                                                                                                                 ccGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCAAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGcc             >  8:350674/1‑148 (MQ=255)
                                                                                                                                   gTACCCCGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCAtcttct                                >  3:460643/1‑127 (MQ=255)
                                                                                                                                    tACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAg         <  3:381854/149‑1 (MQ=255)
                                                                                                                                       cccGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGTAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGaa          >  6:278979/1‑145 (MQ=255)
                                                                                                                                       cccGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGTAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGaa          <  5:278979/145‑1 (MQ=255)
                                                                                                                                       cccGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAgggg      <  7:311099/149‑1 (MQ=255)
                                                                                                                                        ccGGAAAGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGc     <  6:324740/149‑1 (MQ=255)
                                                                                                                                         cGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTCCCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGAtt                <  1:333867/137‑1 (MQ=255)
                                                                                                                                           gaaggaaCCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTAc  <  8:386940/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCAAGTATTAATTGAAGCGGTAAAGCATCATCCCCGCAGCAAATATGTCTCTTTGCTGACGAACTGGTGCGGCGAGCACTCCTCGCAAGAGGCACGACGTTTATTCAGCGAAGCCGGGCTGCCGACCTACCGTACCCCGGAAGGAACCATCACTGCTTTTATGCATATGGTGGAGTACCGGCGTAATCAGAAGCAACTACGCGAAACGCCGGCGTTGCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTAC  >  NC_000913/2721101‑2721388

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: