Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 548,703 A→G N116S (AAC→AGC)  ylbE → DUF1116 domain‑containing protein YlbE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913548,7030AG100.0% 238.1 / NA 64N116S (AAC→AGC) ylbEDUF1116 domain‑containing protein YlbE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/28);  total (36/28)

GGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGTATTAAGCG  >  NC_000913/548568‑548837
                                                                                                                                       |                                                                                                                                       
gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTg                                                                                                                                     >  2:475383/1‑140 (MQ=255)
 ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGa                                                                                                                                    <  2:273346/140‑1 (MQ=255)
              gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCg                                                                                                                       >  2:557522/1‑140 (MQ=255)
                 cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGc                                                                                                                    <  1:330559/140‑1 (MQ=255)
                  aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCa                                                                                                                           >  1:615885/1‑132 (MQ=255)
                  aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCa                                                                                                                           <  2:615885/132‑1 (MQ=255)
                         ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCGGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCa                                                                                                            >  2:886632/1‑140 (MQ=255)
                          gCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCaa                                                                                                           <  1:557522/140‑1 (MQ=255)
                           cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAAc                                                                                                          <  2:982681/140‑1 (MQ=255)
                              tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTc                                                                                                       >  2:651429/1‑140 (MQ=255)
                                 tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAAc                                                                                                    >  1:6863/1‑140 (MQ=255)
                                  ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACg                                                                                                   <  1:405366/140‑1 (MQ=255)
                                    aGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGGGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGaa                                                                                                 >  1:710121/1‑140 (MQ=255)
                                    aGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGaa                                                                                                 >  1:1071536/1‑140 (MQ=255)
                                        gggggAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAAc                                                                                                          <  1:982779/127‑1 (MQ=255)
                                        gggggAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAAc                                                                                                          >  2:982779/1‑127 (MQ=255)
                                        gggggAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTa                                                                                             >  1:678037/1‑140 (MQ=255)
                                        gggggAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTa                                                                                             >  2:1001966/1‑140 (MQ=255)
                                        gggggAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTa                                                                                             >  1:441715/1‑140 (MQ=255)
                                         ggggAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTAt                                                                                            >  2:527205/1‑140 (MQ=255)
                                                    cTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGa                                                                                 >  1:319916/1‑140 (MQ=255)
                                                       cATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGc                                                                              >  1:650052/1‑140 (MQ=255)
                                                              tGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTgg                                                                      >  2:973061/1‑140 (MQ=255)
                                                              tGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTgg                                                                      >  2:12208/1‑140 (MQ=255)
                                                                tCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTG‑‑ACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCt                                                                  <  2:39863/140‑1 (MQ=255)
                                                                 cACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGgcgc                                                                   >  1:983490/1‑140 (MQ=255)
                                                                      tgtgAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACg                                                              >  1:210794/1‑140 (MQ=255)
                                                                      tgtgAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACg                                                              >  1:759689/1‑140 (MQ=255)
                                                                        tgAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGc                                                            >  1:689835/1‑140 (MQ=255)
                                                                          aaTGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTTGGCGCTTACGGCg                                                           >  1:416189/1‑140 (MQ=255)
                                                                          aaTGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGa                                                          >  1:869649/1‑140 (MQ=255)
                                                                           aTGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGgcg                                                                    <  1:1072587/129‑1 (MQ=255)
                                                                           aTGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGgcg                                                                    >  2:1072587/1‑129 (MQ=255)
                                                                           aTGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGaa                                                         <  2:755127/140‑1 (MQ=255)
                                                                           aTGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGaa                                                         <  2:678037/140‑1 (MQ=255)
                                                                              ccGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGAt                                                      >  2:163884/1‑140 (MQ=255)
                                                                               cGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATg                                                     >  2:376339/1‑140 (MQ=255)
                                                                                 tCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTc                                                   >  1:692095/1‑140 (MQ=255)
                                                                                 tCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTc                                                   >  1:1103296/1‑140 (MQ=255)
                                                                                 tCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTc                                                   <  1:651429/140‑1 (MQ=255)
                                                                                 tCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTc                                                   >  2:75280/1‑140 (MQ=255)
                                                                                        aaTGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTc                                            <  2:1030484/140‑1 (MQ=255)
                                                                                              cGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATc                                      <  2:210794/140‑1 (MQ=255)
                                                                                                gTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGc                                    <  2:869649/140‑1 (MQ=255)
                                                                                                      cTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATg                              <  2:675368/140‑1 (MQ=255)
                                                                                                        tCTGCCAGTATGCCGATGCTGGTGGTTGGGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATgcg                            <  2:710121/140‑1 (MQ=255)
                                                                                                        tCTGCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATgcg                            <  1:1001966/140‑1 (MQ=255)
                                                                                                           gCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGa                         >  1:32410/1‑140 (MQ=255)
                                                                                                           gCCAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGa                         >  1:207492/1‑140 (MQ=255)
                                                                                                            ccAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAt                        <  2:983490/140‑1 (MQ=255)
                                                                                                            ccAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAt                        <  2:32410/140‑1 (MQ=255)
                                                                                                             cAGTATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAtg                       <  1:653823/140‑1 (MQ=255)
                                                                                                                tATGCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATgtgt                    >  2:572806/1‑140 (MQ=255)
                                                                                                                   gCCGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTaa                 <  1:12208/140‑1 (MQ=255)
                                                                                                                     cGATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATg               <  2:759689/140‑1 (MQ=255)
                                                                                                                      gATGCTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGc              <  1:163884/140‑1 (MQ=255)
                                                                                                                          cTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGt          >  1:1029701/1‑140 (MQ=255)
                                                                                                                          cTGGTGGTTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGAGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGt          >  1:652041/1‑140 (MQ=255)
                                                                                                                             gtggtTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTaa                 <  2:569023/130‑1 (MQ=255)
                                                                                                                             gtggtTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTaa                 >  1:569023/1‑130 (MQ=255)
                                                                                                                             gtggtTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGTAtt       <  2:6863/140‑1 (MQ=255)
                                                                                                                             gtggtTGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGTAtt       >  2:919328/1‑140 (MQ=255)
                                                                                                                                  tGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGTATTAAgcg  <  2:1071536/140‑1 (MQ=255)
                                                                                                                                  tGAGAGCGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGTATTAAgcg  <  1:49451/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGATGTGTTAATGCCAGTATTAAGCG  >  NC_000913/548568‑548837

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