Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,561,667 G→A 5.5% R53C (CGT→TGT)  glpE ← thiosulfate:cyanide sulfurtransferase (rhodanese)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,561,6670GA 5.5% 236.4 / 3.4 73R53C (CGT→TGT) glpEthiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base A (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAA  >  NC_000913/3561520‑3561814
                                                                                                                                                   |                                                                                                                                                   
cGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGc                                                                                                                                                    <  2:370052/149‑1 (MQ=255)
 gggAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCa                                                                                                                                                   >  2:120443/1‑149 (MQ=255)
  ggAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCAt                                                                                                                                                  >  3:205430/1‑149 (MQ=255)
     aaCTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATaaa                                                                                                                                               >  4:34594/1‑149 (MQ=255)
      aCTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATaaaa                                                                                                                                              <  3:438332/149‑1 (MQ=255)
        tGACGTTGCCAGGCTGCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAgc                                                                                                                                            <  4:161950/149‑1 (MQ=255)
            gTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGccc                                                                                                                                        >  1:279210/1‑149 (MQ=255)
              tGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAg                                                                                                                                      <  1:235018/149‑1 (MQ=255)
               gCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGc                                                                                                                                     <  4:205430/149‑1 (MQ=255)
                ccAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCg                                                                                                                                    <  3:34594/149‑1 (MQ=255)
                 cAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCa                                                                                                                                                   <  4:312242/133‑1 (MQ=255)
                 cAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCa                                                                                                                                                   >  3:312242/1‑133 (MQ=255)
                 cAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCgt                                                                                                                                   <  4:174984/149‑1 (MQ=255)
                  aGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCgtg                                                                                                                                  <  1:120443/149‑1 (MQ=255)
                  aGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACACATAAAAGCGCCCAGCgtg                                                                                                                                  >  4:320478/1‑149 (MQ=255)
                  aGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACACATAAAAGCGCCCAGCgtg                                                                                                                                  >  2:43095/1‑149 (MQ=255)
                       tCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGtt                                                                                                                             <  3:189689/149‑1 (MQ=255)
                        cAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTg                                                                                                                            <  2:301434/149‑1 (MQ=255)
                         aaaGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACACATAAAAGCGCCCAGCGTGTCGTTgg                                                                                                                           <  3:320478/149‑1 (MQ=255)
                         aaaGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACACATAAAAGCGCCCAGCGTGTCGTTgg                                                                                                                           <  1:43095/149‑1 (MQ=255)
                           aGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACTCATAAAAGCGCCCAGCGTGTCGTTGGtt                                                                                                                         <  3:407806/149‑1 (MQ=255)
                                  gTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGa                                                                                                                  >  3:364478/1‑149 (MQ=255)
                                       tGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGcc                                                                                                             >  1:243821/1‑149 (MQ=255)
                                              gACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCAGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCg                                                                                                      >  2:44435/1‑149 (MQ=255)
                                               aCCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGc                                                                                                     <  4:294975/149‑1 (MQ=255)
                                                ccACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCa                                                                                                    >  3:357618/1‑149 (MQ=255)
                                                 cacaTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCAt                                                                                                   >  3:275485/1‑149 (MQ=255)
                                                 cacaTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCAt                                                                                                   <  4:120902/149‑1 (MQ=255)
                                                 cacaTCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCAt                                                                                                   >  1:408868/1‑149 (MQ=255)
                                                           cccTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGc                                                                                         <  1:44435/149‑1 (MQ=255)
                                                             cTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGa                                                                                       >  3:53586/1‑149 (MQ=255)
                                                               gTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGaaa                                                                                     <  4:357618/149‑1 (MQ=255)
                                                                   cagcagATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAActct                                                                                 <  3:179500/149‑1 (MQ=255)
                                                                     gcagATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCtgt                                                                               <  2:279210/149‑1 (MQ=255)
                                                                      cagATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCtgtg                                                                              >  3:196010/1‑149 (MQ=255)
                                                                       agATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTgg                                                                             >  1:3575/1‑149 (MQ=255)
                                                                               cgccgcgccTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGaa                                                                     <  4:53586/149‑1 (MQ=255)
                                                                                  cgcgccTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAatat                                                                  >  2:73568/1‑149 (MQ=255)
                                                                                     gccTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGa                                                               >  3:74674/1‑149 (MQ=255)
                                                                                     gccTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGa                                                               <  3:296732/149‑1 (MQ=255)
                                                                                        tttGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCa                                                            >  1:284996/1‑149 (MQ=255)
                                                                                        tttGCTGCTATTGCCGTGATAACACATCACCATCACCGGACTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCa                                                            >  4:276680/1‑149 (MQ=255)
                                                                                             tgctATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCAcc                                                       >  3:142898/1‑149 (MQ=255)
                                                                                             tgctATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCAcc                                                       >  4:276908/1‑149 (MQ=255)
                                                                                                     ccGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCttt                                               >  4:423831/1‑149 (MQ=255)
                                                                                                     ccGTGATAACACACCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCTAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCttt                                               <  2:408868/149‑1 (MQ=255)
                                                                                                      cGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAgc                                                                                                                                            <  2:53416/55‑1 (MQ=255)
                                                                                                      cGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAgc                                                                                                                                            >  1:53416/1‑55 (MQ=255)
                                                                                                         gATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTcc                                           <  4:142898/149‑1 (MQ=255)
                                                                                                          aTAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCt                                          <  2:284996/149‑1 (MQ=255)
                                                                                                            aaCACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGc                                        >  2:453864/1‑149 (MQ=255)
                                                                                                            aaCACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGc                                        >  1:83374/1‑149 (MQ=255)
                                                                                                              cacaTCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCaa                                      >  3:438313/1‑149 (MQ=255)
                                                                                                                 atcaccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATc                                                                          >  3:470442/1‑110 (MQ=255)
                                                                                                                 atcaccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATc                                                                          <  4:470442/110‑1 (MQ=255)
                                                                                                                 atcaccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACtt                                   >  2:293044/1‑149 (MQ=255)
                                                                                                                  tcaccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTc                                  <  4:275485/149‑1 (MQ=255)
                                                                                                                    accatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTg                                >  2:282342/1‑149 (MQ=255)
                                                                                                                    accatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTg                                <  1:293044/149‑1 (MQ=255)
                                                                                                                     ccatcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGGCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTgg                               <  2:243821/149‑1 (MQ=255)
                                                                                                                      cgtcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGaa                                                                     >  4:326532/3‑110 (MQ=255)
                                                                                                                      catcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGaa                                                                     <  3:326532/110‑1 (MQ=255)
                                                                                                                       atcaccGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTg                             >  3:362588/1‑149 (MQ=255)
                                                                                                                        tcaccGGAGTGTCAAAGTCGTTATCGCGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTgc                            <  3:12528/149‑1 (MQ=255)
                                                                                                                          accGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTgcgc                          >  4:423283/1‑149 (MQ=255)
                                                                                                                                     aaaGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGt               >  4:292053/1‑149 (MQ=255)
                                                                                                                                        gTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTaa            <  3:292053/149‑1 (MQ=255)
                                                                                                                                           gTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATAc         >  4:328627/1‑149 (MQ=255)
                                                                                                                                           gTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATAc         >  1:60483/1‑149 (MQ=255)
                                                                                                                                                 aCGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATcg                                                                         >  2:157102/1‑79 (MQ=255)
                                                                                                                                                 aCGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATcg                                                                         <  1:157102/79‑1 (MQ=255)
                                                                                                                                                  cGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGaa  <  4:344492/149‑1 (MQ=255)
                                                                                                                                                  cGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGaa  >  3:344492/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGGAAACTGACGTTGCCAGGCTTCAAAGCCGCCGTCAATGCTATAGACCACATCGTAGCCCTGTTGCAGCAGATACTGCGCCGCGCCTTTGCTGCTATTGCCGTGATAACACATCACCATCACCGGAGTGTCAAAGTCGTTATCACGCATAAAAGCGCCCAGCGTGTCGTTGGTTAAATGGAAAGCCTGCACCGCATGTCCCATTGCGAAACTCTGTGGATCGCGAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAA  >  NC_000913/3561520‑3561814

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: