Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 370,620 A→G H115R (CAC→CGC)  mhpB → 3‑carboxyethylcatechol 2,3‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913370,6200AG100.0% 218.3 / NA 61H115R (CAC→CGC) mhpB3‑carboxyethylcatechol 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/30);  total (31/30)

GGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCACGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACG  >  NC_000913/370484‑370755
                                                                                                                                        |                                                                                                                                       
ggAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCgg                                                                                                                                      >  2:849083/1‑140 (MQ=255)
ggAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCgg                                                                                                                                      >  2:407670/1‑140 (MQ=255)
   gCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCATCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGtt                                                                                                                                   >  2:756039/1‑140 (MQ=255)
     gACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCg                                                                                                                                 >  2:100916/1‑140 (MQ=255)
        ggCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGccc                                                                                                                              <  2:499904/140‑1 (MQ=255)
          cAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAg                                                                                                                            <  1:546175/140‑1 (MQ=255)
           aaTTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGc                                                                                                                           <  2:798544/140‑1 (MQ=255)
                gTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTg                                                                                                                      <  1:100916/140‑1 (MQ=255)
                    tttCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGt                                                                                                                  >  2:74649/1‑140 (MQ=255)
                        ggCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCt                                                                                                              >  1:658764/1‑140 (MQ=255)
                             tGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCg                                                                                                         <  1:74649/140‑1 (MQ=255)
                             tGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCg                                                                                                         <  2:449120/140‑1 (MQ=255)
                               cGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGt                                                                                                       <  1:639506/140‑1 (MQ=255)
                                 gCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTgg                                                                                                     <  2:495160/140‑1 (MQ=255)
                                 gCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTgg                                                                                                     <  1:470167/140‑1 (MQ=255)
                                        agCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGAt                                                                                              >  1:137009/1‑140 (MQ=255)
                                         gCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCggc                                                                                                       <  2:67567/130‑2 (MQ=255)
                                         gCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCggc                                                                                                       >  1:67567/1‑129 (MQ=255)
                                         gCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATa                                                                                             <  1:756039/140‑1 (MQ=255)
                                          cTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATaa                                                                                            <  1:273637/140‑1 (MQ=255)
                                            gCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTgg                                                                                                     >  2:763656/1‑129 (MQ=255)
                                            gCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTgg                                                                                                     <  1:763656/129‑1 (MQ=255)
                                            gCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAgg                                                                                          <  2:658764/140‑1 (MQ=255)
                                            gCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGAAGTTCCTGCTCGGTGGGCTGGATAAgg                                                                                          >  2:146190/1‑140 (MQ=255)
                                            gCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCA‑‑‑GGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGc                                                                                       <  1:39355/140‑1 (MQ=255)
                                               cGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGc                                                                                       <  1:803413/140‑1 (MQ=255)
                                                  gCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAg                                                                                    >  2:565871/1‑140 (MQ=255)
                                                  gCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAg                                                                                    >  1:756568/1‑140 (MQ=255)
                                                  gCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAg                                                                                    <  1:565871/140‑1 (MQ=255)
                                                          aGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCt                                                                            <  2:656149/140‑1 (MQ=255)
                                                          aGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCt                                                                            <  1:505776/140‑1 (MQ=255)
                                                              ggCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTgtgt                                                                        >  1:836300/1‑140 (MQ=255)
                                                                    ggCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTtcatca                                                                  >  2:368131/1‑140 (MQ=255)
                                                                     gCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCaa                                                                 <  2:137009/140‑1 (MQ=255)
                                                                         gtgCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGgt                                                             >  1:98863/1‑140 (MQ=255)
                                                                               aTGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGcc                                                       >  1:164004/1‑140 (MQ=255)
                                                                                tGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCa                                                      <  1:849083/140‑1 (MQ=255)
                                                                                       aTGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGAAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCt                                               <  1:146190/140‑1 (MQ=255)
                                                                                             agCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCgg                                         >  2:182305/1‑140 (MQ=255)
                                                                                              gCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGt                                        >  1:649894/1‑140 (MQ=255)
                                                                                                             ccGTTTCTTACTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGc                         >  2:561755/1‑140 (MQ=255)
                                                                                                                       cTGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTg               <  1:182305/140‑1 (MQ=255)
                                                                                                                        tGTATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGCTAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGa              <  2:756568/140‑1 (MQ=255)
                                                                                                                          tATGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAg            <  2:164004/140‑1 (MQ=255)
                                                                                                                            tGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGc                         >  1:765107/1‑125 (MQ=255)
                                                                                                                            tGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGcc          >  2:84973/1‑140 (MQ=255)
                                                                                                                            tGCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGcc          >  1:209609/1‑140 (MQ=255)
                                                                                                                            tGCAGGTGGACCGCGGGTTCGCCCACCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGGCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGc                         <  2:765107/125‑1 (MQ=255)
                                                                                                                             gCAGGTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCa         >  1:511548/1‑140 (MQ=255)
                                                                                                                                ggTGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTg      <  2:836300/140‑1 (MQ=255)
                                                                                                                                  tGGACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGa    <  1:822552/140‑1 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGAc   <  1:393269/139‑1 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGAc   >  2:393269/1‑139 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  >  1:603582/1‑140 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  >  1:385656/1‑140 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  >  1:344049/1‑140 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  >  2:618821/1‑140 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  <  2:511548/140‑1 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  >  1:690737/1‑140 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  >  2:78003/1‑140 (MQ=255)
                                                                                                                                    gACCGCGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACg  >  2:38198/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GGAGCGACGGCAATTGGTGATTTCGGCAGTGCGGCAGGAGAGCTGCCCGTGCCTGTGGAGCTGGCGGAGGCCTGTGCGCATGCCGTCATGAAGAGCGGGATCGATCTTGCCGTTTCTTACTGTATGCAGGTGGACCACGGGTTCGCCCAGCCGCTGGAGTTCCTGCTCGGTGGGCTGGATAAGGTGCCAGTTCTGCCTGTGTTCATCAACGGTGTCGCCACGCCGCTGCCCGGTTTCCAGCGTACCCGCATGTTGGGTGAAGCCATTGGACG  >  NC_000913/370484‑370755

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: ‑