Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,720,495 C→T 5.3% G181G (GGC→GGT pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,720,4950CT 5.3% 241.3 / 2.4 75G181G (GGC→GGTpkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base T (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATT  >  NC_000913/2720347‑2720642
                                                                                                                                                    |                                                                                                                                                   
ccTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGc                                                                                                                                                     >  3:109051/1‑149 (MQ=255)
ccTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGc                                                                                                                                                     <  3:71309/149‑1 (MQ=255)
 cTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGc                                                                                                                                                     >  4:182042/1‑148 (MQ=255)
 cTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACCGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGc                                                                                                                                                     >  8:35351/1‑148 (MQ=255)
       aGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGcc                                                                                                                                              <  2:261787/149‑1 (MQ=255)
        ggTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCt                                                                                                                                             >  2:302509/1‑149 (MQ=255)
         gTCTGAAGGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTg                                                                                                                                            <  8:196181/149‑1 (MQ=255)
          tCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTgg                                                                                                                                           <  4:109051/149‑1 (MQ=255)
            tGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGa                                                                                                                                          >  8:198479/1‑148 (MQ=255)
              aaTGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGata                                                                                                                                        <  4:305271/148‑1 (MQ=255)
               aTGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATc                                                                                                                                      <  2:212337/149‑1 (MQ=255)
                tGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATc                                                                                                                                      >  7:353718/1‑148 (MQ=255)
                   cAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGTGACAGCCTGGATATCg                                                                                                                                     <  5:147726/146‑1 (MQ=255)
                   cAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGTGACAGCCTGGATATCg                                                                                                                                     >  6:147726/1‑146 (MQ=255)
                       ttttCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTtga                                                                                                                              >  5:202486/1‑149 (MQ=255)
                         ttCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCACTCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTtgat                                                                                                                             <  2:94952/148‑1 (MQ=255)
                         ttCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTtgat                                                                                                                             <  8:353718/148‑1 (MQ=255)
                         ttCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTtga                                                                                                                              <  5:168801/147‑1 (MQ=255)
                             ccTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAAtt                                                                                                                        <  6:202486/149‑1 (MQ=255)
                                 ttccGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGc                                                                                                                      <  1:101286/145‑1 (MQ=255)
                                 tgCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCtt                                                                                                                    <  6:340017/149‑1 (MQ=255)
                                 tgCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCtt                                                                                                                    <  1:302509/149‑1 (MQ=255)
                                   ccGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGa                                                                                                                  <  4:336261/149‑1 (MQ=255)
                                     gATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTg                                                                                                                                            >  7:237121/1‑121 (MQ=255)
                                     gATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTg                                                                                                                                            <  8:237121/121‑1 (MQ=255)
                                          aaCGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTg                                                                                                           >  2:137284/1‑149 (MQ=255)
                                           aCGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTgg                                                                                                          <  1:137284/149‑1 (MQ=255)
                                                    gCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGa                                                                                                   >  6:5450/1‑147 (MQ=255)
                                                     cTGGCGTTTATTTCGCAATCGGCGGCCGTCTCCAACACCATGCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACa                                                                                                 <  8:99485/148‑1 (MQ=255)
                                                      tGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAAtt                                                                                                                        >  2:247381/1‑124 (MQ=255)
                                                      tGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAAtt                                                                                                                        <  1:247381/124‑1 (MQ=255)
                                                          gTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCTACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTTCTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTAGCGTATCTGGCACGCGACAGTaaca                                                                                           >  4:67228/1‑147 (MQ=255)
                                                               tttCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAgc                                                                                      <  5:5450/149‑1 (MQ=255)
                                                                    cAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTATCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGAGGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCAt                                                                                 >  8:63631/1‑149 (MQ=255)
                                                                     aaTCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGTCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATc                                                                                <  3:182042/149‑1 (MQ=255)
                                                                     aaTCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATc                                                                                >  7:211469/1‑149 (MQ=255)
                                                                            tGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAg                                                                                       >  1:208632/1‑135 (MQ=255)
                                                                            tGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAg                                                                                       <  2:208632/135‑1 (MQ=255)
                                                                              ccGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCTCCATCCTGCTCTAt                                                                       >  8:125379/1‑149 (MQ=255)
                                                                               cGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTAtc                                                                      <  8:211469/149‑1 (MQ=255)
                                                                                    ccAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGaa                                                                 >  2:276207/1‑149 (MQ=255)
                                                                                     cAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGaa                                                                 >  1:53443/1‑148 (MQ=255)
                                                                                     cAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGaa                                                                 <  8:84577/148‑1 (MQ=255)
                                                                                      aaCACCATCCTCGACCGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACa                                                               <  7:35351/149‑1 (MQ=255)
                                                                                           cATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGtt                                                            >  3:107409/1‑147 (MQ=255)
                                                                                           cATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGtt                                                            >  3:3092/1‑147 (MQ=255)
                                                                                             tCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGGTTGACTATCTGGCACGCGACa                                                                                                 <  1:22255/108‑1 (MQ=255)
                                                                                             tCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGGTTGACTATCTGGCACGCGACa                                                                                                 >  2:22255/1‑108 (MQ=255)
                                                                                               cTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATcc                                                                               >  4:58468/1‑124 (MQ=255)
                                                                                                tCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAgc                                                                                      <  3:120719/116‑1 (MQ=255)
                                                                                                tCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAgc                                                                                      >  4:120719/1‑116 (MQ=255)
                                                                                                tCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTg                                                                             <  2:88896/125‑1 (MQ=255)
                                                                                                tCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTg                                                                             >  1:88896/1‑125 (MQ=255)
                                                                                                    cTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTAACTGGCACGCGACAGTAAAACCAGCCCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgc                                                 >  2:78163/1‑149 (MQ=255)
                                                                                                     tGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGa                                                                                                   >  1:182534/1‑98 (MQ=255)
                                                                                                     tGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGa                                                                                                   <  2:182534/98‑1 (MQ=255)
                                                                                                            aaCAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCttt                                         <  4:3092/149‑1 (MQ=255)
                                                                                                            aaCAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCttt                                         >  1:231605/1‑149 (MQ=255)
                                                                                                                 cGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTc                                    <  4:107409/149‑1 (MQ=255)
                                                                                                                  gTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcg                                   <  2:53443/149‑1 (MQ=255)
                                                                                                                  gTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcg                                   >  1:153809/1‑149 (MQ=255)
                                                                                                                  gTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcg                                   >  1:104627/1‑149 (MQ=255)
                                                                                                                   tAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcgg                                  <  5:163338/149‑1 (MQ=255)
                                                                                                                       aTGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCggcggc                              >  6:270236/1‑149 (MQ=255)
                                                                                                                               ttCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGt                       <  1:276207/148‑1 (MQ=255)
                                                                                                                                  cTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCttt                                         >  8:342860/1‑127 (MQ=255)
                                                                                                                                  cTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCttt                                         <  7:342860/127‑1 (MQ=255)
                                                                                                                                   tACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTc                  >  4:31474/1‑149 (MQ=255)
                                                                                                                                     cTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTcgc                <  2:104627/149‑1 (MQ=255)
                                                                                                                                      tttATTGCGCTCGGCGACAGCCTGGATATCGACGTAGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAAGAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGCCCCGTAGTGCCTcgcg               >  1:315836/1‑149 (MQ=255)
                                                                                                                                         aTTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgc                                                   >  2:94616/1‑110 (MQ=255)
                                                                                                                                         aTTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgc                                                   <  1:94616/110‑1 (MQ=255)
                                                                                                                                             cgcTCTGCGACAGGCTGGATATCGACGTTGATGAATTGCTTGAGTATCTTGCTCGCGACAGTAAAACCTGCGCCATCCTGCTCTATCTCTAACAGTTAAGCGACGCGCGTCGCTTTGTTTCGGCGTCCCGTAGTGCCTCGCGTAATaaa        <  2:315836/149‑1 (MQ=255)
                                                                                                                                              gcTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAAc       <  2:231605/149‑1 (MQ=255)
                                                                                                                                                tCGGTGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACTGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATaaa        <  2:142931/146‑1 (MQ=255)
                                                                                                                                                tCGGTGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACGGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATaaa        >  1:142931/1‑146 (MQ=255)
                                                                                                                                                  ggCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTc                  >  4:163668/1‑134 (MQ=255)
                                                                                                                                                  ggCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTc                  <  3:163668/134‑1 (MQ=255)
                                                                                                                                                   gCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGTCAGTAAGACCTGCGCGATCCTGCTCTATCTCGGAGTGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCGCGTAGTGCCTCGCGTAATAAACGGagt  >  2:253504/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGATT  >  NC_000913/2720347‑2720642

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: