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,2170CA97.2% 192.0 / NA 73W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base A (22/47);  minor base T (1/1);  total (24/49)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913809,2180AG97.3% 249.1 / ‑2.0 73W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (23/48);  total (24/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTC  >  NC_000913/808975‑809458
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                 
gTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGGGACATGGCATCAATCTGCGACTTCATCGCCTCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCCCCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACCAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAATCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATgg                                                                                                                                                                                                                                           <  2:199838/251‑1 (MQ=255)
  ttGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCTCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGc                                                                                                                                                                                                                                          <  1:374417/250‑1 (MQ=255)
              tCAATGGCTGGCTCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGtt                                                                                                                                                                                                                                       <  2:178243/241‑1 (MQ=255)
              tCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGtt                                                                                                                                                                                                                                       >  1:178243/1‑241 (MQ=255)
               cAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTATATGGGTGCAGGATATGGCGTTGGTCAAAGGc                                                                                                                                                                                                                             <  1:560351/250‑1 (MQ=255)
                aaTGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCa                                                                                                                                                                                                                            >  1:504072/1‑250 (MQ=255)
                     ctggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGACATGGATGTGTATGGGTGCAGCATATGGCGTTGGGCAAAGGcc                                                                                                                                                                                                                            >  2:603229/1‑244 (MQ=255)
                         cgcATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCg                                                                                                                                                                                                                  >  1:667506/1‑251 (MQ=255)
                          gcATGGGTGATTCCGCCAAACATACCATGCGACATGGCAACAAGTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATCCCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTGCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGt                                                                                                                                                                                                                 <  2:309002/251‑1 (MQ=255)
                          gcATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGt                                                                                                                                                                                                                 <  1:635887/251‑1 (MQ=255)
                                  gATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGGCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATa                                                                                                                                                                                                          <  1:635881/250‑1 (MQ=255)
                                         ccAAACATGCCATGCGGCATGGCATCAATTTGCGACATCATCGCCGCACTAAGCTGCGGGTGATTGTAGCCGGGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAACCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACtt                                                                                                                                                                                                  <  2:685588/251‑1 (MQ=255)
                                          cAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTg                                                                                                                                                                                                 <  1:531857/251‑1 (MQ=255)
                                           aaaCATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGt                                                                                                                                                                                                <  1:786725/251‑1 (MQ=255)
                                            aaCATCACATGCGCCATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTa                                                                                                                                                                                               <  2:809587/251‑1 (MQ=255)
                                               attacaggcgcaaTGGCATCACTTTGCGACTTCATCGCCGCATGACTCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACATGCGGCGCTCACCACCGGATAAACCGGCAGAGGGTAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGCCAAATGCATGATCGTCCGTTGTCATAATCGACTTGTAAAc                                                                                                                                                                                            <  2:186770/239‑1 (MQ=255)
                                               aTCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAAc                                                                                                                                                                                            <  1:603229/251‑1 (MQ=255)
                                                tCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAAcc                                                                                                                                                                                           >  2:280148/1‑251 (MQ=255)
                                                ccacaTGCTACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGACCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAAcc                                                                                                                                                                                           <  2:578551/250‑1 (MQ=255)
                                                      gCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGGCAACCAGGCGTCTGCCGTCAGACAAAAGCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTg                                                                                                                                                                                     <  2:864856/251‑1 (MQ=255)
                                                        gACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGaa                                                                                                                                                                                   <  1:193091/251‑1 (MQ=255)
                                                         aCATGGAATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCACCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGaaa                                                                                                                                                                                  <  1:368932/251‑1 (MQ=255)
                                                          cATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGAGCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGaaaa                                                                                                                                                                                 <  1:233700/251‑1 (MQ=255)
                                                          cATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGGATGGGTGCTGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGaaa                                                                                                                                                                                  >  2:351614/1‑250 (MQ=255)
                                                           aTGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGGCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGaaaa                                                                                                                                                                                 >  2:689143/1‑250 (MQ=255)
                                                             ggCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGGTGTCATAATCGACTTGTAAACCAAATTGAAAAGAt                                                                                                                                                                              <  1:689143/251‑1 (MQ=255)
                                                                  cAATTTGCGCCTTCTCGGCCGCATTAAGCTGCGGGTGATTGTAGCCGCGGAGCGCCGCCCGCCAGGACGACACACCGTCATCCAGGCGTCTGCCGTCAGACAAAAGCAGCTCGCATCCTGCGGCGCTCCCCTCCGGAGAAACCTGCAGAGGGTAGGTCAGGGATGTGTATGTGTGCAGGATATGGCGTTGGTCAAAGGCAAGACAGTCCGTTTTCACAACCGACTTGTAAATCAAATTGAGAAGATTTAgg                                                                                                                                                                         <  2:564139/251‑1 (MQ=255)
                                                                   aaTTGGCGACTCAATCGCCCCATTAAGCGGCGGGTGTTGTTAGCCGGGGATCCCCGCCCACCAGGACTAAACACCGTCACCCAGGCGTCTGCCGTCAGACCAAATCCGCTCGCCACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGTAGGCCATGGATGCGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCTTTGGCTTAACCGACCTGTAAACCAAATCGAAAAGATTTAGGt                                                                                                                                                                        <  2:345178/251‑1 (MQ=255)
                                                                    aTTTGCGACTTCCTCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACACAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGACCGTTGTCCTAATCGACTTGTAAACCAAATTGAAAAGATTTAGGtt                                                                                                                                                                       >  1:268766/1‑251 (MQ=255)
                                                                              tcagcgccgcATGAAGCTGCGGGTGATATTAGCCGTGGTCCGCCGCCCGCCAGGGCGACATACCGTGAGCCAGGCTACTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGTAGTTGGACAAAGGCAAGATCGTTCGTTGTCATAATCGACGTGTAAACCAAATTTAAAAGATTTAGGTTTACAAGTCt                                                                                                                                                              <  2:789079/246‑1 (MQ=255)
                                                                               cATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTa                                                                                                                                                             <  1:633035/250‑1 (MQ=255)
                                                                                  cgccgcATTAAGCTGCGGGTGATTGTAGCCGTGGCTCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGCTCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAATATTTAGGTTTACAAGTCTACAcc                                                                                                                                                         >  2:188845/1‑251 (MQ=255)
                                                                                  cgccgcATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTacac                                                                                                                                                          <  1:677856/250‑1 (MQ=255)
                                                                                      cgccttAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATTCCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGaa                                                                                                                                                      <  2:466877/247‑1 (MQ=255)
                                                                                       gcATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCTACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAACGGCAAGATCGTCCGCTGTACTAATCGACTTGTAACACAAATTGAAAAGATTAAGGTTTACAAGACTACACCTAAt                                                                                                                                                     >  2:272460/1‑251 (MQ=255)
                                                                                        gcttAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGGCGACATTCCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAAtt                                                                                                                                                    <  2:257527/249‑1 (MQ=255)
                                                                                            aaGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTaaca                                                                                                                                                <  1:520368/251‑1 (MQ=255)
                                                                                            aaGCTGCGGGTGATTGGAGCCGGGGATCGCCGCCCACCAGGACGACACACCGTCAACCAGGCGTCTGCCGTCAGACAAAAGCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGGTTACAAGTCTACACCGAATTaac                                                                                                                                                 <  1:348079/250‑1 (MQ=255)
                                                                                                 gCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAacaaca                                                                                                                                             <  1:442782/249‑1 (MQ=255)
                                                                                                    ggTGACTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACaaaaa                                                                                                                                         >  2:665455/1‑250 (MQ=255)
                                                                                                     gTGATAGTCGCCGTGGATCGCAGCCCACCAGGACGGCATACTGTACACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCAACACCGGATAAACCGGCAGAGGGGAGGTCATGGGTGTGTATGGGTGCAGGATATGGCATTGGTCAAAGACAAGATCGACCGCTGTCATAATCGACTTGTAAACCCATTTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACacaaac                                                                                                                                        >  2:433802/1‑249 (MQ=255)
                                                                                                           gTAGCCGTGGATCGCCGCCCACCAGGACGACATACCTTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGtt                                                                                                                                  <  2:758858/250‑1 (MQ=255)
                                                                                                              aGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTg                                                                                                                               <  1:665455/251‑1 (MQ=255)
                                                                                                                    ggATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGaa                                                                                                                          <  1:603251/250‑1 (MQ=255)
                                                                                                                         gccgccCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACCAGTCTACACCGACTTAACAACAAAAAACACGTTTTGGAGAAGcccc                                                                                                                     >  2:599168/1‑250 (MQ=255)
                                                                                                                           cgccCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGATTTGGAGAAGCCCCATg                                                                                                                  >  2:710434/1‑251 (MQ=255)
                                                                                                                             cccACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGc                                                                                                                >  1:192662/1‑251 (MQ=255)
                                                                                                                                  caGGACGCCATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGCGGTCATGGATGTGTATGGGTGCAGGATATGGCGATGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAGGATTTAGGTTTACAAGTCTACACCGAATTAACAAAAAAAAAACCGTTTTGGAGAAGCCACATGGCTcccc                                                                                                           >  2:9388/1‑248 (MQ=255)
                                                                                                                                        gacATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCAc                                                                                                     <  1:254512/251‑1 (MQ=255)
                                                                                                                                           aTACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCt                                                                                                  <  1:710434/251‑1 (MQ=255)
                                                                                                                                                        ggCGTCTGCCGTCAGACAAAAACAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGACGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCa                                                                                     <  1:426588/251‑1 (MQ=255)
                                                                                                                                                                  gTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAAt                                                                                                                                                     >  2:307394/1‑177 (MQ=255)
                                                                                                                                                                  gTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAAt                                                                                                                                                     <  1:307394/177‑1 (MQ=255)
                                                                                                                                                                   gcAGACAAAATCAGCTCGCAACCTGCGGCGCGCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGAGGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTTTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAATAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCCCGCTGGACATTGTCGCAAGTCACAGAat                                                                          <  1:793745/250‑1 (MQ=255)
                                                                                                                                                                    cAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTAGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAat                                                                          <  1:87920/250‑1 (MQ=255)
                                                                                                                                                                    cAGACAAAATCAGCTCGAAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAGTTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACACCAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAatt                                                                         >  2:404360/1‑251 (MQ=255)
                                                                                                                                                                           aaTCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGaa                                                                  <  2:114152/251‑1 (MQ=255)
                                                                                                                                                                              cAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGaaaaa                                                               <  1:295371/251‑1 (MQ=255)
                                                                                                                                                                              cAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGaaaa                                                                >  1:699250/1‑250 (MQ=255)
                                                                                                                                                                                           tctgcgcTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTAGGGGTGCAGGATATGGCGTTTGTCAAAGGCAAGATCGTCCGTTGGCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAAATAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCCCGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGAt                                                   <  2:805789/247‑1 (MQ=255)
                                                                                                                                                                                           tCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGAt                                                   <  1:151942/250‑1 (MQ=255)
                                                                                                                                                                                             ggCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTTCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctg                                                <  1:537718/251‑1 (MQ=255)
                                                                                                                                                                                             cgcggccaccacCGGAATAACCGGCAAGAGGGAGGGTATTGATGGGTATGGGTGCAGGATATGGGGTTGGGCAAAGGCAAGAACGTCCGTTGTCAAAATCCACTTGTAAACCAAATTGAAAAAATTTAGGTTTACAAGGCTACCCCCAATTAACAACAAAAAACCCGGTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctg                                                <  1:587371/245‑1 (MQ=255)
                                                                                                                                                                                                     ccaccGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAatta                                                                        >  1:81295/1‑219 (MQ=255)
                                                                                                                                                                                                     ccaccGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAatta                                                                        <  2:81295/219‑1 (MQ=255)
                                                                                                                                                                                                     ccaccGGATAAACCGGAAGAGGGGAGGTCATGGATGTGTATGGGTGCTGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCAAAATCGACGTGTAAACCAAATGGAAAAGATTTAGGTTTACAAGTCTACACCGAAATAACAACAAATAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCGGGACAGTGTCGCAAGTCACCGAATGATTTGAAAACCTGTTGCTGGATCTGCTGTTTGa                                        <  2:200200/251‑1 (MQ=255)
                                                                                                                                                                                                         cGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGAAGGATATGGCGTTGGTCCAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGCGAAGCCACAAGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAgcg                                    >  1:408648/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGATGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGcagcagtt                            >  2:274550/1‑249 (MQ=255)
                                                                                                                                                                                                                       gagGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCg                      >  2:140702/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggggAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCTACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACAACGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGcc                    <  2:179815/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCg           <  1:643489/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGCAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGCATTATTTGAAAAAACGTTGCTGGATCTGCTGTTTGAAGCGCAACAGGTGCATAGCCAGCATTTCg           >  2:673472/1‑250 (MQ=255)
                                                                                                                                                                                                                                             aTGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTc  <  1:532823/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                 
GTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTC  >  NC_000913/808975‑809458

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