Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,561,718 T→C 16.8% G35G (GGA→GGG)  glpE ← thiosulfate:cyanide sulfurtransferase (rhodanese)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,561,7180TC16.8% 145.5 / 4.5 63G35G (GGA→GGG) glpEthiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base T (38/13);  new base C (3/9);  total (41/22)
Fisher's exact test for biased strand distribution p-value = 2.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCT  >  NC_000913/3561483‑3561952
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
aCAGTATAAAGCGTTCCGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTcc                                                                                                                                                                                                                                               <  2:713740/238‑1 (MQ=255)
    tataAAGCGTTACGCGCCGTACGCCGCCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCAtt                                                                                                                                                                                                                                           >  2:173781/1‑238 (MQ=255)
      taaacggttccccgcCGTCCCCCCCCTCTGGGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTACAGACCACATCGTAGCCCTG‑TTGCAGCAGTTACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  1:10814/227‑1 (MQ=255)
      taAAGGGTTACGCGCCGTACGCCACCCCTGGTGGAATCTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTAGAGCCCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCTCGCATAAAAGCGCCGAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  2:107416/239‑1 (MQ=255)
      taAAGCGTTCCGCGCCGTACGCCCCCTCGGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  2:574374/239‑1 (MQ=255)
      taAAGCGTTACGGGCCGTCCGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  2:670359/239‑1 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAACGCCGCCGTCAATGCTATAGACCACATCGTCGCCCTG‑TTTCAGCAGCTACTGCCCCGCGCCTTTGCTTCTTTTTCCGTGTTACCTAATCACCTTCACCGGCCGGTCAAAGTCGTTATCACGCATAAAGCCGCCCACCGGTTCGTCG‑TTTAAATGGAAAGCCTCACCCCCATGTCACATTGCGCAACTC‑t‑g                                                                                                                                                                                                                              >  1:71780/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGTTAACACATCACCATCACCGGAGTGTCAAAGGCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCAAGTCCCATTGCGCAACGC‑t‑g                                                                                                                                                                                                                              >  2:147569/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGTAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑t‑g                                                                                                                                                                                                                              >  2:890462/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑t‑g                                                                                                                                                                                                                              >  2:573547/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCCCATGTCCCATTCCGAAActc‑tt                                                                                                                                                                                                                               >  2:240988/1‑238 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGAAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑t‑g                                                                                                                                                                                                                              >  2:227589/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTATTCACGCATAAAAACGCCAAGCGTGTCGTTG‑GTTAACTGGAAAACCGGCACCGCATGTCCCATAGAGAAAatcgg                                                                                                                                                                                                                                >  1:215870/1‑237 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACACCACCATCACCGGAGTGTCAAAGACGTTGTCGCGCCTCAAACCACCCCGCGTGTCGGTG‑GTTAACAGGAAAGCCTTCACCACAGCCCCCATAGCGaaaagcgt                                                                                                                                                                                                                                >  1:766999/1‑234 (MQ=255)
              tACGCGCCGTACGCCACCTCTGC‑GGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAACCCTGCACCGCATGCCCCATTGCGAAAccc‑ttt                                                                                                                                                                                                                              >  2:389643/1‑237 (MQ=255)
                cgcgCCGTACGCCCCCCGTTCGGGAAACTGAGGTGGCCAGGCTTCAAAGCCGCCGCCAATGCTAGAGACCACATCTTAGCCCGG‑TTGCAGCAGATACGGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACGCATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑t‑gt‑g                                                                                                                                                                                                                           <  1:627004/239‑1 (MQ=255)
                     cGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGCAAGCCTGCACCGCATGTCCCATTGCGCAACTC‑T‑GT‑GGATc                                                                                                                                                                                                                       >  2:383632/1‑238 (MQ=255)
                           cccccTCCGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGGCCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑Aata                                                                                                                                                                                                               <  2:609727/236‑1 (MQ=255)
                               ctctGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGCCCCATTGCGAAACCC‑T‑GT‑GGTTCGCG‑AAAATca                                                                                                                                                                                                            >  2:529906/1‑237 (MQ=255)
                                 ctGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGAcc                                                                                                                                                                                                         >  2:462183/1‑239 (MQ=255)
                                                   gCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCCCG‑AATATCGACCAG‑CACCGCCTCTTTTTc                                                                                                                                                                                       >  2:386495/1‑238 (MQ=255)
                                                                            gCTATACACCCCATC‑TTGCCCTGTTTACAGAAGATCCTGCGCCGCGCTTTGGCTTCTATTGCCGTTAAAACACACCCCCATCCCCGGGGTGTAAAAGGCGTTACCCCAAAAAAAAGCGCCCGGCGTGTCGTGG‑GTTAAATGAAAAGCCTGCCCCCCATTCCCCATTGAGAATCTC‑T‑GT‑GGAGCGCG‑AATACCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGc                                                                                                                                                               <  1:503430/237‑1 (MQ=255)
                                                                               ataGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGt                                                                                                                                                          >  2:556968/1‑239 (MQ=255)
                                                                                       caTCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGt                                                                                                                                                                         >  2:655225/1‑216 (MQ=255)
                                                                                        aTCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑Aata                                                                                                                                                                                                               >  2:620281/1‑178 (MQ=255)
                                                                                        aTCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑Aata                                                                                                                                                                                                               <  1:620281/178‑1 (MQ=255)
                                                                                        aTCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCACCTTCTGGGGCGCGCCGGCaaacgtaaaaaaac                                                                                                                                                 >  1:618141/1‑227 (MQ=255)
                                                                                        aTCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATa                                                                                                                                                     >  2:503430/1‑235 (MQ=255)
                                                                                        aTCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAACCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCCCTTTTTCCTGCCACTTCTGTTGCGCGCCGGCACCGTTAATAcct                                                                                                                                                  >  1:541746/1‑236 (MQ=255)
                                                                                               cccTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGGGCGCGCCGGCAACGTTAATACATTCGAACTg                                                                                                                                          >  2:244955/1‑239 (MQ=255)
                                                                                                     ttGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGCCCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGACCCa                                                                                                                                     >  2:846846/1‑239 (MQ=255)
                                                                                                     ttGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATc                                                                                                                                       >  2:707305/1‑237 (MQ=255)
                                                                                                                        gcgccTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGCGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑TTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GAACGGCG‑AATATCGCCAAG‑ACCCGCCCCTTTTCCCTGCAACCTCGGGTGCGCGCTGGCCCACTTAATAACTTCGAACTGCGCCACTCGCATTCCTTTCTTTac                                                                                                                  >  1:572910/1‑239 (MQ=255)
                                                                                                                               tggtgttatttccttaaaaaCACACCACCGTCACCGGAGTGTAAAAGTCTGTATCCCGAAAAAAAGCGCCAAGGGTGGCGTGG‑GTTAAATGGAAAGCCTGCACCGCAGGCCCCATTGGGAACTGT‑T‑GG‑GGATGGGG‑AATACCGACCAG‑CCCCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAg                                                                                                           <  1:173781/221‑1 (MQ=255)
                                                                                                                                      aTTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GG‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCaa                                                                                                    >  2:809872/1‑239 (MQ=255)
                                                                                                                                            gTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTCAAAACAAGTGGGCAAATTTa                                                                                               >  2:198542/1‑238 (MQ=255)
                                                                                                                                               aTAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATTTACCAG‑CACCGCCTCGTTTTCTTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCTTCTCTCTTTTTTTACAAAAAAATGGGCAAATTTAccac                                                                                           >  1:163023/1‑237 (MQ=255)
                                                                                                                                               aTAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACAAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGc                                                                                           >  2:164156/1‑239 (MQ=255)
                                                                                                                                                aaacacctccaccgccccGGTGTGTGAAATGCTTTACCACACAAAAAAGCCCCCAG‑GGGTGGGTGTGTAAAATGGAAACCCCTCCCCGCTTGCCCCTTTGCG‑AACTG‑TGGG‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGca                                                                                          <  1:707305/223‑1 (MQ=255)
                                                                                                                                                        caccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTAc                                                                                  >  2:772150/1‑239 (MQ=255)
                                                                                                                                                        caccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCAACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTTGTGCGCGCCGGCAACGTTAATACATTCGAACTGATCCATGAGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTAc                                                                                  >  2:419166/1‑239 (MQ=255)
                                                                                                                                                                     gtgtCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGGGCGCGTCGGCAACGTTAATACATTCCAACTGATCCCTGCGTCTCTCTTTTTTTACAAACAAGAGGGCAAATTTACCCCAAAGTTTACGTCGacagggaag                                                                     >  1:731222/1‑231 (MQ=255)
                                                                                                                                                                        tCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCCCTCTTTCTTTCCAACAAAGGGGGCAAATTTACCGCCAAATTTAACTGGAACCGGCaaaaa                                                                  >  1:53381/1‑235 (MQ=255)
                                                                                                                                                                          aaaGTCGTTACCCGGCAAAAAGCCGCCCAGGGTGTCGTTG‑GTTAAGGGGAAACCCCGCCCCGCATGTCCCATTGGAAACCTT‑T‑TG‑GGATCGGA‑AATTTCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATaaa                                                                <  2:618141/239‑1 (MQ=255)
                                                                                                                                                                               ttttaccccgaaTAAAAGC‑CCCAGGGTGGTGTTT‑GTTAAAAGGAAAACCCTCCCCCGATGTCCCCATGCAAAAC‑C‑T‑GTGGGGGTGCGAAAAATCCACCAC‑CCCCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  1:556968/228‑1 (MQ=255)
                                                                                                                                                                               ctttctcccgtaaaaaaGCGCCCGGGGTGTGTTTG‑GTAAAATGGAAAGCCTCGCCCGTATTCCCCATGGC‑AAACTCTT‑GG‑GGATGGGA‑AATATCGCCCGG‑CCCCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  1:207631/226‑1 (MQ=255)
                                                                                                                                                                               ctttatccggaataaagccccccAGGATGCCGTTG‑GTAAAATGGAACCCCTGCACCGCATGCCCAATTGCGAAACTC‑T‑GT‑GGATCGCA‑AATTTCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATGAACGCCa                                                           <  2:124111/219‑1 (MQ=255)
                                                                                                                                                                               ctttaccgcGCAAAAAGCCGCCCAGCGGTTCGTGG‑GTAAATTGGAAAGCCGCCCCCGAAGGCCCCATTGCGAATCTT‑T‑GT‑GGACCGCG‑AATATCGCCCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  2:541746/231‑1 (MQ=255)
                                                                                                                                                                               ctttATCACGCAAAAACGCGCCCAGCGTGTCTTGT‑GTAAAGGGGAAAGCCTGCACCGCATTCCCCTTTCCAAACCTC‑T‑GT‑GGTTCGCA‑ATTAGCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  2:454114/237‑1 (MQ=255)
                                                                                                                                                                               cGTTACCCCCCAAAAAAGGCCCCGGCGTTTCTGTG‑GTTAAAGGGAAAGCCCGCACCGCATGCCCCATTGCAAACCTC‑T‑GG‑GGATCGGA‑AATACCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  2:47660/239‑1 (MQ=255)
                                                                                                                                                                                    tcaccaaaaaaggccgcaggCGTGTGGTTG‑GTTAAAGAGAAACCCGGCAACGCATGCCCAATCGGAAAACTG‑T‑GT‑GAGTAGGG‑AATTCCGGACGC‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATg                                                      <  2:766999/220‑1 (MQ=255)
                                                                                                                                                                                       cGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTta                                                   >  2:201637/1‑239 (MQ=255)
                                                                                                                                                                                                   cccAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATTTtata                                                 >  2:566986/1‑229 (MQ=255)
                                                                                                                                                                                                   cccAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGGTTACGTCGAAGCGGCAGATAACCGCCATAATGTTATACATATCActc                                       >  2:50655/1‑239 (MQ=255)
                                                                                                                                                                                                            tCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtt                              >  2:871629/1‑239 (MQ=255)
                                                                                                                                                                                                                gg‑gTAAAAGGGAACCCCTGTCCCGACTGTCCCATTGCAAAACTC‑G‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtttttt                          <  2:818267/238‑1 (MQ=255)
                                                                                                                                                                                                                 ttgtaaaaGGGAAACCCTCCCCCCAATCTCCCTTTGCGAAACT‑‑T‑GTGGGAGCCAG‑AATCGCCGCAAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtttttt                          <  1:846846/234‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACCTATCACCCTAAAATGTTTTTTCAATg                     >  2:766724/1‑239 (MQ=255)
                                                                                                                                                                                                                             aaGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCCCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATAAATATACCCCTAAAATGTTTTTTAAATTTTACATa              >  1:250391/1‑239 (MQ=255)
                                                                                                                                                                                                                                       cccATGTCCCATTGGAAAACCT‑T‑GT‑GGTTCCGG‑AATTGCG‑CCAGAGACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGAtt    <  2:665387/237‑1 (MQ=255)
                                                                                                                                                                                                                                         cATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCCCTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCt  >  2:804697/1‑239 (MQ=255)
                                                                                                                                                                                                                                         cATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCCCTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCt  >  2:893010/1‑239 (MQ=255)
                                                                                                                                                                                                                                         cATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCt  >  2:447514/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
ACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTG‑TTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTG‑GTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTC‑T‑GT‑GGATCGCG‑AATATCGACCAG‑CACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCT  >  NC_000913/3561483‑3561952

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