Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,491,268 T→A 5.1% D228V (GAT→GTT)  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,2680TA 5.1% 271.3 / 4.9 78D228V (GAT→GTT) oxcoxalyl CoA decarboxylase, ThDP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (28/46);  new base A (2/2);  total (30/48)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

ACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTC  >  NC_000913/2491121‑2491416
                                                                                                                                                   |                                                                                                                                                     
gcaaGCATGAAAACGTCATCATTTGTCAGGGCAAAAGAACGCGCAGCTGCCGCAGAAAGTGGATGGGTATCTTCAAGGCTCCTACCTGACAGAGAACTTTTCCGGAATGGAATCTGATCCCTTTCAATAAATTCACGGACAGGTtcatc                                                                                                                                                      <  2:622996/148‑1 (MQ=255)
acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCTCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                      >  2:592566/1‑149 (MQ=255)
acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                      >  3:652421/1‑149 (MQ=255)
acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                      <  4:652421/149‑1 (MQ=255)
acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                      <  1:286521/149‑1 (MQ=255)
acaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                      >  4:440794/1‑149 (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                                                                                                                      >  2:183520/1‑144 (MQ=255)
                                    gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATATGCCGCGCCTTTGCCAAgg                                                                                                                      <  1:183520/144‑1 (MQ=255)
                                    gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                 >  1:536842/1‑149 (MQ=255)
                                    gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                      >  1:653908/1‑124 (MQ=38)
                                    gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                      >  1:371487/1‑124 (MQ=38)
                                    gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                      <  2:371487/128‑5 (MQ=38)
                                    gaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACAAGCTGATGAACAGCTTCGTGAAtttat                                                                                                                                      <  2:653908/128‑5 (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:536842/149‑1 (MQ=255)
                                                      gAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCttt                                                                                               <  2:272730/149‑1 (MQ=255)
                                                       aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                              <  1:69956/149‑1 (MQ=255)
                                                       aaaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtttt                                                                                              <  4:559678/149‑1 (MQ=255)
                                                        aaGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTg                                                                                             >  3:323069/1‑149 (MQ=255)
                                                             gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                        <  3:397136/149‑1 (MQ=255)
                                                             gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                        >  1:24366/1‑149 (MQ=255)
                                                             gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaa                                                                                        <  1:487732/149‑1 (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:616533/149‑1 (MQ=255)
                                                                gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                     >  1:81180/1‑149 (MQ=255)
                                                                  gTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                                   >  2:401282/1‑149 (MQ=255)
                                                                       ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                             <  1:262996/118‑1 (MQ=255)
                                                                       ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGAtaattaat                                                                                                             >  2:262996/1‑118 (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                                                           <  2:516374/149‑1 (MQ=255)
                                                                                          aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCgg                                                           >  1:164550/1‑149 (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:551711/149‑1 (MQ=255)
                                                                                                                          ttCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCaa                           <  1:401282/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                                                                    >  2:442993/1‑94 (MQ=255)
                                                                                                                                        gAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgac                                                                    <  1:442993/94‑1 (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)
                                                                                                                                                   tcaGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTc  >  1:566834/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTTGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTC  >  NC_000913/2491121‑2491416

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: