Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 809,217 2 bp→AG coding (40‑41/1290 nt) bioA ← adenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913809,2170CA98.2% 188.4 / ‑4.5 55coding (41/1290 nt)bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (21/33);  minor base . (1/0);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913809,2180AG98.2% 187.5 / NA 55W14R (TGG→CGG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/33);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTC  >  NC_000913/808978‑809458
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                   
tGCGGCACAGCTCAATGGCTGGCGCATGGGTGATCCCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCg                                                                                                                                                                                                                                           <  2:183254/250‑1 (MQ=255)
             aaTGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTGCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCaa                                                                                                                                                                                                                             <  2:366061/251‑1 (MQ=255)
                      cgcATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCg                                                                                                                                                                                                                    >  1:438896/1‑251 (MQ=255)
                       gcATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGt                                                                                                                                                                                                                   >  1:155031/1‑251 (MQ=255)
                           gggTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑c                                                                                                                                                                                                              >  1:177075/1‑251 (MQ=255)
                               gATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCAGTAAGGTGCGGGTGATTGTAGCCGGGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTGACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTTTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATaa                                                                                                                                                                                                          <  2:537020/251‑1 (MQ=255)
                                aTACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAAt                                                                                                                                                                                                         >  1:475093/1‑251 (MQ=255)
                                  ccccgcAAACATCTAATGCGACATTGCATCAATGTGCGACTTCATCGCCGCAATAAGCTGCGGGTTATTGGAGCCGTGGACCGCCGCCCGCCAGGACGACATCCCGTCAACCAGGCGTCTTCCGTCAGACAAAATCAGCTGGCGATGTGCGGCGCTCACCACCGGATAAACCGGCAGAGGGGATGTCATGCATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGTCAAGATCGTCCGTTGT‑CATAATc                                                                                                                                                                                                        <  2:293023/245‑1 (MQ=255)
                                  aCCGCCAAACATCGCAGGCGACATGGCATCAATTTGCGACTTCATCGCCTCATTAATCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCATCCAGGCGTCTGCCGTCAGCCATAATCAGCTCGCATCCTTCGGCGCTCACCACCGGATAAACCGTCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑GTATGGCGTTGGTCAAAGGCACGATCGTCCGTTGT‑CATAATCg                                                                                                                                                                                                       <  2:222005/251‑1 (MQ=255)
                                     gCCAATCATCACATGCGACAGGGCATCAAGTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCGGCCGTCAGACAAAATCAGCGCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACt                                                                                                                                                                                                    <  2:253572/251‑1 (MQ=255)
                                       cAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTg                                                                                                                                                                                                  <  2:107804/251‑1 (MQ=255)
                                         aaCATCACATGCAACATGGCATCAATTTGCGACTTCATCGCAGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACCACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCCCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGCTGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATAGTCCGTTGT‑CATAATCGCCTTGt                                                                                                                                                                                                 >  2:244761/1‑250 (MQ=255)
                                           cATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑aaa                                                                                                                                                                                             <  1:412539/251‑1 (MQ=255)
                                             tCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCCGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAAcc                                                                                                                                                                                           >  1:433505/1‑251 (MQ=255)
                                              cacaTGCGACATGGCAGCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGAGAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGGCTTGT‑AAAcc                                                                                                                                                                                           <  2:155031/250‑1 (MQ=255)
                                                      aCATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGaaa                                                                                                                                                                                  <  1:479446/251‑1 (MQ=255)
                                                            cATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGCTTGAAGCCGTGGATCACCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGACAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGGTGTGTATGGGTGCAGG‑GTATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGAttt                                                                                                                                                                            >  1:514959/1‑251 (MQ=255)
                                                                   ttGCGACTTCACCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGTATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGG‑AAACCAAATTGAAAAGAGTTAGGttt                                                                                                                                                                      <  2:57366/250‑1 (MQ=255)
                                                                      cGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACACAATCAGCTCGCAACCTTCGGCGCTCAACCCCGGCTAAACAGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATCTGGAGTTGGTCAAAGGCAAGATCGTCCGTTGC‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACa                                                                                                                                                                   >  2:206359/1‑250 (MQ=255)
                                                                       gACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAg                                                                                                                                                                 <  2:438896/251‑1 (MQ=255)
                                                                        aCTTCAGCGCCGCATAACGCGGCGGGTGATTGCATCCCTGGCTCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACCAAAGGAGCTCGCAACCTTCGGCGCTCACCAGCGGCTAAAGCGGCAGTGGGGAGGTAATGGATGTGTAGGGGTGC‑AGTATATGGGGCTGACAAAAGGCAAGATCGTCCGTTGA‑CATAATCGACTTGT‑AAACCAAATTGAAAAGAGTTAGGTTTaaaatt                                                                                                                                                                >  2:392698/1‑246 (MQ=255)
                                                                            cagcgccgcATTAAGCGGCTGGTGATGGTAGCCGCGGATCGCCGCCCAACAGGACGTCAGACCGGCAACCAGGGGTCTGCCGTCAGTCAAAATCAGCTCGCAACCTTCGGAGCACACCACCGGATAAACCGGCAGAGGGGAGGTCAGGTATGTGTATGGGTGCAGG‑ATAAGGCGTTGGGAAAATGTAAGATCGTCCGTTGT‑TATAATCGACT‑GTGATACCAAATTGAAAAGATTTATGTTTACAATTCTac                                                                                                                                                            <  2:432554/248‑1 (MQ=255)
                                                                               cgccgcATTAAGCTGCTGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAATGCAAGATCGTCCGGTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTacac                                                                                                                                                          <  1:22666/250‑1 (MQ=255)
                                                                                   gcATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAAt                                                                                                                                                     <  1:342022/251‑1 (MQ=255)
                                                                                                      tGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGt                                                                                                                                   <  2:83503/250‑1 (MQ=255)
                                                                                                         aGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGtttt                                                                                                                                <  1:107439/250‑1 (MQ=255)
                                                                                                            cGTGGATCGCCGCCCACCAGGACGACATACCGGCAACCAGGCGTCTTCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGGATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTgtag                                                                                                                            <  1:372680/251‑4 (MQ=255)
                                                                                                              tGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGgaga                                                                                                                           >  1:344983/1‑250 (MQ=255)
                                                                                                                     ccgccCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCATCCTGCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCAt                                                                                                                   <  1:279815/251‑1 (MQ=255)
                                                                                                                        cccACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGc                                                                                                                >  1:187638/1‑251 (MQ=255)
                                                                                                                           accaGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTACCAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTc                                                                                                              >  2:18629/1‑250 (MQ=255)
                                                                                                                                  cgacATACCGTCAAACAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAAACTTCGGCGCTCACCAACGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGGTGT‑CATAATCGACTTGT‑AAACCACATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAAAAAAAAACACGTTTTGGAGAAGCCACATGGCTCACCGcca                                                                                                       >  2:61989/1‑249 (MQ=255)
                                                                                                                                        cccGTCAACCAGGCGTCTGACGTCAGACAAAATCAGCTCGCACCCTTCGGCTCTCACCACCGGATAAACCGGCATAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCTAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCTCGCTgg                                                                                                <  2:408063/250‑1 (MQ=255)
                                                                                                                                                ccAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTc                                                                                        <  1:105939/251‑1 (MQ=255)
                                                                                                                                                ccAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCGCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTc                                                                                        <  2:475093/251‑1 (MQ=255)
                                                                                                                                                     cGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAg                                                                                   <  1:338816/251‑1 (MQ=255)
                                                                                                                                                                  aCAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCaa                                                                                                                                                                                                                             >  1:368869/1‑102 (MQ=255)
                                                                                                                                                                  aCAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCaa                                                                                                                                                                                                                             <  2:368869/102‑1 (MQ=255)
                                                                                                                                                                       aTCAGCTCGCAACCTTCGGCGCTCACCACCTGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGaaa                                                                 <  1:506362/251‑1 (MQ=255)
                                                                                                                                                                                catccTTCGGCGCTCACCACCGTATAATCCGGCAGAGGGGAGGGCATGGATGTGTATGGGTGCAGG‑ATATGGCTTTGGTCAAAGGCAAGATCGTCCGGTGT‑CATAATCGACTTGG‑AAACCAAATTGAAAAGATTGAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAATCGTTGc                                                        <  2:194092/248‑1 (MQ=255)
                                                                                                                                                                                   ccTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTgg                                                     <  2:85264/251‑1 (MQ=255)
                                                                                                                                                                                        ggCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctg                                                <  1:224134/251‑1 (MQ=255)
                                                                                                                                                                                                  accGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATct                                                 >  2:76770/1‑240 (MQ=255)
                                                                                                                                                                                                  accGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATct                                                 <  1:76770/240‑1 (MQ=255)
                                                                                                                                                                                                           aCCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAACAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAgg                             >  1:372194/1‑251 (MQ=255)
                                                                                                                                                                                                              ggCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGc                          >  1:341965/1‑251 (MQ=255)
                                                                                                                                                                                                                 agagGGGAGGGCATGGATGTGTATGTGTGCAGG‑ATATGGCGGTGGTTAAAGGCAAGATCGTCCGGTGT‑CATAATCGACTTGG‑AAAACAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGTCCACGCTGGACATTGTCGAAAGTCACAGAATTATTTGAAAAACCGTTGCTGGAGCTGCTGTGTGAAGCGCAGTAGGTggggc                       <  2:524009/251‑5 (MQ=255)
                                                                                                                                                                                                                            cATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGAGTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTc            <  2:177075/251‑1 (MQ=255)
                                                                                                                                                                                                                              tGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAg                                                                             >  1:139130/1‑184 (MQ=255)
                                                                                                                                                                                                                              tGGATGTGTATGGGTGCAGG‑ATATGGCGTCGGTCAAAGGCAAGATCGTCCGT‑GTGCATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACACCAAAAAACCCTTTTGGGAGAAGCCCCATGGCTCCCCGCCCACGCTGGACATTGTCGCATGTCACAg                                                                             <  2:139130/184‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGTGTATGGGTGCAGG‑ATACGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGCCTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACAGAATTAACAACAAAAAACAAGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATcc       >  2:445486/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tgtgTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCt      >  1:297002/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTc  >  1:113446/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTc  <  2:341965/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTc  >  1:161927/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                   
TGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGT‑CATAATCGACTTGT‑AAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTC  >  NC_000913/808978‑809458

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