Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,508,092 C→A 5.6% E167* (GAG→TAG)  php ← phosphotriesterase homology protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,0920CA 5.6% 248.5 / 5.1 71E167* (GAG→TAG) phpphosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base A (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

AACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAG  >  NC_000913/3507945‑3508239
                                                                                                                                                   |                                                                                                                                                   
aaCTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTcc                                                                                                                                                    >  4:381334/1‑149 (MQ=255)
   tGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGc                                                                                                                                                 >  4:333319/1‑149 (MQ=255)
   tGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGc                                                                                                                                                 <  4:63991/149‑1 (MQ=255)
    gTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGcc                                                                                                                                                <  4:293845/149‑1 (MQ=255)
        ttGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCAt                                                                                                                                            >  1:532241/1‑149 (MQ=255)
        ttGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCAt                                                                                                                                            >  2:28799/1‑149 (MQ=255)
               tGGTGTCGAACTGCACGTAGGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTg                                                                                                                                     <  2:496274/149‑1 (MQ=255)
                   gtCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAacg                                                                                                                                 >  2:115699/1‑149 (MQ=255)
                    tCGAACTGCACTTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAacga                                                                                                                                >  3:490884/1‑149 (MQ=255)
                    tCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAacga                                                                                                                                >  4:523856/1‑149 (MQ=255)
                     cGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAacgac                                                                                                                               <  4:445644/149‑1 (MQ=255)
                     cGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCCTGGGCTTGTAGCAACCCCAGTTGCTCCAGCCCCATCGTGCTGAacgac                                                                                                                               >  2:297597/1‑149 (MQ=255)
                          tGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAacgacg                                                                                                                              >  4:480219/1‑145 (MQ=255)
                          tGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAacgacg                                                                                                                              <  3:480219/145‑1 (MQ=255)
                              cGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTg                                                                                                                      <  2:518497/149‑1 (MQ=255)
                                 aCGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGgag                                                                                                                   <  3:381334/149‑1 (MQ=255)
                                      ccGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCgg                                                                                                              <  3:404527/149‑1 (MQ=255)
                                      ccGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCgg                                                                                                              <  1:496602/149‑1 (MQ=255)
                                      ccGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCgg                                                                                                              <  3:187154/149‑1 (MQ=255)
                                             cgatcATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCg                                                                                                       <  1:139237/149‑1 (MQ=255)
                                              gatcATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAGAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCgg                                                                                                      >  4:388133/1‑149 (MQ=255)
                                               atcatcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGt                                                                                                     >  3:290007/1‑149 (MQ=255)
                                                  atcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTg                                                                                                  >  4:553376/1‑149 (MQ=255)
                                                  atcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTACAGCCCCATCGTGCTGAACGACGTATGCGTgg                                                                                                                     >  2:280972/1‑130 (MQ=255)
                                                  atcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTACAGCCCCATCGTGCTGAACGACGTATGCGTgg                                                                                                                     <  3:466149/130‑1 (MQ=255)
                                                  atcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTACAGCCCCATCGTGCTGAACGACGTATGCGTgg                                                                                                                     <  1:280972/130‑1 (MQ=255)
                                                  atcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTACAGCCCCATCGTGCTGAACGACGTATGCGTgg                                                                                                                     >  4:466149/1‑130 (MQ=255)
                                                   tcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATc                                                                                                                                           >  3:83180/1‑107 (MQ=255)
                                                   tcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATc                                                                                                                                           <  4:83180/107‑1 (MQ=255)
                                                   tcTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTgg                                                                                                 >  4:548319/1‑149 (MQ=255)
                                                     ttCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGtt                                                                                               <  2:26322/149‑1 (MQ=255)
                                                       cAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTAt                                                                                             <  1:252643/149‑1 (MQ=255)
                                                         aaaTGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATgc                                                                                           <  1:342810/149‑1 (MQ=255)
                                                          aaTGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATgcg                                                                                          >  1:114291/1‑149 (MQ=255)
                                                          aaTGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATgcg                                                                                          <  1:131695/149‑1 (MQ=255)
                                                                      gTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCgg                                                                              >  4:508763/1‑149 (MQ=255)
                                                                      gTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCgg                                                                              >  4:452759/1‑149 (MQ=255)
                                                                          tCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCaat                                                                          >  4:336401/1‑149 (MQ=255)
                                                                            tttCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCaataa                                                                        >  4:82314/1‑149 (MQ=255)
                                                                              tCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCaataaat                                                                      <  3:516138/149‑1 (MQ=255)
                                                                                 gATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATAcc                                                                   <  2:175560/149‑1 (MQ=255)
                                                                                 gATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCGTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCGGGTTATGCGCCAGCGCAGCGGCAATAAATAcc                                                                   >  2:81901/1‑149 (MQ=255)
                                                                                   tCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCtt                                                                 <  3:336401/149‑1 (MQ=255)
                                                                                       aGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTctct                                                             <  4:502912/149‑1 (MQ=255)
                                                                                          caCCAACGGTGCCGCGCGAAAGATCAACCCAGTGGGCTTGTAGCACCGCCAGTTGCCCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTATTcc                                                          <  2:389499/148‑1 (MQ=255)
                                                                                              aaCGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAgcg                                                      <  1:429633/149‑1 (MQ=255)
                                                                                              aaCGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAgcg                                                      >  4:493702/1‑149 (MQ=255)
                                                                                                   tGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTa                                                 <  3:82314/149‑1 (MQ=255)
                                                                                                     aCGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAAt                                               <  1:83015/149‑1 (MQ=255)
                                                                                                      cgcgcgAAAGATCAACACCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATc                                              <  1:28799/149‑1 (MQ=255)
                                                                                                       gcgcgAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCt                                             >  4:261573/1‑149 (MQ=255)
                                                                                                         gcgAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCttt                                           <  3:553376/149‑1 (MQ=255)
                                                                                                          cgAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTc                                          <  3:466811/149‑1 (MQ=255)
                                                                                                          cgAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTc                                          >  2:531441/1‑149 (MQ=255)
                                                                                                            aaaGATCAACCCCGTGGGCTTGTAGCAACGCCGGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCt                                        <  2:304288/149‑1 (MQ=255)
                                                                                                             aaGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCtt                                       >  4:90034/1‑149 (MQ=255)
                                                                                                               gATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCg                                     >  2:242090/1‑149 (MQ=255)
                                                                                                                aTCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGc                                    <  4:290007/149‑1 (MQ=255)
                                                                                                                aTCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGc                                    >  4:192675/1‑149 (MQ=255)
                                                                                                                aTCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGc                                    <  1:13506/149‑1 (MQ=255)
                                                                                                                 tCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCt                                   <  3:333319/149‑1 (MQ=255)
                                                                                                                    aCCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGt                                >  1:414368/1‑149 (MQ=255)
                                                                                                                     ccccGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTg                               <  2:532241/149‑1 (MQ=255)
                                                                                                                        cGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCg                            >  1:482177/1‑149 (MQ=255)
                                                                                                                               ttGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCcg                     >  2:528840/1‑149 (MQ=255)
                                                                                                                               ttGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCcg                     <  1:242090/149‑1 (MQ=255)
                                                                                                                                    gCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCgatg                >  1:261217/1‑149 (MQ=255)
                                                                                                                                         gCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATccc           >  2:258777/1‑149 (MQ=255)
                                                                                                                                             gTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCt       <  3:164315/149‑1 (MQ=255)
                                                                                                                                             gTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCt       <  1:128437/149‑1 (MQ=255)
                                                                                                                                                  tCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAg  >  2:120654/1‑149 (MQ=255)
                                                                                                                                                  tCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAg  >  1:399851/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACTGTTCTTGCCGATGGTGTCGAACTGCACGTACGCGCCGAGATCGATCATCTTCAAAATGTTGTCGAGGTTGTCTTTCAGATCGCAGTGACCAACGGTGACGCGCGAAAGATCAACCCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAG  >  NC_000913/3507945‑3508239

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