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→GGGglpEthiosulfate: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

ACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCT  >  NC_000913/3561483‑3561952
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
aCAGTATAAAGCGTTCCGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTcc                                                                                                                                                                                                                                               <  2:713740/238‑1 (MQ=255)
    tataAAGCGTTACGCGCCGTACGCCGCCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCAtt                                                                                                                                                                                                                                           >  2:173781/1‑238 (MQ=255)
      taaacggttccccgcCGTCCCCCCCCTCTGGGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTACAGACCACATCGTAGCCCTGTTGCAGCAGTTACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  1:10814/227‑1 (MQ=255)
      taAAGGGTTACGCGCCGTACGCCACCCCTGGTGGAATCTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTAGAGCCCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCTCGCATAAAAGCGCCGAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  2:107416/239‑1 (MQ=255)
      taAAGCGTTCCGCGCCGTACGCCCCCTCGGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  2:574374/239‑1 (MQ=255)
      taAAGCGTTACGGGCCGTCCGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCg                                                                                                                                                                                                                                        <  2:670359/239‑1 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAACGCCGCCGTCAATGCTATAGACCACATCGTCGCCCTGTTTCAGCAGCTACTGCCCCGCGCCTTTGCTTCTTTTTCCGTGTTACCTAATCACCTTCACCGGCCGGTCAAAGTCGTTATCACGCATAAAGCCGCCCACCGGTTCGTCGTTTAAATGGAAAGCCTCACCCCCATGTCACATTGCGCAACTCtg                                                                                                                                                                                                                              >  1:71780/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGTTAACACATCACCATCACCGGAGTGTCAAAGGCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCAAGTCCCATTGCGCAACGCtg                                                                                                                                                                                                                              >  2:147569/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGTAAAGCCTGCACCGCATGTCCCATTGCGAAACTCtg                                                                                                                                                                                                                              >  2:890462/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCtg                                                                                                                                                                                                                              >  2:573547/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCCCATGTCCCATTCCGAAActctt                                                                                                                                                                                                                               >  2:240988/1‑238 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGAAAAGCCTGCACCGCATGTCCCATTGCGAAACTCtg                                                                                                                                                                                                                              >  2:227589/1‑239 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTATTCACGCATAAAAACGCCAAGCGTGTCGTTGGTTAACTGGAAAACCGGCACCGCATGTCCCATAGAGAAAatcgg                                                                                                                                                                                                                                >  1:215870/1‑237 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACACCACCATCACCGGAGTGTCAAAGACGTTGTCGCGCCTCAAACCACCCCGCGTGTCGGTGGTTAACAGGAAAGCCTTCACCACAGCCCCCATAGCGaaaagcgt                                                                                                                                                                                                                                >  1:766999/1‑234 (MQ=255)
              tACGCGCCGTACGCCACCTCTGCGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAACCCTGCACCGCATGCCCCATTGCGAAAcccttt                                                                                                                                                                                                                              >  2:389643/1‑237 (MQ=255)
                cgcgCCGTACGCCCCCCGTTCGGGAAACTGAGGTGGCCAGGCTTCAAAGCCGCCGCCAATGCTAGAGACCACATCTTAGCCCGGTTGCAGCAGATACGGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACGCATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCtgtg                                                                                                                                                                                                                           <  1:627004/239‑1 (MQ=255)
                     cGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGCAAGCCTGCACCGCATGTCCCATTGCGCAACTCTGTGGATc                                                                                                                                                                                                                       >  2:383632/1‑238 (MQ=255)
                           cccccTCCGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGGCCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAata                                                                                                                                                                                                               <  2:609727/236‑1 (MQ=255)
                               ctctGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGCCCCATTGCGAAACCCTGTGGTTCGCGAAAATca                                                                                                                                                                                                            >  2:529906/1‑237 (MQ=255)
                                 ctGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGAcc                                                                                                                                                                                                         >  2:462183/1‑239 (MQ=255)
                                                   gCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCCCGAATATCGACCAGCACCGCCTCTTTTTc                                                                                                                                                                                       >  2:386495/1‑238 (MQ=255)
                                                                            gCTATACACCCCATCTTGCCCTGTTTACAGAAGATCCTGCGCCGCGCTTTGGCTTCTATTGCCGTTAAAACACACCCCCATCCCCGGGGTGTAAAAGGCGTTACCCCAAAAAAAAGCGCCCGGCGTGTCGTGGGTTAAATGAAAAGCCTGCCCCCCATTCCCCATTGAGAATCTCTGTGGAGCGCGAATACCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGc                                                                                                                                                               <  1:503430/237‑1 (MQ=255)
                                                                               ataGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGt                                                                                                                                                          >  2:556968/1‑239 (MQ=255)
                                                                                       caTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGt                                                                                                                                                                         >  2:655225/1‑216 (MQ=255)
                                                                                        aTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAata                                                                                                                                                                                                               >  2:620281/1‑178 (MQ=255)
                                                                                        aTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAata                                                                                                                                                                                                               <  1:620281/178‑1 (MQ=255)
                                                                                        aTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCACCTTCTGGGGCGCGCCGGCaaacgtaaaaaaac                                                                                                                                                 >  1:618141/1‑227 (MQ=255)
                                                                                        aTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATa                                                                                                                                                     >  2:503430/1‑235 (MQ=255)
                                                                                        aTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAACCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCCCTTTTTCCTGCCACTTCTGTTGCGCGCCGGCACCGTTAATAcct                                                                                                                                                  >  1:541746/1‑236 (MQ=255)
                                                                                               cccTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGGGCGCGCCGGCAACGTTAATACATTCGAACTg                                                                                                                                          >  2:244955/1‑239 (MQ=255)
                                                                                                     ttGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGCCCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGACCCa                                                                                                                                     >  2:846846/1‑239 (MQ=255)
                                                                                                     ttGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATc                                                                                                                                       >  2:707305/1‑237 (MQ=255)
                                                                                                                        gcgccTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGCGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGTTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGAACGGCGAATATCGCCAAGACCCGCCCCTTTTCCCTGCAACCTCGGGTGCGCGCTGGCCCACTTAATAACTTCGAACTGCGCCACTCGCATTCCTTTCTTTac                                                                                                                  >  1:572910/1‑239 (MQ=255)
                                                                                                                               tggtgttatttccttaaaaaCACACCACCGTCACCGGAGTGTAAAAGTCTGTATCCCGAAAAAAAGCGCCAAGGGTGGCGTGGGTTAAATGGAAAGCCTGCACCGCAGGCCCCATTGGGAACTGTTGGGGATGGGGAATACCGACCAGCCCCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAg                                                                                                           <  1:173781/221‑1 (MQ=255)
                                                                                                                                      aTTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGGGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCaa                                                                                                    >  2:809872/1‑239 (MQ=255)
                                                                                                                                            gTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTCAAAACAAGTGGGCAAATTTa                                                                                               >  2:198542/1‑238 (MQ=255)
                                                                                                                                               aTAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATTTACCAGCACCGCCTCGTTTTCTTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCTTCTCTCTTTTTTTACAAAAAAATGGGCAAATTTAccac                                                                                           >  1:163023/1‑237 (MQ=255)
                                                                                                                                               aTAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACAAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGc                                                                                           >  2:164156/1‑239 (MQ=255)
                                                                                                                                                aaacacctccaccgccccGGTGTGTGAAATGCTTTACCACACAAAAAAGCCCCCAGGGGTGGGTGTGTAAAATGGAAACCCCTCCCCGCTTGCCCCTTTGCGAACTGTGGGGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGca                                                                                          <  1:707305/223‑1 (MQ=255)
                                                                                                                                                        caccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTAc                                                                                  >  2:772150/1‑239 (MQ=255)
                                                                                                                                                        caccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCAACCAGCACCGCCTCTTTTTCCTGCAACTTCTTGTGCGCGCCGGCAACGTTAATACATTCGAACTGATCCATGAGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTAc                                                                                  >  2:419166/1‑239 (MQ=255)
                                                                                                                                                                     gtgtCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGGGCGCGTCGGCAACGTTAATACATTCCAACTGATCCCTGCGTCTCTCTTTTTTTACAAACAAGAGGGCAAATTTACCCCAAAGTTTACGTCGacagggaag                                                                     >  1:731222/1‑231 (MQ=255)
                                                                                                                                                                        tCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCCCTCTTTCTTTCCAACAAAGGGGGCAAATTTACCGCCAAATTTAACTGGAACCGGCaaaaa                                                                  >  1:53381/1‑235 (MQ=255)
                                                                                                                                                                          aaaGTCGTTACCCGGCAAAAAGCCGCCCAGGGTGTCGTTGGTTAAGGGGAAACCCCGCCCCGCATGTCCCATTGGAAACCTTTTGGGATCGGAAATTTCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATaaa                                                                <  2:618141/239‑1 (MQ=255)
                                                                                                                                                                               ttttaccccgaaTAAAAGCCCCAGGGTGGTGTTTGTTAAAAGGAAAACCCTCCCCCGATGTCCCCATGCAAAACCTGTGGGGGTGCGAAAAATCCACCACCCCCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  1:556968/228‑1 (MQ=255)
                                                                                                                                                                               ctttctcccgtaaaaaaGCGCCCGGGGTGTGTTTGGTAAAATGGAAAGCCTCGCCCGTATTCCCCATGGCAAACTCTTGGGGATGGGAAATATCGCCCGGCCCCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  1:207631/226‑1 (MQ=255)
                                                                                                                                                                               ctttatccggaataaagccccccAGGATGCCGTTGGTAAAATGGAACCCCTGCACCGCATGCCCAATTGCGAAACTCTGTGGATCGCAAATTTCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATGAACGCCa                                                           <  2:124111/219‑1 (MQ=255)
                                                                                                                                                                               ctttaccgcGCAAAAAGCCGCCCAGCGGTTCGTGGGTAAATTGGAAAGCCGCCCCCGAAGGCCCCATTGCGAATCTTTGTGGACCGCGAATATCGCCCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  2:541746/231‑1 (MQ=255)
                                                                                                                                                                               ctttATCACGCAAAAACGCGCCCAGCGTGTCTTGTGTAAAGGGGAAAGCCTGCACCGCATTCCCCTTTCCAAACCTCTGTGGTTCGCAATTAGCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  2:454114/237‑1 (MQ=255)
                                                                                                                                                                               cGTTACCCCCCAAAAAAGGCCCCGGCGTTTCTGTGGTTAAAGGGAAAGCCCGCACCGCATGCCCCATTGCAAACCTCTGGGGATCGGAAATACCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCa                                                           <  2:47660/239‑1 (MQ=255)
                                                                                                                                                                                    tcaccaaaaaaggccgcaggCGTGTGGTTGGTTAAAGAGAAACCCGGCAACGCATGCCCAATCGGAAAACTGTGTGAGTAGGGAATTCCGGACGCCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATg                                                      <  2:766999/220‑1 (MQ=255)
                                                                                                                                                                                       cGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTta                                                   >  2:201637/1‑239 (MQ=255)
                                                                                                                                                                                                   cccAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATTTtata                                                 >  2:566986/1‑229 (MQ=255)
                                                                                                                                                                                                   cccAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGGTTACGTCGAAGCGGCAGATAACCGCCATAATGTTATACATATCActc                                       >  2:50655/1‑239 (MQ=255)
                                                                                                                                                                                                            tCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtt                              >  2:871629/1‑239 (MQ=255)
                                                                                                                                                                                                                gggTAAAAGGGAACCCCTGTCCCGACTGTCCCATTGCAAAACTCGGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtttttt                          <  2:818267/238‑1 (MQ=255)
                                                                                                                                                                                                                 ttgtaaaaGGGAAACCCTCCCCCCAATCTCCCTTTGCGAAACTTGTGGGAGCCAGAATCGCCGCAAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtttttt                          <  1:846846/234‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACCTATCACCCTAAAATGTTTTTTCAATg                     >  2:766724/1‑239 (MQ=255)
                                                                                                                                                                                                                             aaGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCCCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATAAATATACCCCTAAAATGTTTTTTAAATTTTACATa              >  1:250391/1‑239 (MQ=255)
                                                                                                                                                                                                                                       cccATGTCCCATTGGAAAACCTTGTGGTTCCGGAATTGCGCCAGAGACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGAtt    <  2:665387/237‑1 (MQ=255)
                                                                                                                                                                                                                                         cATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCCCTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCt  >  2:804697/1‑239 (MQ=255)
                                                                                                                                                                                                                                         cATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCCCTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCt  >  2:893010/1‑239 (MQ=255)
                                                                                                                                                                                                                                         cATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCt  >  2:447514/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
ACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCT  >  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: