Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,722,202 T→G D750E (GAT→GAG pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,722,2020TG97.3% 259.1 / ‑2.4 74D750E (GAT→GAGpkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (41/31);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGC  >  NC_000913/2722056‑2722344
                                                                                                                                                  |                                                                                                                                              
gagaGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCt                                                                                                                                              >  7:235989/1‑149 (MQ=255)
gagaGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCt                                                                                                                                              >  8:61340/1‑149 (MQ=255)
       ggCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTttatta                                                                                                                                        <  8:222435/148‑1 (MQ=255)
        gCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTAcc                                                                                                                                      >  4:286567/1‑149 (MQ=255)
           ggCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGc                                                                                                                                    >  3:227136/1‑148 (MQ=255)
                  cgcgCTTATCCGCATCTGCTGTAAGATTGGGTATAATTGAATATCGGTGAACGCTGCTTGTTCCGCCCTATTTTGCCATAAGATTAGCCTGAACTTCAGCAATTCATTTAGCGAGTCACCAAAGAAGAGCTTTATTTCCGCTACTTTa                                                                                                                             <  5:54161/148‑1 (MQ=255)
                   gcCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGc                                                                                                                           <  6:208329/149‑1 (MQ=255)
                      cTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCgag                                                                                                                        <  1:257576/149‑1 (MQ=255)
                         aTCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGAt                                                                                                                      <  5:17397/148‑1 (MQ=255)
                            cGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAAc                                                                                                                  <  8:235989/149‑1 (MQ=255)
                            cGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAAc                                                                                                                  >  1:166836/1‑149 (MQ=255)
                            cGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAAc                                                                                                                  <  8:107947/149‑1 (MQ=255)
                              cATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACg                                                                                                                 >  1:325038/1‑148 (MQ=255)
                               aTCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGaa                                                                                                               >  3:359719/1‑149 (MQ=255)
                                  aGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTTCCAGAAGATGATCCACAACTTCAGCAAATCATTTCGCGAGTCACCAAAGAAGAGATTTATTAACTCTACTTTAGCGAGATCAACTAAtt                                                                                                             >  4:218250/1‑148 (MQ=255)
                                          gaaTGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACACATg                                                                                                     >  5:327078/1‑148 (MQ=255)
                                           aaTGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTa                                                                                                                                  >  7:405355/1‑118 (MQ=255)
                                           aaTGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTa                                                                                                                                  <  8:405355/118‑1 (MQ=255)
                                            aTGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCtt                                                                                                                                             >  3:10453/1‑106 (MQ=255)
                                            aTGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCtt                                                                                                                                             <  4:10453/106‑1 (MQ=255)
                                               ggTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGAtt                                                                                               <  4:343098/149‑1 (MQ=255)
                                                gTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGAtt                                                                                               >  2:215179/1‑148 (MQ=255)
                                                gTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGAtt                                                                                               >  4:336564/1‑148 (MQ=255)
                                                    aaTTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGcc                                                                                          <  1:176568/149‑1 (MQ=255)
                                                     aTTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGcc                                                                                          >  8:28975/1‑148 (MQ=255)
                                                      ttGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCaa                                                                                        >  7:407758/1‑149 (MQ=255)
                                                        gAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTTCCAGAAGATTAGCCACATCTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAAc                                                                                       <  7:207597/148‑1 (MQ=255)
                                                        gAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACa                                                                                      >  5:259940/1‑149 (MQ=255)
                                                           aaaCGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGa                                                                                   >  3:336553/1‑149 (MQ=255)
                                                             aCGGTGCACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGaca                                                                                 <  2:295832/149‑1 (MQ=255)
                                                             aCGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGaca                                                                                 <  4:359719/149‑1 (MQ=255)
                                                                gTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGaca                                                                                 >  5:19926/1‑146 (MQ=255)
                                                                 tGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACCAATTTACCCATTAAGATTTAGCCAACATGaca                                                                                 >  7:397299/1‑145 (MQ=255)
                                                                   aaCGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGaca                                                                                 >  2:4569/1‑143 (MQ=255)
                                                                   aaCGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAt                                                                            <  8:13285/148‑1 (MQ=255)
                                                                    aCGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCg                                                                          <  7:61340/149‑1 (MQ=255)
                                                                          ctTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACCAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTAc                                                                     >  6:391027/1‑148 (MQ=255)
                                                                                 cGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTGCCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCggg                                                              >  1:279640/1‑148 (MQ=255)
                                                                                 cGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCggg                                                              >  4:241341/1‑148 (MQ=255)
                                                                                  gCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGa                                                             >  7:406251/1‑148 (MQ=255)
                                                                                   cccGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGaaa                                                           >  6:146850/1‑149 (MQ=255)
                                                                                     ctaTTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGAGCGGGAAATg                                                         <  5:186355/147‑1 (MQ=255)
                                                                                      gATTTTGACAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATgg                                                        <  3:218250/149‑1 (MQ=255)
                                                                                        ttttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCg                                                      <  5:208259/149‑1 (MQ=255)
                                                                                          ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGt                                                     <  7:76350/148‑1 (MQ=255)
                                                                                           tGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGtt                                                    <  6:190061/148‑1 (MQ=255)
                                                                                            gCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGttt                                                   <  1:4569/148‑1 (MQ=255)
                                                                                              cAGAAGATGAGCCACAACTTCAGCAATTGATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTg                                                  >  5:11480/1‑147 (MQ=255)
                                                                                              cAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTg                                                  >  4:361010/1‑147 (MQ=255)
                                                                                              cAGAAGATGAGACACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGATTTATTACCGCTACTTTAGAGAGATCAACGAATTTACCCATGTAGATTTAGCCAACATGACAAAGATCGACTACGATCGGGAAATGGAGTTTGTa                                                <  1:368614/149‑1 (MQ=255)
                                                                                                 aagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCg                                             >  1:9758/1‑149 (MQ=255)
                                                                                                  agaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCaa                                                                                        <  1:366963/105‑1 (MQ=255)
                                                                                                  agaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCaa                                                                                        >  2:366963/1‑105 (MQ=255)
                                                                                                   gaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCgg                                            >  1:222349/1‑148 (MQ=255)
                                                                                                        gCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTacgac                                      <  4:227136/149‑1 (MQ=255)
                                                                                                         ccACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTacgacg                                     >  4:146058/1‑149 (MQ=255)
                                                                                                           caaaCTTAAGCAATTAATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATAAACGAATTTACCCAAGAAGATTTAGCCAACATGACACAGATCGACTAAGATCGCGAAATGCCGTTTGTAGCGGTACGACGTa                                   >  6:154830/3‑149 (MQ=255)
                                                                                                           acaACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGAGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCCGTACGACGTa                                   >  5:55247/1‑149 (MQ=255)
                                                                                                                  cAGCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACTAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                              <  5:165769/147‑1 (MQ=255)
                                                                                                                    gCAATTCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCaaa                          <  6:259940/149‑1 (MQ=255)
                                                                                                                        ttCATTTCGCGAGTCACCAAAGAAGAGCTTTATTACAGCTACTTTAGAGAGATCAACGAATTTACCCATGAAGATTTAGCAAAAATGACACAGAACGACTACGATCGGGAAATGGCGATTGTAGCGGTACGACGTATTGATCAAACgg                       >  1:69301/1‑148 (MQ=255)
                                                                                                                           aTTTCGCGAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGaa                     >  7:153353/1‑147 (MQ=255)
                                                                                                                                gcgAGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATc               <  6:19926/148‑1 (MQ=255)
                                                                                                                                   aGTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCg           <  3:286567/149‑1 (MQ=255)
                                                                                                                                    gTCACCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTTGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCgg          >  7:79759/1‑149 (MQ=255)
                                                                                                                                       aCCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                              >  4:352323/1‑126 (MQ=255)
                                                                                                                                       aCCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                              <  3:352323/126‑1 (MQ=255)
                                                                                                                                       aCCAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt       >  5:64224/1‑149 (MQ=255)
                                                                                                                                         cAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAAt                                                          <  4:332294/96‑1 (MQ=255)
                                                                                                                                         cAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAAt                                                          >  3:332294/1‑96 (MQ=255)
                                                                                                                                         cAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGACGAGATCCTCGGCGTCa     >  2:193959/1‑149 (MQ=255)
                                                                                                                                         cAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa     >  3:264688/1‑149 (MQ=255)
                                                                                                                                         cAAAGAAGAGCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa     <  5:146850/149‑1 (MQ=255)
                                                                                                                                            ataagaGCTATATTACCGCTACTTTAGCGAGATCAACGAATTTACCAATGAAGATTTAGACAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcgc  >  6:376893/3‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGTGAACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGC  >  NC_000913/2722056‑2722344

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