Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,491,274:1 +C 5.1% coding (677/1695 nt) oxc ← oxalyl CoA decarboxylase, ThDP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,491,2741.C 5.1% 355.2 / 13.3 79coding (677/1695 nt)oxcoxalyl CoA decarboxylase, ThDP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (29/46);  new base C (2/2);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCTT  >  NC_000913/2491128‑2491421
                                                                                                                                                   |                                                                                                                                                   
tGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gt                                                                                                                                                   <  4:86422/149‑1 (MQ=255)
    aaCGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAAt                                                                                                                                               <  4:574085/149‑1 (MQ=255)
           gCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgccgc                                                                                                                                        >  4:222689/1‑149 (MQ=255)
               ttGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGCGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGTGCCt                                                                                                                                    <  4:235960/149‑1 (MQ=255)
                tGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCtt                                                                                                                                   <  1:48890/149‑1 (MQ=255)
                tGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCtt                                                                                                                                   <  1:534194/149‑1 (MQ=255)
                 gCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCttt                                                                                                                                  >  4:397136/1‑149 (MQ=255)
                   cAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGc                                                                                                                                <  4:616525/149‑1 (MQ=255)
                        cAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCGATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAATCTGTTTATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAgg                                                                                                                           <  2:392245/149‑1 (MQ=255)
                         aaaCGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGa                                                                                                                          >  2:551711/1‑149 (MQ=255)
                             gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATATGCCGCGCCTTTGCCAAgg                                                                                                                           <  1:183520/144‑1 (MQ=255)
                             gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATATGCCGCGCCTTTGCCAAgg                                                                                                                           >  2:183520/1‑144 (MQ=255)
                             gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                      >  1:536842/1‑149 (MQ=255)
                             gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                           <  2:371487/128‑5 (MQ=38)
                             gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                           >  1:653908/1‑124 (MQ=38)
                             gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                           <  2:653908/128‑5 (MQ=38)
                             gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                           >  1:371487/1‑124 (MQ=38)
                                cgcgCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaattaa                                                                                                                   <  3:348647/149‑1 (MQ=255)
                                     aGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                                                                              >  4:675273/1‑149 (MQ=255)
                                            gCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGc                                                                                                       <  2:291296/149‑1 (MQ=255)
                                               gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                                                                                    <  2:536842/149‑1 (MQ=255)
                                               gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                                                                                    <  2:272730/149‑1 (MQ=255)
                                                aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                                   <  1:69956/149‑1 (MQ=255)
                                                aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                                   <  4:559678/149‑1 (MQ=255)
                                                 aaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTg                                                                                                  >  3:323069/1‑149 (MQ=255)
                                                      gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                             <  1:487732/149‑1 (MQ=255)
                                                      gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                             <  3:397136/149‑1 (MQ=255)
                                                      gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                             >  1:24366/1‑149 (MQ=255)
                                                       aTGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                                                                            <  4:323069/149‑1 (MQ=255)
                                                       aTGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCGT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                                                                            <  1:592566/149‑1 (MQ=255)
                                                       aTGCGTATCTCCAAGGATCCCTTTCGCCATAGACATTGGCAGCCATGGAATCTGAGCACTTTCAATAAATTCTCGCAGCTGTTAATCAGCTC‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAGTGGCCGTTCAGCTTTTGCTaaaa                                                                                            <  3:675273/149‑1 (MQ=255)
                                                        tGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAg                                                                                           <  2:24366/149‑1 (MQ=255)
                                                         gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                          <  1:616533/149‑1 (MQ=255)
                                                         gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                          >  1:81180/1‑149 (MQ=255)
                                                           gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                                        >  2:401282/1‑149 (MQ=255)
                                                                ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                                  <  1:262996/118‑1 (MQ=255)
                                                                ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                                  >  2:262996/1‑118 (MQ=255)
                                                                ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTg                                                                                   >  3:267921/1‑149 (MQ=255)
                                                                 tCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                                                                                  >  4:240089/1‑149 (MQ=255)
                                                                     ggATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTa                                                                              <  4:81710/149‑1 (MQ=255)
                                                                       aTCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGt                                                                            >  4:410045/1‑149 (MQ=255)
                                                                          ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                      >  2:459765/1‑104 (MQ=255)
                                                                          ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                      <  1:459765/104‑1 (MQ=255)
                                                                          ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                                         >  2:364982/1‑149 (MQ=255)
                                                                           cTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact                                                                        >  2:412561/1‑149 (MQ=255)
                                                                            tttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAGATTCCCGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGTCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg                                                                       <  2:410182/149‑1 (MQ=255)
                                                                            tttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg                                                                       >  4:76067/1‑149 (MQ=255)
                                                                             ttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactga                                                                      <  1:644350/149‑1 (MQ=255)
                                                                                  cATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCg                                                                 <  3:436281/149‑1 (MQ=255)
                                                                                   aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg                                                                <  2:516374/149‑1 (MQ=255)
                                                                                   aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg                                                                >  1:164550/1‑149 (MQ=255)
                                                                                         aTTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATgg                                                          <  4:133155/149‑1 (MQ=255)
                                                                                             gCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAGGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCaat                                                      <  3:222689/149‑1 (MQ=255)
                                                                                               aGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTc                                                                  >  1:265489/1‑135 (MQ=255)
                                                                                               aGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTc                                                                  <  2:265489/135‑1 (MQ=255)
                                                                                               aGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCaataa                                                    <  4:3958/149‑1 (MQ=255)
                                                                                                    tggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTg                                               >  1:89278/1‑149 (MQ=255)
                                                                                                     ggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTgg                                              <  2:164550/149‑1 (MQ=255)
                                                                                                        atCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGa                                           >  3:421514/1‑149 (MQ=255)
                                                                                                         tCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGAc                                          <  2:320367/149‑1 (MQ=255)
                                                                                                          cTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCCAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACg                                         <  4:126637/149‑1 (MQ=255)
                                                                                                              gCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGAtt                                     <  4:652915/149‑1 (MQ=255)
                                                                                                                 cTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTc                                  <  2:550000/149‑1 (MQ=255)
                                                                                                                   ttCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCaa                                <  1:551711/149‑1 (MQ=255)
                                                                                                                   ttCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCaa                                <  1:401282/149‑1 (MQ=255)
                                                                                                                         aaaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                          <  1:412561/149‑1 (MQ=255)
                                                                                                                         aaaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                          <  4:267921/149‑1 (MQ=255)
                                                                                                                            ttCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACaa                       <  3:440794/149‑1 (MQ=255)
                                                                                                                              cccGAAGCGGTTCATCAGCTT‑GTGAATACGCCTCCCCTGTGCGAAGGATAATTACTGCCCGTTCAGCGTGTGCTACAAGCGAAATTTGGATCGTGACTAGCTTCGCACATGCCACTAATAATGTCGCCGTTTTCTCATCTTCAACACTc                     <  2:249710/147‑1 (MQ=255)
                                                                                                                               aCGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCg                    >  4:87825/1‑149 (MQ=255)
                                                                                                                                 gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                                         >  2:442993/1‑94 (MQ=255)
                                                                                                                                 gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                                         <  1:442993/94‑1 (MQ=255)
                                                                                                                                 gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg                  >  4:49343/1‑149 (MQ=255)
                                                                                                                                 gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg                  <  3:76067/149‑1 (MQ=255)
                                                                                                                                 gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg                  >  3:300371/1‑149 (MQ=255)
                                                                                                                                   aGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGt                <  3:240089/149‑1 (MQ=255)
                                                                                                                                       gTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAc            >  3:521084/1‑149 (MQ=255)
                                                                                                                                            tcaGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTc       >  1:566834/1‑149 (MQ=255)
                                                                                                                                              aGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGt     >  3:571212/1‑149 (MQ=255)
                                                                                                                                               gCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTc    <  2:81180/149‑1 (MQ=255)
                                                                                                                                                 tt‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCtt  <  2:566834/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCTT  >  NC_000913/2491128‑2491421

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