Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,175,546 T→A 5.0% Q154L (CAG→CTG)  yqiA ← acyl CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,175,5460TA 5.0% 280.1 / 4.9 80Q154L (CAG→CTG) yqiAacyl CoA esterase
Reads supporting (aligned to +/- strand):  ref base T (35/41);  new base A (2/2);  total (37/43)
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.27e-01

TTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCT  >  NC_000913/3175399‑3175693
                                                                                                                                                   |                                                                                                                                                   
tttAAGTTAGTAATTCGAGTTGGTCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTg                                                                                                                                                    >  1:379603/1‑149 (MQ=255)
tttAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTg                                                                                                                                                    <  3:494581/149‑1 (MQ=255)
       tAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTaa                                                                                                                                             >  3:191284/1‑149 (MQ=255)
        aGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAAt                                                                                                                                            >  4:86797/1‑149 (MQ=255)
         gTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATc                                                                                                                                           >  1:51046/1‑149 (MQ=255)
           aaTTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCa                                                                                                                                         <  2:80950/149‑1 (MQ=255)
                   ttGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATATTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCa                                                                                                                                 >  4:241098/1‑149 (MQ=255)
                    tGATCGTCAGAGATGGTGCAGACCAAGTAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCAt                                                                                                                                <  3:432960/149‑1 (MQ=255)
                    tGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCAt                                                                                                                                <  2:592902/149‑1 (MQ=255)
                      aTCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCAtct                                                                                                                              <  1:484750/149‑1 (MQ=255)
                      aTCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCAtct                                                                                                                              >  4:378299/1‑149 (MQ=255)
                         gTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTccc                                                                                                                           >  1:285049/1‑149 (MQ=255)
                         gTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCAGGCGGtgg                                                                                                                                             >  1:524332/1‑129 (MQ=255)
                         gTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCAGGCGGtgg                                                                                                                                             <  2:524332/131‑3 (MQ=255)
                          tCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCg                                                                                                                          >  3:245771/1‑149 (MQ=255)
                                     gCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTgcagc                                                                                                               >  2:120329/1‑149 (MQ=255)
                                      cAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTgcagca                                                                                                              <  4:195628/149‑1 (MQ=255)
                                             aGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGAt                                                                                                       <  4:308930/149‑1 (MQ=255)
                                              gAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATc                                                                                                      >  2:137057/1‑149 (MQ=255)
                                                aaaTCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCgg                                                                                                                                                >  3:558698/1‑105 (MQ=255)
                                                aaaTCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCgg                                                                                                                                                <  4:558698/105‑1 (MQ=255)
                                                aaaTCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCaa                                                                                                    >  3:82864/1‑149 (MQ=255)
                                                 aaTCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCaaa                                                                                                   <  4:213637/149‑1 (MQ=255)
                                                 aaTCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCaaa                                                                                                   >  2:405116/1‑149 (MQ=255)
                                                  aTCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAAt                                                                                                  >  1:270917/1‑149 (MQ=255)
                                                     gacgaTCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCa                                                                                                                                         <  4:350942/107‑1 (MQ=255)
                                                     gacgaTCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCa                                                                                                                                         >  3:350942/1‑107 (MQ=255)
                                                      acgaTCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCgg                                                                                              <  2:51046/149‑1 (MQ=255)
                                                         aTCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGgcgc                                                                                           <  2:20293/149‑1 (MQ=255)
                                                          tCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCt                                                                                          >  2:585381/1‑149 (MQ=255)
                                                                tGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGc                                                                                    >  4:53968/1‑149 (MQ=255)
                                                                 gAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCg                                                                                   <  4:191284/149‑1 (MQ=255)
                                                                        tCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAAt                                                                            >  4:413372/1‑149 (MQ=255)
                                                                          ttCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCt                                                                          >  3:449061/1‑149 (MQ=255)
                                                                          ttCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCt                                                                          >  1:623087/1‑149 (MQ=255)
                                                                              aaGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCAt                                                                      >  2:157824/1‑149 (MQ=255)
                                                                              aaGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCAt                                                                      <  2:76483/149‑1 (MQ=255)
                                                                                  ccGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCgg                                                                                  >  2:620942/1‑133 (MQ=255)
                                                                                  ccGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCgg                                                                                  <  1:620942/133‑1 (MQ=255)
                                                                                  ccGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACt                                                                  >  1:620477/1‑149 (MQ=255)
                                                                                      gAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTaa                                                              <  3:285286/149‑1 (MQ=255)
                                                                                         tGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGAt                                                           >  3:295381/1‑149 (MQ=255)
                                                                                         tGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGAt                                                           <  1:35761/149‑1 (MQ=255)
                                                                                            gTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGt                                                        <  4:278917/149‑1 (MQ=255)
                                                                                            gTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGt                                                        >  4:491999/1‑149 (MQ=255)
                                                                                                  gccgccTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAatat                                                  <  1:405116/149‑1 (MQ=255)
                                                                                                   ccgccTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATg                                                 <  3:241098/149‑1 (MQ=255)
                                                                                                     gccTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGc                                               <  3:86797/149‑1 (MQ=255)
                                                                                                     gccTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGc                                               <  2:623087/149‑1 (MQ=255)
                                                                                                     gccTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGc                                               >  2:431396/1‑149 (MQ=255)
                                                                                                      ccTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCg                                              <  4:245771/149‑1 (MQ=255)
                                                                                                       cTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGt                                             >  3:242759/1‑149 (MQ=255)
                                                                                                       cTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGt                                             <  3:219679/149‑1 (MQ=255)
                                                                                                        ttCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGgcg                                                                                            <  4:154676/101‑1 (MQ=255)
                                                                                                        tccTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACAGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGgcg                                                                                            >  3:154676/3‑101 (MQ=255)
                                                                                                         tCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGa                                           <  2:379603/149‑1 (MQ=255)
                                                                                                            aTGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGActc                                        >  1:574539/1‑149 (MQ=255)
                                                                                                            aTGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGActc                                        <  3:178700/149‑1 (MQ=255)
                                                                                                             tGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGActct                                       <  3:495916/149‑1 (MQ=255)
                                                                                                               aCAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCAtc                                                                                                                               >  1:78473/1‑59 (MQ=255)
                                                                                                               aCAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCAtc                                                                                                                               <  2:78473/59‑1 (MQ=255)
                                                                                                               aCAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAg                                     >  3:364403/1‑149 (MQ=255)
                                                                                                                 aGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGca                                   <  4:208024/149‑1 (MQ=255)
                                                                                                                   tCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTccc                                                                                                                           <  3:276523/59‑1 (MQ=255)
                                                                                                                   tCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTccc                                                                                                                           >  4:276523/1‑59 (MQ=255)
                                                                                                                   tCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGcaca                                 <  3:491999/149‑1 (MQ=255)
                                                                                                                   tCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGcaca                                 <  2:285049/149‑1 (MQ=255)
                                                                                                                     tGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCCATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACata                               <  4:295381/149‑1 (MQ=255)
                                                                                                                            aGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATtgctgc                        <  2:598440/149‑1 (MQ=255)
                                                                                                                                   gtagtaCGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGCCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGgtgt                 >  3:37819/1‑149 (MQ=255)
                                                                                                                                    tagtaCGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTa                <  4:332401/149‑1 (MQ=255)
                                                                                                                                    tagtaCGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTa                <  2:184119/149‑1 (MQ=255)
                                                                                                                                     agtaCGCCACCGCCTGGCGGTTATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAg               <  4:449061/149‑1 (MQ=255)
                                                                                                                                     agtaCGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAg               <  1:157824/149‑1 (MQ=255)
                                                                                                                                      gtaCGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCTGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAgg              >  3:60179/1‑149 (MQ=255)
                                                                                                                                         gattaCCGCCAGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATc                                                                           <  3:604486/81‑1 (MQ=255)
                                                                                                                                         gattaCCGCCAGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATc                                                                           >  4:604486/5‑85 (MQ=255)
                                                                                                                                          gCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAGGGGtt          <  1:120329/149‑1 (MQ=255)
                                                                                                                                             aCCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAGGGGTtctc       >  4:574156/1‑149 (MQ=255)
                                                                                                                                                  cTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCt  <  1:137057/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGACGATCGGGTTGAAATAATCTTCGAAGCCCGTGAATGCGTGGTTGCCGCCTTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCT  >  NC_000913/3175399‑3175693

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: