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

CCAACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAC  >  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)
          tGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTtgt                                                                                                                                         <  4:86422/149‑1 (MQ=255)
              aaCGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAAt                                                                                                                                     <  4:574085/149‑1 (MQ=255)
                     gCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATAcgccgc                                                                                                                              >  4:222689/1‑149 (MQ=255)
                         ttGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGCGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGTGCCt                                                                                                                          <  4:235960/149‑1 (MQ=255)
                          tGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCtt                                                                                                                         <  1:48890/149‑1 (MQ=255)
                          tGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCtt                                                                                                                         <  1:534194/149‑1 (MQ=255)
                           gCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCttt                                                                                                                        >  4:397136/1‑149 (MQ=255)
                             cAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGc                                                                                                                      <  4:616525/149‑1 (MQ=255)
                                  cAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCGATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAATCTGTTTATCAGCTTGTGAATACGCCGCGCCTTTGCCAAgg                                                                                                                 <  2:392245/149‑1 (MQ=255)
                                   aaaCGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGa                                                                                                                >  2:551711/1‑149 (MQ=255)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATATGCCGCGCCTTTGCCAAgg                                                                                                                 <  1:183520/144‑1 (MQ=255)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATATGCCGCGCCTTTGCCAAgg                                                                                                                 >  2:183520/1‑144 (MQ=255)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                            >  1:536842/1‑149 (MQ=255)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                 >  1:371487/1‑124 (MQ=38)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                 <  2:371487/128‑5 (MQ=38)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                 <  2:653908/128‑5 (MQ=38)
                                       gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                 >  1:653908/1‑124 (MQ=38)
                                          cgcgCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaa                                                                                                         <  3:348647/149‑1 (MQ=255)
                                               aGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                                                                    >  4:675273/1‑149 (MQ=255)
                                                      gCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGc                                                                                             <  2:291296/149‑1 (MQ=255)
                                                         gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                                                                          <  2:272730/149‑1 (MQ=255)
                                                         gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                                                                          <  2:536842/149‑1 (MQ=255)
                                                          aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                         <  4:559678/149‑1 (MQ=255)
                                                          aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                         <  1:69956/149‑1 (MQ=255)
                                                           aaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTg                                                                                        >  3:323069/1‑149 (MQ=255)
                                                                gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                   <  1:487732/149‑1 (MQ=255)
                                                                gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                   <  3:397136/149‑1 (MQ=255)
                                                                gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                   >  1:24366/1‑149 (MQ=255)
                                                                 aTGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                                                                  <  4:323069/149‑1 (MQ=255)
                                                                 aTGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCGTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                                                                  <  1:592566/149‑1 (MQ=255)
                                                                 aTGCGTATCTCCAAGGATCCCTTTCGCCATAGACATTGGCAGCCATGGAATCTGAGCACTTTCAATAAATTCTCGCAGCTGTTAATCAGCTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAGTGGCCGTTCAGCTTTTGCTaaaa                                                                                  <  3:675273/149‑1 (MQ=255)
                                                                  tGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAg                                                                                 <  2:24366/149‑1 (MQ=255)
                                                                   gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                >  1:81180/1‑149 (MQ=255)
                                                                   gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                <  1:616533/149‑1 (MQ=255)
                                                                     gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                              >  2:401282/1‑149 (MQ=255)
                                                                          ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                        >  2:262996/1‑118 (MQ=255)
                                                                          ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                        <  1:262996/118‑1 (MQ=255)
                                                                          ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTg                                                                         >  3:267921/1‑149 (MQ=255)
                                                                           tCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                                                                        >  4:240089/1‑149 (MQ=255)
                                                                               ggATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTa                                                                    <  4:81710/149‑1 (MQ=255)
                                                                                 aTCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGt                                                                  >  4:410045/1‑149 (MQ=255)
                                                                                    ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                            >  2:459765/1‑104 (MQ=255)
                                                                                    ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                            <  1:459765/104‑1 (MQ=255)
                                                                                    ccTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                               >  2:364982/1‑149 (MQ=255)
                                                                                     cTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgact                                                              >  2:412561/1‑149 (MQ=255)
                                                                                      tttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAGATTCCCGAAGCTGTTCATCAGCTTGTGAATACGCCGCGTCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg                                                             <  2:410182/149‑1 (MQ=255)
                                                                                      tttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg                                                             >  4:76067/1‑149 (MQ=255)
                                                                                       ttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactga                                                            <  1:644350/149‑1 (MQ=255)
                                                                                            cATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCg                                                       <  3:436281/149‑1 (MQ=255)
                                                                                             aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg                                                      >  1:164550/1‑149 (MQ=255)
                                                                                             aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg                                                      <  2:516374/149‑1 (MQ=255)
                                                                                                   aTTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATgg                                                <  4:133155/149‑1 (MQ=255)
                                                                                                       gCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAGGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCaat                                            <  3:222689/149‑1 (MQ=255)
                                                                                                         aGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTc                                                        >  1:265489/1‑135 (MQ=255)
                                                                                                         aGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTc                                                        <  2:265489/135‑1 (MQ=255)
                                                                                                         aGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCaataa                                          <  4:3958/149‑1 (MQ=255)
                                                                                                              tggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTg                                     >  1:89278/1‑149 (MQ=255)
                                                                                                               ggaatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTgg                                    <  2:164550/149‑1 (MQ=255)
                                                                                                                  atCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGa                                 >  3:421514/1‑149 (MQ=255)
                                                                                                                   tCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGAc                                <  2:320367/149‑1 (MQ=255)
                                                                                                                    cTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCCAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACg                               <  4:126637/149‑1 (MQ=255)
                                                                                                                        gCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGAtt                           <  4:652915/149‑1 (MQ=255)
                                                                                                                           cTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTc                        <  2:550000/149‑1 (MQ=255)
                                                                                                                             ttCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCaa                      <  1:401282/149‑1 (MQ=255)
                                                                                                                             ttCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCaa                      <  1:551711/149‑1 (MQ=255)
                                                                                                                                   aaaTTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                <  4:267921/149‑1 (MQ=255)
                                                                                                                                   aaaTTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                <  1:412561/149‑1 (MQ=255)
                                                                                                                                      ttCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACaa             <  3:440794/149‑1 (MQ=255)
                                                                                                                                        cccGAAGCGGTTCATCAGCTTGTGAATACGCCTCCCCTGTGCGAAGGATAATTACTGCCCGTTCAGCGTGTGCTACAAGCGAAATTTGGATCGTGACTAGCTTCGCACATGCCACTAATAATGTCGCCGTTTTCTCATCTTCAACACTc           <  2:249710/147‑1 (MQ=255)
                                                                                                                                         aCGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCg          >  4:87825/1‑149 (MQ=255)
                                                                                                                                           gAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                               <  1:442993/94‑1 (MQ=255)
                                                                                                                                           gAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                               >  2:442993/1‑94 (MQ=255)
                                                                                                                                           gAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg        >  4:49343/1‑149 (MQ=255)
                                                                                                                                           gAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg        <  3:76067/149‑1 (MQ=255)
                                                                                                                                           gAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg        >  3:300371/1‑149 (MQ=255)
                                                                                                                                             aGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGt      <  3:240089/149‑1 (MQ=255)
                                                                                                                                                 gTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAc  >  3:521084/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCAACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAC  >  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: