Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,196,495 T→C intergenic (+122/+372) icd → / ← ymfD isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,4950TC100.0% 245.1 / NA 66intergenic (+122/+372)icd/ymfDisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (44/22);  total (44/22)

GTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAAGTTAGTTATGTATTCACTGTTAT  >  NC_000913/1196348‑1196641
                                                                                                                                                   |                                                                                                                                                  
gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAgg                                                                                                                                                 >  2:657938/1‑151 (MQ=255)
      cgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGcc                                                                                                                                           <  2:249470/151‑1 (MQ=255)
      cgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGcc                                                                                                                                           >  2:144852/1‑151 (MQ=255)
        atcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCgg                                                                                                                                         >  2:77494/1‑151 (MQ=255)
         tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGAGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCAGGt                                                                                                                                        >  1:37124/1‑151 (MQ=255)
             cGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAgg                                                                                                                                    >  2:96293/1‑151 (MQ=255)
               aaaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAggtg                                                                                                                                  >  1:137214/1‑151 (MQ=255)
               aaaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAggtg                                                                                                                                  <  1:587664/151‑1 (MQ=255)
                  aCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAggtggt                                                                                                                                >  1:118564/1‑150 (MQ=255)
                        aaTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTct                                                                                                                          >  1:222886/1‑150 (MQ=255)
                        aaTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTct                                                                                                                          <  2:19854/150‑1 (MQ=255)
                        aaTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTct                                                                                                                          >  1:57915/1‑150 (MQ=255)
                                      ttAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTccc                                                                                                            >  2:635500/1‑150 (MQ=255)
                                      ttAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTccc                                                                                                            >  1:662764/1‑150 (MQ=255)
                                       tAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCcac                                                                                                          >  1:85758/1‑151 (MQ=255)
                                        aaaTTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCcac                                                                                                          >  2:235679/1‑150 (MQ=255)
                                         aaTTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACAt                                                                                                        >  2:322878/1‑151 (MQ=255)
                                          aTTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAGAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACAt                                                                                                        >  2:397475/1‑150 (MQ=255)
                                          aTTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATa                                                                                                       <  1:136269/151‑1 (MQ=255)
                                            ttattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACAt                                                                                                        >  1:482342/1‑148 (MQ=255)
                                            ttattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACAt                                                                                                        >  1:387669/1‑148 (MQ=255)
                                            ttattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACAt                                                                                                        >  1:64931/1‑148 (MQ=255)
                                             tattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACAt                                                                                                        >  1:457916/1‑147 (MQ=255)
                                                taACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATAtt                                                                                                 >  2:271864/1‑151 (MQ=255)
                                                 aaCGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTc                                                                                                <  2:593810/151‑1 (MQ=255)
                                                 aaCGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTc                                                                                                >  2:140300/1‑151 (MQ=255)
                                                  aCGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTcc                                                                                               >  2:610687/1‑151 (MQ=255)
                                                    gggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCtt                                                                                             <  2:264149/151‑1 (MQ=255)
                                                        gCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCt                                                                                         <  2:668353/151‑1 (MQ=255)
                                                         cGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCtt                                                                                        >  1:290663/1‑151 (MQ=255)
                                                           tAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCtttt                                                                                      >  2:433842/1‑151 (MQ=255)
                                                            aaCGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCttttt                                                                                     >  1:661484/1‑151 (MQ=255)
                                                             aCGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCttttt                                                                                     >  2:164778/1‑150 (MQ=255)
                                                               gCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTttatt                                                                                  >  2:667931/1‑151 (MQ=255)
                                                                    cGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGTGGCAGGCTGGCCGGTCAGGTTGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGc                                                                             <  1:65117/151‑1 (MQ=255)
                                                                       tgttTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTg                                                                           <  1:598870/150‑1 (MQ=255)
                                                                        gttTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGc                                                                         >  2:321031/1‑151 (MQ=255)
                                                                            tttGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGa                                                                      <  2:306046/150‑1 (MQ=255)
                                                                                          aaCGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAatat                                                         <  2:293650/149‑1 (MQ=255)
                                                                                           aCGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTa                                                      <  2:616943/151‑1 (MQ=255)
                                                                                            cGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTaa                                                     <  2:614905/151‑1 (MQ=255)
                                                                                            cGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTa                                                      >  1:69761/1‑150 (MQ=255)
                                                                                            cGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTa                                                      >  1:27117/1‑150 (MQ=255)
                                                                                              gTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAgg                                                   >  2:209224/1‑151 (MQ=255)
                                                                                                tATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTc                                                 >  2:188645/1‑151 (MQ=255)
                                                                                                  tCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCtt                                               >  1:129944/1‑151 (MQ=255)
                                                                                                   cAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTc                                              >  1:140112/1‑151 (MQ=255)
                                                                                                    aaaaTTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTc                                              >  1:347513/1‑150 (MQ=255)
                                                                                                    aaaaTTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCa                                             <  1:448393/151‑1 (MQ=255)
                                                                                                       aTTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCAtt                                           >  1:86525/1‑150 (MQ=255)
                                                                                                       aTTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTg                                          >  2:30626/1‑151 (MQ=255)
                                                                                                        ttttATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTg                                          >  1:233874/1‑150 (MQ=255)
                                                                                                          ttATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGAt                                        <  1:66113/150‑1 (MQ=255)
                                                                                                             tCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTa                                      >  1:22862/1‑149 (MQ=255)
                                                                                                                  aaaaGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATc                               <  2:358600/151‑1 (MQ=255)
                                                                                                                   aaaGTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATcc                              >  1:590151/1‑151 (MQ=255)
                                                                                                                      gTTATCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCa                           <  1:49265/151‑1 (MQ=255)
                                                                                                                          tCAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATACCTTAGCTTTTTATTATTGCTGGCGGACGCACGTTAATATTGAAGGTCTTCATTGATTAAGACATCCCCaaaa                        <  1:196319/150‑2 (MQ=255)
                                                                                                                           cAAAACCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAagtt                      <  1:595631/151‑1 (MQ=255)
                                                                                                                             aaaCCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAagtt                      >  1:393328/1‑149 (MQ=255)
                                                                                                                               aCCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAagttagtt                  >  2:673816/1‑151 (MQ=255)
                                                                                                                               aCCCCTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGCTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGGCTTCATTGATTAAGACATCCCCAAagttagtt                  >  1:54385/1‑151 (MQ=255)
                                                                                                                                  ccTCGGTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGAGTAAGACATCCCCAAAGTTAGTTATg               <  2:350085/151‑1 (MQ=255)
                                                                                                                                      ggTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAAGTTAGTTATGTAtt           <  1:140296/151‑1 (MQ=255)
                                                                                                                                      ggTAGTTTTGGGGCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAAGTTAGTTATGTAtt           <  1:434005/151‑1 (MQ=255)
                                                                                                                                               ggggCAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAAGTTAGTTATGTATTCACTGttat  >  2:627049/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTTTATTATTGCTGGCGGACGCTCGTTAATATTTAAGGTCTTCATTGATTAAGACATCCCCAAAGTTAGTTATGTATTCACTGTTAT  >  NC_000913/1196348‑1196641

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