Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,491,266 C→G 5.3% E229Q (GAA→CAA)  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,2660CG 5.3% 261.3 / 5.7 78E229Q (GAA→CAA) oxcoxalyl CoA decarboxylase, ThDP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (28/44);  new base G (2/2);  total (30/48)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CCAACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAC  >  NC_000913/2491118‑2491411
                                                                                                                                                    |                                                                                                                                                  
ccAACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtc                                                                                                                                                    <  1:106148/149‑1 (MQ=255)
   gcaaGCATGAAAACGTCATCATTTGTCAGGGCAAAAGAACGCGCAGCTGCCGCAGAAAGTGGATGGGTATCTTCAAGGCTCCTACCTGACAGAGAACTTTTCCGGAATGGAATCTGATCCCTTTCAATAAATTCACGGACAGGTtcatc                                                                                                                                                 <  2:622996/148‑1 (MQ=255)
   acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCTCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                 >  2:592566/1‑149 (MQ=255)
   acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                 >  4:440794/1‑149 (MQ=255)
   acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                 >  3:652421/1‑149 (MQ=255)
   acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                 <  1:286521/149‑1 (MQ=255)
   acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                 <  4:652421/149‑1 (MQ=255)
       gCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCt                                                                                                                                             >  1:314765/1‑149 (MQ=255)
          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                                                                                                                                 >  1:371487/1‑124 (MQ=38)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                 <  2:371487/128‑5 (MQ=38)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                 <  2:653908/128‑5 (MQ=38)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                 >  1:653908/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:272730/149‑1 (MQ=255)
                                                         gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                                                                          <  2:536842/149‑1 (MQ=255)
                                                          aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                         <  4:559678/149‑1 (MQ=255)
                                                          aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                         <  1:69956/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:81180/1‑149 (MQ=255)
                                                                   gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                <  1:616533/149‑1 (MQ=255)
                                                                     gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                              >  2:401282/1‑149 (MQ=255)
                                                                          ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                        >  2:262996/1‑118 (MQ=255)
                                                                          ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                        <  1:262996/118‑1 (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                                                      >  1:164550/1‑149 (MQ=255)
                                                                                             aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg                                                      <  2:516374/149‑1 (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:401282/149‑1 (MQ=255)
                                                                                                                             ttCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCaa                      <  1:551711/149‑1 (MQ=255)
                                                                                                                                   aaaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                <  4:267921/149‑1 (MQ=255)
                                                                                                                                   aaaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                <  1:412561/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                                                               <  1:442993/94‑1 (MQ=255)
                                                                                                                                           gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                               >  2:442993/1‑94 (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)
                                                                                                                                                    |                                                                                                                                                  
CCAACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAC  >  NC_000913/2491118‑2491411

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