Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,948,631 A→G L447L (TTG→CTG)  yifB ← putative magnesium chelatase YifB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,948,6310AG100.0% 252.9 / NA 69L447L (TTG→CTG) yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/34);  total (35/34)

TCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAACCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCT  >  NC_000913/3948494‑3948766
                                                                                                                                         |                                                                                                                                       
tCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGcc                                                                                                                                       >  1:92896/1‑140 (MQ=255)
    tAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACAt                                                                                                                                   <  2:900616/140‑1 (MQ=255)
          aCTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATc                                                                                                                             >  1:936661/1‑140 (MQ=255)
            tGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCtt                                                                                                                           >  1:149808/1‑140 (MQ=255)
            tGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCtt                                                                                                                           <  1:225476/140‑1 (MQ=255)
              cctcctGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCg                                                                                                                         >  1:689098/1‑140 (MQ=255)
                   tGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCtctc                                                                                                                    >  1:151831/1‑140 (MQ=255)
                               cGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGaa                                                                                                        <  2:789094/140‑1 (MQ=255)
                                     aTAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCg                                                                                                  <  2:689098/140‑1 (MQ=255)
                                          gTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTAttt                                                                                             <  2:760945/140‑1 (MQ=255)
                                          gTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTAttt                                                                                             <  2:328751/140‑1 (MQ=255)
                                          gTCAGACTGATCAATAGCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCACATCTTCGCTCTCAAGTTTGCAGAATTGGCGTAttt                                                                                             <  2:718767/140‑1 (MQ=255)
                                           tCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTc                                                                                            >  1:527577/1‑140 (MQ=255)
                                             aGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCg                                                                                                  >  1:103093/1‑132 (MQ=255)
                                             aGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCg                                                                                                  <  2:103093/132‑1 (MQ=255)
                                             aGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCg                                                                                          >  1:258276/1‑140 (MQ=255)
                                             aGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCg                                                                                          >  1:258277/1‑140 (MQ=255)
                                             aGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCg                                                                                          <  1:374348/140‑1 (MQ=255)
                                             aGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCg                                                                                          <  2:258276/140‑1 (MQ=255)
                                             aGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCg                                                                                          <  2:258277/140‑1 (MQ=255)
                                               aCTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGa                                                                                        <  2:30286/140‑1 (MQ=255)
                                               aCTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGa                                                                                        >  1:181786/1‑140 (MQ=255)
                                               aCTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGa                                                                                        <  2:479168/140‑1 (MQ=255)
                                                   aTCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCg                                                                                                  >  1:217113/1‑126 (MQ=255)
                                                   aTCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCg                                                                                                  <  2:217113/126‑1 (MQ=255)
                                                   aTCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTAt                                                                                    >  1:284440/1‑140 (MQ=255)
                                                   aTCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTAt                                                                                    <  1:858894/140‑1 (MQ=255)
                                                        tatCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATccagc                                                                               <  2:181786/140‑1 (MQ=255)
                                                        tatCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATccagc                                                                               <  2:92896/140‑1 (MQ=255)
                                                           cAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATccagccag                                                                            <  2:936661/140‑1 (MQ=255)
                                                             gCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGc                                                                          >  1:768811/1‑140 (MQ=255)
                                                             gCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGc                                                                          >  1:617028/1‑140 (MQ=255)
                                                             gCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGc                                                                          >  1:311463/1‑140 (MQ=255)
                                                             gCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGc                                                                          >  2:532784/1‑140 (MQ=255)
                                                             gCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGc                                                                          >  2:187803/1‑140 (MQ=255)
                                                              cAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCa                                                                         <  1:187803/140‑1 (MQ=255)
                                                              cAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCa                                                                         <  1:458919/140‑1 (MQ=255)
                                                               aaTGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCAt                                                                        >  2:875658/1‑140 (MQ=255)
                                                                aTGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCAtt                                                                       <  2:527577/140‑1 (MQ=255)
                                                                  ggTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTACAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCa                                                                     >  2:387933/1‑140 (MQ=255)
                                                                      cGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGttt                                                                 >  1:152283/1‑140 (MQ=255)
                                                                        tGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCAt                                                                        <  2:377914/131‑1 (MQ=255)
                                                                        tGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCAt                                                                        >  1:377914/1‑131 (MQ=255)
                                                                           aaCTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCt                                                            <  1:532784/140‑1 (MQ=255)
                                                                            aCTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTg                                                           <  2:151831/140‑1 (MQ=255)
                                                                                          cTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTg                                                           <  2:193175/126‑1 (MQ=255)
                                                                                          cTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTg                                                           >  1:193175/1‑126 (MQ=255)
                                                                                          cTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGgc                                             >  2:358992/1‑140 (MQ=255)
                                                                                          cTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGgc                                             >  2:807957/1‑140 (MQ=255)
                                                                                           tGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGgcg                                            <  1:358992/140‑1 (MQ=255)
                                                                                           tGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGgcg                                            <  1:807957/140‑1 (MQ=255)
                                                                                            gCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGTACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGgcgc                                           >  1:327542/1‑140 (MQ=255)
                                                                                            gCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGc                                                       >  1:364059/1‑128 (MQ=255)
                                                                                            gCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGc                                                       <  2:364059/128‑1 (MQ=255)
                                                                                            gCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGgcgc                                           >  1:571607/1‑140 (MQ=255)
                                                                                            gCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGgcgc                                           >  1:201529/1‑140 (MQ=255)
                                                                                                           aTAACCCCAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCAtt                            <  1:399112/140‑1 (MQ=255)
                                                                                                                cccAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGTACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTAcgcg                       <  2:327542/140‑1 (MQ=255)
                                                                                                                cccAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTAcgcg                       <  2:311463/140‑1 (MQ=255)
                                                                                                                 ccAGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCgt                      >  2:395751/1‑140 (MQ=255)
                                                                                                                   aGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCgttg                    >  2:641214/1‑140 (MQ=255)
                                                                                                                   aGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCgttg                    >  2:929863/1‑140 (MQ=255)
                                                                                                                   aGATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCgttg                    >  2:455632/1‑140 (MQ=255)
                                                                                                                    gATGGATCAGTGTTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCgttgt                   >  2:836321/1‑140 (MQ=255)
                                                                                                                              tgtTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTAcg                         <  2:690154/124‑1 (MQ=255)
                                                                                                                              tgtTCCTTCCAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTAcg                         >  1:690154/1‑124 (MQ=255)
                                                                                                                                ttccttccAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGgcg       <  1:542695/140‑1 (MQ=255)
                                                                                                                                 tccttccAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGgcgc      <  2:617028/140‑1 (MQ=255)
                                                                                                                                     tccAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCgctgct  <  1:929863/140‑1 (MQ=255)
                                                                                                                                     tccAGCCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCgctgct  >  2:196641/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TCGATAGCTAACTGCCTCCTGCAAATGCTGACGTGTGATAATGTCAGACTGATCAATATCAGCAATGGTTCGTGCAACTTTCAATAACCGCTGCCAGGCACGAATCGATAACCCCAGATGGATCAGTGTTCCTTCCAACCACATCGCATCTTCGCTCTCAAGTTTGCAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCT  >  NC_000913/3948494‑3948766

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