Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 211,231 T→A 11.5% L518H (CTC→CAC)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913211,2310TA11.5% 147.6 / 3.9 58L518H (CTC→CAC) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (4/3);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCGGATCAAA  >  NC_000913/211083‑211379
                                                                                                                                                    |                                                                                                                                                    
ccTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCt                                                                                                                                                      <  2:528519/149‑1 (MQ=255)
ccTGGCGAACAGCGTCACGGCGTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCt                                                                                                                                                      <  2:723416/149‑1 (MQ=255)
 cTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTc                                                                                                                                                     <  2:774926/149‑1 (MQ=255)
 cTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTAGTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTc                                                                                                                                                     <  2:204691/149‑1 (MQ=255)
  tGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGCTGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcg                                                                                                                                                    >  1:506122/1‑149 (MQ=255)
   ggCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCCCTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGTATCCCGGCGGCGATGGGAGCAAAATTCCAcga                                                                                                                                                   >  1:19864/1‑149 (MQ=255)
   ggCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAACTTCCAcga                                                                                                                                                   >  1:21886/1‑149 (MQ=255)
   ggCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcga                                                                                                                                                   <  2:8895/149‑1 (MQ=255)
   ggCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcga                                                                                                                                                   <  1:270598/149‑1 (MQ=255)
   ggCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATAGGTTTCTCGATCCGGCACAAGACACGAATTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGAAGCAAAAATCCAcga                                                                                                                                                   >  1:547577/1‑149 (MQ=255)
   ggCGAACAGTGGCACGGCTTTAACGATGCGAATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCAcga                                                                                                                                                   >  1:417276/1‑149 (MQ=255)
    gCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGGGGGGATCCCGGCGGCGCGGGTAGCAAAATTCCTcgac                                                                                                                                                  >  1:182023/1‑149 (MQ=255)
    gCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcgac                                                                                                                                                  <  1:340430/149‑1 (MQ=255)
     cGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCAGGCGGCGATGGTAGCAAAATTCCTcgacg                                                                                                                                                 >  2:788670/1‑149 (MQ=255)
      gAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcgacga                                                                                                                                                >  2:740430/1‑149 (MQ=255)
       gaccggggCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATAACGTTGTCGCGGGTAGCAAAATTCCTCGACGaa                                                                                                                                               <  1:794358/143‑1 (MQ=255)
            tGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTgg                                                                                                                                          >  2:801312/1‑149 (MQ=255)
              gCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGGCATCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCCAAATTCCTCGACGAACGTGGGa                                                                                                                                        >  1:483932/1‑149 (MQ=255)
                 cGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                                                                                                                     >  1:332000/1‑149 (MQ=255)
                    cTTGAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGTGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtag                                                                                                                                  <  2:608084/149‑1 (MQ=255)
                     tttAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtagt                                                                                                                                 <  1:623800/149‑1 (MQ=255)
                     tttAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGAACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGGCGAACGGGGGacagtcga                                                                                                                                 >  2:183520/1‑142 (MQ=255)
                        aaCGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCAATATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTaga                                                                                                                              >  2:378096/1‑149 (MQ=255)
                           gATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGAAGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                                                                                                                           >  1:788481/1‑149 (MQ=255)
                                ggATGCCGATCATACGTTCCTCGCTCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGCGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGCCGAACGTGGGATCGTAGTAGAGAAAACCg                                                                                                                      <  1:660233/149‑1 (MQ=255)
                                 gATGCCGATTATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACAGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCgg                                                                                                                     <  2:546762/149‑1 (MQ=255)
                                  ctGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGc                                                                                                                    <  1:155857/148‑1 (MQ=255)
                                      cGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACGGGGGATGGACGTGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAGAATTCGTCAGCGAACGTGGGATGGTGGTAGAGAACGCCGGCCGtt                                                                                                                >  2:1264/1‑149 (MQ=255)
                                       gATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                                                                                                               >  1:495664/1‑149 (MQ=255)
                                       gATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                                                                                                               >  1:495694/1‑149 (MQ=255)
                                          cATATGTTTCTCGATCCGGTTCAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa                                                                                                            <  1:556130/149‑1 (MQ=255)
                                               gTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTGGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc                                                                                                       >  1:507268/1‑149 (MQ=255)
                                                  tctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGt                                                                                                    <  1:132989/149‑1 (MQ=255)
                                                    tcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGttt                                                                                                  <  1:424693/149‑1 (MQ=255)
                                                       attcGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGGACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                                               <  2:71775/146‑1 (MQ=255)
                                                                  gTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTCTAACCTGCTGTTTCTCTTTATTATTgc                                                                                    >  2:58895/1‑148 (MQ=255)
                                                                   tCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAATACCGCCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGc                                                                                   >  2:364955/1‑149 (MQ=255)
                                                                            tGACACCGGGGATGGACGAGCAGGGTAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGTTCGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaa                                                                          <  1:539740/149‑1 (MQ=255)
                                                                                 ccGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCaa                                                                     >  1:660930/1‑149 (MQ=255)
                                                                                    gggATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAAGGCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                                                  >  1:423718/1‑149 (MQ=255)
                                                                                                 aGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAGAATTCCTCGGCGAACGTGGGAACGTAGTAGAGAAAACCGGCCCTTATAAACTGCTGTTTCTCTTTAGCATTGGCATCGATAAAACCAAAGCAATGGGATTATTg                                                     >  2:527774/1‑149 (MQ=255)
                                                                                                   ggCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCg                                                   >  2:411263/1‑149 (MQ=255)
                                                                                                    gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCGGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGACTTATTGCgc                                                  >  2:240091/1‑148 (MQ=255)
                                                                                                     cAATATGCGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATTACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                                                 <  1:766453/149‑1 (MQ=255)
                                                                                                     cAAGATGGGCGGGGAGGGGATTCCGGCGGCGCGGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                                                 <  1:476622/149‑1 (MQ=255)
                                                                                                      aaTATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATACCCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTgg                                                <  1:191166/149‑1 (MQ=255)
                                                                                                       atatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTggg                                               >  2:660909/1‑149 (MQ=255)
                                                                                                              cGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCCATGGGATTATTGCGTGGGTTGAggg                                        >  1:345390/1‑149 (MQ=255)
                                                                                                              cGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACgg                                        >  1:46984/1‑149 (MQ=255)
                                                                                                                     gggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCaa                                 <  1:378096/149‑1 (MQ=255)
                                                                                                                     gggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAAAGGGATTCTTGCGTGGGTTGACGGAATTCaa                                 >  2:40849/1‑149 (MQ=255)
                                                                                                                           cctgcggcgCTGGTAGCAAAATTCCGCGACGAACGTGGGATCGTAGTAGAGAAAACCTGTCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATGGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGctc                           <  2:21756/146‑1 (MQ=255)
                                                                                                                           ccgccgGCGCTGGTAGCAAAATTCCACGACGAACGTGGGATCGTAGTAGAGAAAACCTGCCCTTATACCCTGCTGTTTCTCTTTAGGATTGGCATCGATAAAACCAAAGCAAGGGGATTATGGCGTGGGTTGACGGAATTCAAACGctc                           <  2:19864/149‑1 (MQ=255)
                                                                                                                           ccGGGGGCGCTGGCAGCAAAATGCCACGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGGATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGctc                           <  2:21886/149‑1 (MQ=255)
                                                                                                                           ccGGCGGCCCTGGTAGCAAAATTCCACGTCGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGctc                           <  2:417276/149‑1 (MQ=255)
                                                                                                                                 gggcgggTAGCAAAAGGCCTCGAGGAACGGCGGGTCGTAGTATAGCAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGCTATCGATAAAACCAAAGCAAGGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGa                     <  1:269028/144‑1 (MQ=255)
                                                                                                                                  cgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAt                    <  1:411263/149‑1 (MQ=255)
                                                                                                                                       gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCaa               >  1:575178/1‑149 (MQ=255)
                                                                                                                                        taacAAAATTCCTCGACGAACGTGTGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGACTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAAc              <  2:305604/146‑1 (MQ=255)
                                                                                                                                           cAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTg           <  2:332000/149‑1 (MQ=255)
                                                                                                                                                    tCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCGGATCaaa  >  1:280987/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCGGATCAAA  >  NC_000913/211083‑211379

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