Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,044 (T)9→8 intergenic (+41/‑172) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,0360T.100.0% 270.0 / NA 66intergenic (+33/‑180)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (31/35);  total (31/35)

GTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAG  >  NC_000913/1369904‑1370158
                                                                                                                                    |                                                                                                                          
gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                      <  2:757923/139‑1 (MQ=255)
gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                      >  2:2054074/1‑139 (MQ=255)
gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCGGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                      <  1:24132/139‑1 (MQ=255)
    gggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                                                                                  >  1:748780/1‑136 (MQ=255)
     ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctgt                                                                                                                 <  2:901252/139‑5 (MQ=255)
        aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAGAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                             <  1:1418953/139‑1 (MQ=255)
          gcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCacac                                                                                                           >  2:1097129/1‑139 (MQ=255)
                    tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 <  1:2054074/139‑1 (MQ=255)
                          gCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                           <  1:1097129/139‑1 (MQ=255)
                          gCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                           <  2:283051/139‑1 (MQ=255)
                               aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcaca                                                                                                            <  2:1216839/117‑1 (MQ=255)
                               aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 <  2:372484/128‑1 (MQ=255)
                               aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 >  1:372484/1‑128 (MQ=255)
                               aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                                                      >  1:907458/1‑139 (MQ=255)
                               aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCAAGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcaca                                                                                                            >  1:1216839/1‑117 (MQ=255)
                                tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                                                     >  1:827519/1‑139 (MQ=255)
                                        aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                             <  1:1634910/139‑1 (MQ=255)
                                         cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                       >  1:765284/1‑128 (MQ=255)
                                         cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                       <  2:765284/128‑1 (MQ=255)
                                         cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                                            >  2:1275325/1‑139 (MQ=255)
                                            aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                      >  1:1843050/1‑95 (MQ=255)
                                            aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                      <  2:1843050/95‑1 (MQ=255)
                                            aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                             >  2:625424/1‑103 (MQ=255)
                                            aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                             <  1:625424/103‑1 (MQ=255)
                                            aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTca                                                                                                              <  2:423112/102‑1 (MQ=255)
                                            aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTca                                                                                                              >  1:423112/1‑102 (MQ=255)
                                             cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 <  1:1264373/114‑1 (MQ=255)
                                             cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 >  2:1264373/1‑114 (MQ=255)
                                             cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                                                        >  2:479027/1‑139 (MQ=255)
                                                  agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGt                                                                                                                            >  2:1489757/1‑83 (MQ=255)
                                                  agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGt                                                                                                                            <  1:1489757/83‑1 (MQ=255)
                                                         cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATccc                                                                                             >  2:1647865/1‑106 (MQ=255)
                                                         cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATccc                                                                                             <  1:1647865/106‑1 (MQ=255)
                                                            tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                                                         >  1:1427033/1‑139 (MQ=255)
                                                               ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                                                      <  1:1275325/139‑1 (MQ=255)
                                                                 tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                                                    <  2:827519/139‑1 (MQ=255)
                                                                   aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                  <  2:276276/123‑1 (MQ=255)
                                                                   aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                  >  1:276276/1‑123 (MQ=255)
                                                                       aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaa                                              <  2:748780/139‑1 (MQ=255)
                                                                        cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                           >  1:238405/1‑93 (MQ=255)
                                                                        cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                           <  2:238405/93‑1 (MQ=255)
                                                                          tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                             >  1:606902/1‑105 (MQ=255)
                                                                          tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                             <  2:606902/105‑1 (MQ=255)
                                                                             cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACAcc                                                                                                          >  1:1798233/1‑73 (MQ=255)
                                                                             cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACAcc                                                                                                          <  2:1798233/73‑1 (MQ=255)
                                                                                 gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                    >  2:249351/1‑139 (MQ=255)
                                                                                 gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                    <  2:907458/139‑1 (MQ=255)
                                                                                  ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc                                   >  2:631601/1‑139 (MQ=255)
                                                                                            aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                             <  1:241938/55‑1 (MQ=255)
                                                                                            aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                             >  2:241938/1‑55 (MQ=255)
                                                                                               gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaaa                      <  1:249351/139‑1 (MQ=255)
                                                                                                  aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                                                     >  2:1113449/1‑42 (MQ=255)
                                                                                                  aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                                                     <  1:1113449/42‑1 (MQ=255)
                                                                                                  aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTa                                                                                                   >  2:1002753/1‑59 (MQ=255)
                                                                                                  aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTa                                                                                                   <  1:1002753/59‑1 (MQ=255)
                                                                                                   aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACAcc                                                                                                          >  2:1301251/1‑51 (MQ=255)
                                                                                                   aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACAcc                                                                                                          <  1:1301251/51‑1 (MQ=255)
                                                                                                      gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt               <  2:1427033/139‑1 (MQ=255)
                                                                                                        aTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGc             <  1:631601/139‑1 (MQ=255)
                                                                                                         ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                             <  2:1132676/42‑1 (MQ=255)
                                                                                                         ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                             >  1:1132676/1‑42 (MQ=255)
                                                                                                         ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 <  1:2935/54‑1 (MQ=255)
                                                                                                         ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 >  2:2935/1‑54 (MQ=255)
                                                                                                         ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                >  1:241951/1‑119 (MQ=255)
                                                                                                         ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                <  2:241951/119‑1 (MQ=255)
                                                                                                         ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCt            >  2:1819383/1‑139 (MQ=255)
                                                                                                               aTGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg      <  1:479027/139‑1 (MQ=255)
                                                                                                                 gAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc    <  2:1855867/139‑1 (MQ=255)
                                                                                                                   gcggcgCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAg  >  2:1572572/1‑139 (MQ=255)
                                                                                                                    cggcgCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 <  1:717035/43‑1 (MQ=255)
                                                                                                                    cggcgCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                                 >  2:717035/1‑43 (MQ=255)
                                                                                                                        gcTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                  >  1:1599523/1‑70 (MQ=255)
                                                                                                                        gcTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                  <  2:1599523/70‑1 (MQ=255)
                                                                                                                           tATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc    <  1:348107/129‑1 (MQ=255)
                                                                                                                           tATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc    >  2:348107/1‑129 (MQ=255)
                                                                                                                              gcgcCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                <  1:1640593/98‑1 (MQ=255)
                                                                                                                              gcgcCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                >  2:1640593/1‑98 (MQ=255)
                                                                                                                                    |                                                                                                                          
GTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAG  >  NC_000913/1369904‑1370158

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